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DIGITAL adaptive CONTROLLERS 
FOR VTOL VEHICLES- -VOLUME II 
Gary L, Hartmann, Gunter Stein, and Stephen G. Pratt 

Honeywell Inc, 

Systems and Research Center 
Minneapolis, Minnesota 

SECTION 1 
INTRODUCTION 

This volume documents the VALT adaptive software developed for NASA 
Langley Research Center under Contract NASl- 14921, Two self-adaptive 
algorithms were evaluated on NASA's VALT (VTOL Approach and Landing 

Technology) simulation. One is based on an implicit model reference 

1 2 
design, the other on an explicit parameter estimation technique. The 

latter was recommended for flight test in the VALT research helicopter. 

The remainder of this volume presents the organization of the software, 

user options, and a nominal set of input data. This volume also contains a 

flow chart and program listing of each algorithm. For a discussion of the 

design theory of these algorithms, the reader is referred to Volume I. 


^Referred to as the Model Reference (MR) algorithm. 

2 

Referred to as the Maximum Likelihood Estimation (MLE) algorithm. 



SECTION 2 

OVERVIEW OF VALT ADAPTIVE SOFTWARE 

The VALT adaptive software consists of separate longitudinal and 
lateral-directional controllers plus an adaptive gain adjustment algorithm 
for the longitudinal controller. The two controllers were designed to accept 
independent commands from a guidance algorithm. The guidance commands 
are the the three attitudes plus vertical velocity. The primary objective of 
the control law is to provide good command- following With minimal cross- 
axis response. The control laws provides differential collective and collec- 
tive servocommands for the longitudinal axis. Cyclic and differential cyclic 
commands are provided for lateral-directional control. 

The Honeywell- supplied software consists of separate subroutines writ- 
ten in standard FORTRAN for the CDC computer facility at NASA LRC. All 
data required by the real-time subroutines are stored in labeled arrays. 

The subroutines are listed in Table 1 and the arrays (common blocks) in 
Table 2. 

The interface between the adaptive control software and NASA's VALT 
simulation on the CDC CYBER computer is contained in subroutine CONTRL3. 
A flow chart of this subroutine is shown in Figure 1. The data interface 
between the NASA simulation and the Honeywell- supplied subroutines is the 
VARDAT array. 

Two different adaptive algorithms were evaluated for the gain adjust- 
ment function. The dashed line in Figure 1 shows the call statement to the 
MR algorithm. The algorithm that was selected adjusts the gain based on 
explicit parameter estimates from a parallel- channel Maximum- Likelihood 
Estimator (PCMLE). Consequently, the dashed line path does not exist in 
the current CONTRL3 subroutine. 
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TABLE l.-LIST OF SUBROUTINES 


Subroutine 


Function 


CONTRL3 


Real-time interface between VALT adaptive 
software and simulation on CDC CYBER 
computer. 


HCONl 


Pitch-axis control law. 


HCON2 


Lateral-directional control law. 


GAINI 


Gain schedule for pitch axis. 


GAIN2 


Gain table for pitch axis. 


PCMLE 

FH 

TSIG 

FILT 

SENS 

ACUM 


Real-time subroutines for MLE algorithm. 


NRTIC 

MODEL 

DISC 

FHIC 

DIAK 

CAL 

MP 


Initialization and filter design software for 
the MLE algorithm. 


INPT 


RFMOD 

INTEGRAL 

Real-time MR algorithm. 

RMIC 

Initialization for MR algorithm. 
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table 2. -COMMON BLOCKS USED BY VALT 
ADAPTIVE SOFTWARE 


Common Block 

Function 

Cl 

Stores four integrator values for control laws. 

VARDAT 

Contains UX and LX arrays. Used for all data 
tranter and mode logic. 

DAT 

SENSP 

IPIC 

MEAS 

1 Data storage used by MLE algorithm. 

SUBR 

Data storage used by MLE initialization. 

SUBRH 

Stores frequency and damping constants for 
highpass filters. 

UTEST 

Stores test signal parameters. 

DEVICE 

Stores logical assignment of card reader and 
line printer. 


MRINT 

MRDAT 


Data storage for MR algorithm. 




































, SECTION 3 

THE VALT ADAPTIVE CONTROL LAWS 
Description 

Separiate longitudinal and lateral- directional controllers are implemented 
in subroutines HCONl and HCON2. Each of the controllers forms servocom- 
mands from a combination of guidance commands and measured responses. 
Program listings are continued in Appendix A, 

The longitudinal controller determines the differential collective com- 
mand (UCl) and the collective command (UC2) from the 6 and V measure- 

z 

ments. Proportional plus integral shaping is applied to each command. The 
two command integrators are denoted XCl and XC2 in the Cl array. 

The functional block diagram of the longitudinal controller is shown in 
Figure 2, which also shows the various gains stored in the TJX array. The 
computation of these gain elements is discussed later. The variables used 
in HCONl are defined in Table 3, 

The lateral- directional controller determines the cyclic (UC3) and dif- 
ferential cyclic commands (UC4). The controller assumes guidance com- 
mands in the form of 0 and \}/ attitude commands. The controller uses pro- 
portional plus integral control for each command and uses measurements of 
p, r, 0, and p as feedback. The two command integrators are XC3 and XC4 
in the Cl array. 

The functional block diagram of the lateral controller is shown in Figure 
3, and the symbols are defined in Table 4, 
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Gain matrix elements 

K = UX(51) 0.8UX(52) 

UX(53) 0.8UX(54) 

Ki = UX(44) 0.8UX(45) 

' UX(49) 0.8UX(50) 

Kxp = 'UX(4D UX(42) UX(43)1 

^ [uX(46) UX(47) UX(48)J 

Figure 2. -B\inctional block diagram for pitch-axis control. 


7 




\ , 
i. 

’r ' '■ 

r 

j::. vr" 

1 ^ TABLE 3. -HCONl VARIABLES 


Variable 

Definition 

UX(1) 

Total command, ft/s 

2 

Total 0 command, rad 

11 

V^, measured vertical velocity, ft/s 

12 

q, measured pitch rate, rad/s 

13 

6, measured pitch attitude, rad 

19 

At, sample rate, sec 

41 

• 

1 

• 

• 

) 10 gain values used by pitch -axis controller. 

50 

) 

XCl 

Integral for differential collective command. 

XC2 

Integral for collective command. 



Subscript c denotes commanded value, 

UC3 = Cyclic command (in.) 

UC4 = Differential cyclic command (in.) 

Gain matrices 

Figure 3. - Functional block diagram for lateral-directional 
control. 
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TABLE 4. -HCON2 VARIABLES 


Variable 

Definition 

UX(3) 

0 command, rad 

4 

xj/ command, rad 

19 

DT, sample time, s 

37 

p, measured roll rate, rad/s 

38 

r, measured yaw rate, rad/s 

39 

0, measured roll attitude, rad 

40 

rj/, measured yaw attitude, rad 


Initialization 


The four integrators in the Cl array must be initialized using trim values 
to prevent transients when the controllers are engaged. The appropriate 
initial values for XCl, XC2, XC3, and XC4 are determined with the follow- 
ing relations from Figures 2 and 3: 


rv 1 fK “1 

Zt PP + 

^tj L J 

'XCll 

-xcd 

[M rvi , 

^ J 

■XC31 

J?C4j 





where the subscript "t” denotes a trim value existing prior to engaging the 
control law and the rates are assumed zero. 
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b 

:n i ; i 


The variable UX(19) must be initialized with the appropriate sample 
time. The remaining program initialization is specific to either the MR or 
the MLE adaptive algorithm and will be discussed in later sections. 

Real-Time Operations 

The VA LT adaptive system is structured to require a call to subroutine 
CONTRL3 at each sample time. This subroutine controls all the signals 
sent to the various subroutines and contains all the logic for exercising the 
various user options discussed below. A program listing is given in Appendix 
A . The listing for the gain computation is given in Appendix B. 

User Inputs and Options 

The interface between the Honeywell adaptive software resides in sub- 
routine CONTRL3. Data are exchanged via the UX array. Various options 
for using the software have been provided by a set of logical variables (LX) 
and five panel switches, FSS(9) through FSS(13), located on the CDC console. 
The various options are defined below: 

• Panel Switches: 

9 - Engages the pitch-axis controller and continually up- 
dates PCMLE high pass filters. 

10 - Positive command. 

11 - Negative command. 

12 - Engages the PCMLE estimator [requires FSS(9) to 

be on], 

13 - Engages the lateral controller. 
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LX A rray ; 


1 - T fixes the PCMLE algorithm on the channel deter- 

mined by starting value of JS. 

- F uses likelihood select logic to determine channel 
used. 

2 - T sets NOEST logic which skips parameter updates 

in PCMLE. 

- F nominal operation. 

3 - T replaces measured v with v estimate from PCMLE 

for gain schedule. 

4 - T uses table lookup in GAIN subroutine. 

- F uses schedule in GAIN2 subroutine. 

5 - T = PRINT prints out selected PCMLE variables. 

This can only be used in a batch mode. 

LX; 

6 - T selects commands from automatic guidance. 

7 - Not used. 

8 - Controls printout in NRTIC subroutine. 

9 - T engages the estimation computation set by switch 12. 

10 - T causes pitch controller to use a fixed set of gains 

from subroutine GAIN2. 

11 - T applies step command via FSS(IO) and FSS(ll) to 

Vj, input. 

12 - T applies step command via FSS(IO) and FSS(ll) to 

6 input . 



13 - T applies step command via FSS(IO) and FSS(ll) to 

0 input, 

14 - T applies step command via FSS(IO) and FSS(ll) to 

input. 

As previously discussed, the logic to select the MR algorithm was 
deleted following selection of the MLE algorithm for flight evaluation. For 
completeness, the MR algorithm is documented in a later section of this 
report. 

Outputs From VALT Adaptive Software 

The primary outputs from the software are the four servocommands 
(UCl, UC2, UC3, UC4) previously discussed. For purposes of checkout, 
eight variables are output on a strip chart for real-time plots. Other selected 
variables are available in the UX array for monitoring from the CDC console. 

The variables plotted are variables from the MLE algorithm and are 
defined in Section 4 (Table 12). 
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SECTION 4 

THE PA RA LLE L CHA NNE L MA XIM UM LIKE LIHOOD 
estimation ALGORITHM 

Description 

This algorithm uses measurements of vertical velocity, pitch rate, and 
pitch attitude together with the collective and differential collective positions 
to estimate the dominant parameters of the VALT helicopter. The estimation 
is based on a maximum -likelihood method using parallel processing of the 
data by Kalman filters. Each filter is a four- state representation of the 
longitudinal dynamics. 

The explicit parameter estimates are used to algebraically adjust the 
feedback and feedforward gains of the pitch-axis control law. 

The estimation software is organized as shown in Figure 4, It consists 
of a non-real-time segment (definition and initialization of Kalman filter 
channels) and the real-time parameter identification. The initialization 
program permits different designs by providing the following flexibility: 

• Variable number and location of Kalman filter channels 

• Variable sample rates 

• Variable number of identified parameters (up to four) 

• Significance tests to monitor data validity 

The computations have been divided among a number of subroutines. 

The subroutines are listed in Table 5 together with their use in either 
initialization or real time. For fast execution, most of the data exchanged 
between subroutines is done via labeled common blocks. Seven blocks are 
used in the estimation software. They are defined in Table 6. 
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Initialization 


Real-time 



Figure 4. -PCMLE software structure. 
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table 5. -PCMLE SUBROUTINES 


Subroutine 

Function 

Non- real-time: 


NRTIC 

Main executive routine for non- real-time opera- 
tion. Reads data to define number and location 
of channels, number of parameters estimated, 
sample rate, etc. Performs all initialization with 
calls to other subroutines. 

MODEL 

Defines the system matrices and sensitivities 
for a discrete four- state model for each channel. 

DISC 

Computes disctete model from a continuous- 
time model. 

FHIC 

Computes filter coefficients and initializes high- 
pass filter applied to measurements. 

DIAK 

Solves for the Kalman filter gains for a discrete 

CAL 

system. It uses an iterative procedure. 

INPT 

Reads card input for starting value of filter gains. 

MP 

Prints out matrix. 

Real-time: 


PCMLE 

Main executive routine for parallel-channel MLE 
real-time computations. 

FH 

Highpass filter applied to measurements. 

TSIG 

Produces test signal and two random numbers for 
simulated sensor noise. 

FILT 

Performs a fourth-order Kalman filter computa- 
tion for each pitch channel. 

SENS 

Performs a sensitivity filter for a given param- 
eter. 

ACUM 

Accumulates likelihood gradients and approximate 
second partials for a Newton- Raphson parameter 
update. 
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TABLE 6. -li^BELED COMMON BLOCKS FOR DATA TRANSFER 


Common 

block 

Function 

yard AT 

Contains UX and LX arrays. Used to input user data and 
output selected real-time parameters. Defined in Tables 
8 and 9, 

DAT 

Main arrays for storing channel model parameters. De- 
fined in MODEL and used in NRTIC, PCMLE, SENS, ACUM, 
and FILT subroutines. The X, F, and D arrays of this 
common block are defined in Tables 3, 4, and 5, 

IPIC 

Contains all the parameter initialization defined by NRTIC 
and transferred to the real-time routine PCMLE. 

MEAS 

Stored current value of high-passed measurements. Initial- 
ized in NRTIC, defined in PCMLE and used in FILT, and 
SENS subroutines. 

SUBR 

Transmits model data to DIAK and CAL for solving a 
Ricatti equation to determine Kalman filter gains. Used in 
NRTIC, PCMLE, MODEL, DiAK, and CA L subroutines. 

SUBRH 

Stores highpass filter constants. Used in NRTIC, PCMLE, 
MODEL, FHIC, and FH subroutines. 

SENSP 

Stores the sensitivity of each channel model variable with 
respect to each of the four parameters estimated. See 
Table 10 for details. Used by NRTIC, PCMLE, MODEL, 
SENS, and ACUM subroutines. 


Definitions of the arrays used in the MLE software are given in Tables 
7, 8, 9, and 10, 


Initialization 


The initialization of PCMLE is performed out of realtime with a call 
to subroutine NRTIC. This subroutine reads a data file containing all the 
common blocks used by the PCMLE software. This data file is generated 
by a separate batch program which is discussed later in this section. 
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table 7. -DEFINITION OF X-ARRAY IN DAT COMMON BLOCK 


X (IC, I) 

IC = 1, 5 for channel index 
1 = 1, 19 for variable 


Index (I) 


Variable 

1 

Vx 


2 


' Current state at time tj^ 

3 



4 

^k 


5 

^x 

^k+1 


6 

^k+L 

t Predicted state at time tj^^^ 

7 

^k+1 


8 

\+l 


9 



10 

"q 

\ Current filter residual 

11 



12 

v*^R ■ 

1 

V 

13 

v*R " 

+ Indet R 


f ^ 1 

14 

sum V 

R 

15 

sum v^R + Indet R 

16 

— 


17 

— 


18 

— 

\ 

19 

— 





table 8, -DEFINITION OF F-ARRAY IN DAT COMMON BLOCK 


F (IC, I) 

IC = 1, 5 filter (channel) index 
1=1, 22 for parameter index 


Index 


KF(1. 1) 
POi’(2, 1) 
KF(3, 1) 
KF(4, 1) 
KF(1,2) 
KF(2,2) 
KF(3,2) 
KF(4, 2) 
KF(1,3) 


10 

KF(2,3) 

11 

KF(3,3) 

12 

KF(4,3) 

13 

Ln Det 

14 

R-l(l, 1) 

15 

R-l(l, 2) 

16 

R-l(l,3) 

17 

R'^(2, 1) 

18 

R-l(2,2) 

19 

R-l(2,3) 

20 

R“l(3, 1) 

21 

R-l(3,2) 

22 

R-l(3,3) 


Variable 


Kalman filter gains 


Inverse of residual co- 
variance matrix 
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TABLE 9. -DEFINITION OF D-ARRAY IN DAT COMMON BLOCK 


D (IC, I) 

IC = 1, 5 for channel index 
1 = 1, 22 for parameter 
For the model; 


X 


n+1 


Ax„ + Bu 
n n 


Index (I) 

Variable 

1 

A(l, 1) 1 


2 

A(1,2) 


3 

A(1,3) 


4 

A(l. 4) 


5 

A (2, 1) 


6 

A (2, 2) 


7 

A (2, 3) 


8 

A (2, 4) 

> Non- zero elements of discrete 

9 

A(3, 1) 

A matrix 

10 

A (3. 2) 


11 

A (3, 3) 


12 

A (3, 4) 


13 

A (4, 1) 


14 

1 A (4, 2) 


15 

A (4, 3) 


16 

A (4, 4) . 


17 

B(l, 1) “ 


18 

B(l, 2) 


19 

B(2, 1) j 

Non- zero elements of discrete 


B(2.2) 

B matrix 


B(3, 1) 


22 

B(3, 2) _ 
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table 10. -DEFINITION OF XS AND DS ARRAY IN SENS? 

COMMON BLOCK 


Array 

Definition 

XS 

(IP. I) 

IP = 1, 4 parameter index 

I =1,19 variable index for sensitivity filter 
quantities 

Same list for las Table 3. 

DS 

(IC, IP, I) 

IC = 1, 5 for channel index 
IP = 1, 4 for parameter index 
I =1, 22 for sensitivity quantities 

in general DS <IC, IP, I) = ^ ,p 

where I is as given in Table 5. 


Real-Time Operation 

The real-time computations are executed with a callto subroutine PCMLE. 
PCMLE is called once per sample time. The sampled values of pitch rate, 
pitch attitude, down velocity, and collective and differential collective posi- 
tion are highpassed, and residuals and likelihood functions are computed for 
each Kalman filter (fixed in parameter space). Gradients and second partials 
are accumulated for a Newton- Raph son parameter update. The remaining 
real-time operations ai^e spread over five subcycles as shown in Figure 4. 
During real-time operation, selected variables can be monitored using the 
UX array. 

The real-time data that are input to the algorithm at each sample time 
are defined in Table 11. The data are transferred from the VALT simulation 
via the UX array. 
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table 11. -REAL-TIME INPUT VARIABLE ASSIGNMENT 


User 

variables 

Description 

M neumonic 
expression 

Units 

UX(ll) 

"—I 

Down velocity 

ZDOTE 


12 

Pitch rate 

QI 


13 

Pitch attitude 

THEI 


14 

Differential collective 

DC(1) 

in. 

15 

Collective 

DC (2) 

in. 

16 

Forward velocity^ 

XDOTE 

ft/s 


3 . 

Not used by controller. Used for evaluation only. 


Selected output variables are available for monitoring performance. 
They are contained in the UX array and are defined in Table 12 . 

A flow chart of the real-time computations is contained in Appendix C. 
A complete program listing appears in Appendix D. 

User Inputs and Options 

Five logical variables (LX) are used in the software as switches to 
control various PCMLE user options. These variables are defined in Table 
13. The nominal settings are LX(9) = True, all others false. 

The other modes are used for checkout and debugging purposes. LX(1) 
= T eliminates the channel switching logic so the fit of any of the Kalman 
filters can be checked against the input data, LX(2) = T skips the param- 
eter estimation steps. LX(5) and LX(8) control line printer output, LX(9) 
is the logical used to engage the PCMLE algorithm for normal operation. 
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TABLE 12, -PCMLE OUTPUT VARIABLES 


User 

variable 

Description 

Mneumonic 

expression 

Units 

UX(20) 

Channel indicator 

0. 25JS-. 5 

■i 

21 

Vg estimate (highpassed) 

X(JS, 2)/20 


22 

q estimate (highpassed) 

X(JS, 3)=:<5. 73 


23 

Scaled- likelihood function, filter 1 

TJ(l)/400 


24 

Sea led -likelihood function, filter 2 

TJ(2)/400 


25 

Sea led -likelihood function, filter 3 

TJ(3)/400 


26 

V estimate 

ZPl 

- 

27 

V accuracy 

SQRT(C1/DET) 

- 

28 

bg-j^ estimate 

ZP2 


29 

b 32 estimate 

ZP3 


30 

agg estimate 

ZP4 


31 

Vx estimate 

X(JS, 1) 

ft/s 

32 

6 estimate 

X(JS,4) 

rad 

33 

Residual ratio 

SIGSQ 

- 


TABLE 13. -USER LOGICA L ASSIGNMENTS^ 


User 

logicals 

Function 

Mneumonic 

expression 

LX(1) 

Disables automatic channel 
changes. 

NOCHC 

2 

Disables Newton- Raph son 
parameter step. 

NOEST 

5 

Enables printout from PCMLE 
(use in batch mode only). 

PRINT 

8 

Disables initialization print- 
out in subroutine model. 


9 

Enables PCMLE algorithm. 

MLE 


^Default values; LX(I) = , F, 
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Generation of PCMLE Data File 


The parameters used by the real-time PCMLE software are computed 
off-line and stored in a data file. This program reads the input data deck 
and defines the specified number of channels at the specified parameter 
values. Each channel is a four- state Kalman filter for the longitudinal dyna- 
mics, The states are forward and down velocity, pitch rate, and the pitch 
attitude. The measurements are down velocity, pitch rate, and pitch atti- 
tude, The initialization program discretizes the Kalman filter model for the 
specified sample rate and solves a Ricatti equation for the steady- state 
filter gains. Parameter sensitivities are computed for each channel using a 
numerical differencing method by individually perturbing each of four 
parameters to be estimated, 

A flow chart of this program is contained in Appendix E, A complete 
program listing is contained in Appendix F, 

i 

The various input parameters are defined in Table 14, A nominal set 
of input data and the resulting initialization file is contained in Appendix G, 

If a stable Kalman filter is not found for the given set of input data, a 
STOP 41 or STOP 31 will be encountered. These are defined in Table 15, 
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TABLE 14. -PCMLE INPUT PARAMETERS 


Input 

Description 

NX 

Number of states in Kalman filters 

NR 

Number of measurements in Kalman filters 

NN 

Number of noise sources in Kalman filters 

NTER 

Number of terms used in series to compute discrete model 

NRM 

Maximum number of measurements 

EE 

Convergence criteria in DIAK for Ricatti solution 

ITER 

Maximum number of iterations in DIAK for Ricatti solution 

DTO 

Nominal sample time (automatically increased if samples 
skipped) 

DT PRINT 

Print interval when LX(5) is true (available in batch mode 
only) 

NC 

Number of channels 

NPO 

Number of parameters estimated 

JSO 

Starting channel location 

WUTO 

Natural frequency of shaping filter on random test signal: 

DUTO 

- Damping of shaping filter on random test signal 

SIGUTO 

- Gain of shaping filter on random test signal 

UTIO 

- Initial state of shaping filter on random test signal 

UT20 

- Initial state of shaping filter on random test signal 


(Above state initialization is useful for generating deter- 
ministic square waves or sine waves with TSIG subroutine. ) 

UTMAX 

Limits size of filtered test signal 

TAUPO 

Time constant of likelihood accumulation filter 

TAUP2 

Time constant of lowpass filter on L, VL, v2l inputs 

WHP 

Frequency of second-order highpass on measurements 

DHP 

Damping of second-order highpass on measurements 

SIGVZ 

RMS noise statistic assumed, for down-velocity measurement 

SIGQ 

RMS noise statistic assumed for pitch- rate measurement 

SIGTH 

RMS noise statistic assumed for pitch-attitude measurement 

SIGXG 

RMS statistic assumed for forward wind gusts 
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TABLE 14. -Concluded 


Input 

Description 

SIGZG 

RMS statistic assumed for vertical wind gusts 

RTJCO 

THRTJCO 

> Threshold parameters controlling channel switch 

RTJZO 

THRTJZO 

/ 

^Threshold parameters controlling ZIMIN selection 

RTJSO 

/ 

Ratio test parameter for significance of likelihood function 

ZPl MAX 

Maximum value of parameter 1 estimate 

ZP2 MAX 

Maximum value of parameter 2 estimate 

ZP3 MAX 

Maximum value of parameter 3 estimate 

ZP4 MAX 

Maximum value of parameter 4 estimate 

ZPl MIN 

Minimum value of parameter 1 estimate 

ZP2 MIN 

Minimum value of parameter 2 estimate 

ZP3 MIN 

Minimum value of parameter 3 estimate 

ZP4 MIN 

Minimum value of parameter 4 estimate 

GSQLO 

Elements of '^^L^ matrix 

ZP 

Matrix defining location of channels (four parameters per 
channel times five channels maximum = 20 parameters 


table 15. -PCMLE STOP CODES 


STOP code 

Condition 

41 

Initialization 

31 

Ricatti equation not converging in subroutine 
CAL during initialization (unstable model). 

Inverse does not exist in subroutine DIAK for 
computing filter gains. Check Data Deck. 
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SECTION 5 

THE MODEL REFERENCE ADAPTIVE ALGORITHM 

Description 

The model reference algorithm adjusts the pitch-axis gains to reduce the 

error vector between the model response for V , q and 9 . The resulting 

z 

gain values are stored inthe UX array for use with subroutine HCONl. All 
the data required by this algorithm are stored in the MRDAT common block. 

Initialization 

This algorithm is initialized with a call to subroutine RMIC, This call 
defines the A and B matrices fora three-state continuous-time model and 
the parameters required by the integration subroutine INTGRAL which are 
stored in common block MRINT. The INTGRAL subroutine Contains three 
integration routines. The first-order Euler method is used to integrate the 
thirteenth -order system composed of three model states plus 10 gains. The 
parameters used in the MRDAT array are defined in Table 16. 

The X-array stores the current state. The components are: 


X(l) = V model 
z 

X(2) = q model 
X(3) =e model 


X(4) 

through 

X(13) 


10 gain values 
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TABLE 16. -DEFINITION OF MODEL REFERENCE VARIABLES 


Variable 

Definition 

X 

Current state 

DX 

Current derivative 

THRl 

Threshold on V error in gain adjustment 
z 

THR2 

Threshold on q error in gain adjustment 

THR3 

Threshold on 6 error in gain adjustment 

Gll 



G12 


- 1 

G21 


Elements of B in gain adjustment 

G22 



T1 

Scaling parameter for the P-matrix 

T2 - Til 

Set to 1 or 0 depending on which of the 
gains are to be adjusted 

EVZ 

V error (ft/s) 
Z 

EQ 

q error (rad /s) 

ETH 

6 error (rad) 

FAC 

) 


FACl 

1 

* Error quantities used for monitoring only 

Pll 



P22 


Elements of negative definite matrix satis- 

P23 


fying Liapunov equation 

P33 

1 




The DX-array stores the current derivative. The equation for the 
third- order model is: 



rvi 

r* 


rvi 

rr 


'u(ir 

d_ 

dt 

z 

q 

= [A] 

Z 

q 

[B] 


_0 . 


-0 . 


.U(2). 


where 

U(l) = V command 
z 

U(2) =0 command 



-2 

0 

0" 


’2 

0" 

A = 

0 

-4 

-7 

B = 

0 

7 


_ 0 

1 

0. 


.0 

0. 


When the controller is added to this model, the closed-loop model has un- 
coupled and 0 responses with the following transfer functions: 

_ 2(s + 1) 

V "" (s + 2) (s + 1) 

"=c 

0 ^ 7(s + 0. 8) 

®c (s + 0. 8) (s^ + 4s + 7) 

The gains consist of a feedback matrix, K^, plus a command input, K^. 
Thus, the commands are determined as: 


■'ucr 

_UC2. 

. N 

rvi 

z 

q 


_9c - »> 



-0 . 
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These matrices are initialized as: 

^ fb -7.6 -271 

^ [0.3 0 oj 

K = fO 

^ [-^.23 0 J 

Real-Time Operation 

The gain adjustment is done with a call to subroutine RFMOD, This 
routine defines the command inputs, measurements, error vector, and 
time derivatives for the gains. The actual integration of the model and the 
gain element is done with subroutine INTGRAL. The integration is performed 
for the sample-time defined for the controller. 

The derivatives of the gain matrices are defined by the following rela- 
tions: 

= -B Pe X 

dtK„ = -B “ 

where 

e = error vector of V , q, 0 errors 

Z 

X = measured V , q, 0 states 
z 

A flow chart of RFMOD is shown in Figure 5. A program listing of 
RMIC, RFMOD, and INTGRAL is contained in Appendix H. 
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Input commands 

AV, A0. 



= model states 
m 


Figure 5. - Flow chart for subroutine REMOD. 



APPENDIX A 

PROGRAM LISTING OF CONTROLLERS 


This appendix consists of the program listings of subroutines CONTRL3 
(Fig. 6), HCONl (Fig. 7). and HON2 (Fig. 8). 


SUBROUTINE C0NTRL3 


73/7R QPT-1 


fTN <>*64<rS)2 


10 


lb 


20 


25 


30 


,35 


RO 


55 


SUBROUTINE CQNTRL3 


THIS SU3ROUTIN6 COMPUTES CONTROLLER COMMANDS* ALL VARIABLES 
ENDING WITH I AAfc INERTIAL MBASUREHENTS. THOSE ENDING WITH 
6 ARE FROM THE GUIDANCE COUPLER.THEN AND PHIN ARE FROM STABOER* 
UNITS OF DCC ARE IN INCHES* CT IS THE ITERATION PERIOD FUR THE 
CONTROL MODULE* ALL PAST VALUES ARE INITIALIZED IN ICONTR* 


COMMON/ EST/LATt# LONE #ZEfXDOTE#YOOTE#ZOOTE*SSVEL WE/ SSTMWE.XOTAOE# 

A yDTADE»ZOTADE#XE> YE 

C0HM0N/AR/DAH( 4)/0LH( <»)#DA1H( A)#DL1H( <»)#0A2H( <i)>DL2H( 

1 DKLH( 4)/DRLlH( A } / DCMIN( ) / DCM AX ( A )/ORLHaX U) # WAH( 4 ) ^ 

2 ZAH(A)/WRH(A)/ZRH( A)#C0NI(6/ A)/C0N2(6/ A)/DTRL>FAR/ 

3 DHH( A)/DH1H( A)/DH2H( A)#VAR2( A)»VAR12( A)>VAR22( A)/ 

A DC(A)/OCN(A)/GAINH(A) 

COHMQN/CQNST/NC»G# RE# RE2>PP> SCHULER# ER ATE# PPldO^PP 16 0T/PP2 
COMHON/CONTROL/DCC (A)#DC1C2#W1G#0PL(A)#PSIYC#PSIY6 
COMKON/FREO/RGT# ALT#RIGT#A0T»GST#GCT#GPT#GRT#CT#PT#HRT#AAT#ENVI# 

1 ENV2#ENV3#ENVA# ENV5#FNV6# EN V7# ENV6# ENV9#XLNT# AN AT# 

2 FKT#aTCT#PTCT/3N0MT#GSLT#DT#TIMEH 
COMMON/GUIDE/UG# VG#WG#XDOTG#YDOTG#ZDOTG#UAG#VAG#WAG#ZRG#PSIHRG 
COMMON/HELCOP/THEH#PHIH#PSIH#DUOOTH#VDOTH#OWDOTH#UH# VH#WH#POOTH# 

1QDOTH#ROOTH#PH#QH#RH#PSIIHR#ALPHAH#0ETAH#FXH#FYH# 
2FZH#PPH#QQH#RRH#X0TAVH 

COMMON/INERT/TH£I#PHII#PSII#UDOTI#VOOTI#WDOTI#PI#01#RI#POI#ODI# 

1 RDI#SIGBRII#$IGRRII#SIGFRII#SIGBRI#SI6RRI#SIGSFRI# 

2 SlG6LAI/SIGRLAI#RRir(3)#RRI(3)#RLAI(3)#BRII(3)# 

3 SFRII(3)#BRI(3)#SFR1(3)#BLA1(3) 
COMMON/NOH/THEN#PHIN#UN#VN#WN#UAN#VAN#WAN 
C0MMaN/SPRNT/DUM0UTU5)#IDUM 

COMMON/TRAJECT/XDDGSN# xooagsn#xdotgsn#xdtagsn# zdogsn#zootgsn# 

A Z0T16SN#Z6SN#X0TNGSN# VNGSN#YOT£GSN#VEG$N#LATGSN# 

8 L0N6SN#XGSN# YGSN#PGSN#QGSN#RGSN#R1GSN#THEGSN# 

C PHIGSN# PSI6SN/PSIHRN#T6$N#XDD6$C#XOOAGSC#XOaTGSC# 

0 XDTAGSC#ZDDG$C#ZOaTGSC#ZGSC#ZOTLGSC#LATGSC#LONGSC# 

E XGSC# YGSC#PGSC#QGSC#RGSC#TH£GSC#PHI6SC#PSIGSC# 

F PSIHRC#T£M6SC#TIHGSC#TIMAGSC#TIMBGSC#TIHCGSC#T6SC# 

G FlGSC#6AINGSC#PMAX#PHINAX#ZOTHAX#XDOnAX#YOOTG$C#YDOT6SN 

COMMON/ VAROAT/UX ( 55 ) # L X ( 15 ) 

common/vel/uad#vad#wad# zbad#xootgs#ydqtgs#zdotgs#psihr#zaa#hter# 

1 zaac#zaag#alphad#bad 

COMMON/WIND/Nw#HW(A) #THW(A>#VELW(A)#QTHCW#QVELCW#QZCW#THTW# veltw# 

1 ztw#sigvrw#sigzrw#ssthw#wvelrw#wzrw#xthw#thcw#xvelw#$svelw# 

2 VELCW/ZCW/S1GTHCm#SIGVCW#S1GZCW#S1GW#THC1W#VELCLW# 

3 VELDTW#WTHCW#WTHC1W#WV£LCW#WVELC1W#WZCW# WZC1W#ZC1W# 

A uw#vw#ww 

COMMON/REALTI M/VARCHHG#INTABLS#IVARBUF(5)#ITYPE#FS$( 16) 

LOGICAL FSS#LX 
DIMENSION 0UT(A>#GN( 2# 10) 

DCC ( I)*-?. 501 3*QI-10.«< THEI-THEN )-0.16116*UG 
OCC(2)»0.99872^DC1C2-0*37068*W6+0.31519*W1G 
DC1C2-0CC (2) 

W1G«WG 

: LATERAL CHANNELS 

: ROLL (CYCLIC) 

DCC C3)"7*5^(PGSC-PI)+15.*(PHI6SC+PMIM-PHII)+0.23^VG 
: YAW (DIFFERENTIAL CYCLIC) 

Figure 6. -Program listing of subroutine CONTRL3. 


78/0A/06. 

0B.36.A7 

RTVALT 

A035 

HQN09 

259 

H0N09 

260 

H0N09 

261 

H0N09 

262 

H0N09 

263 

H0N09 

26A 

H0NQ9 

265 

H0N09 

266 

H0N09 

267 

H0N09 

266 

H0N09 

269 

HQN09 

2 70 

H0N09 

271 

H0N09 

272 

H0N09 

273 

H0N09 

27A 

H0N09 

275 

H0N09 

276 

H0NO9 

277 

H0N09 

273 

H0N09 

2 79 

H0N09 

260 

H0N09 

201 

H0N09 

232 

H0N09 

2B3 

H0N09 

26A 

HON09 

285 

H0N09 

286 

H0N09 

287 

H0N09 

266 

H0N09 

289 

H0N09 

290 

H0N09 

291 

H0N09 

292 

H0N09 

293 

H0N09 

29A 

H0N09 

295 

H0N09 

296 

H0N09 

297 

H0N09 

298 

H3N09 

299 

HO NO 9 

300 

H0N09 

301 

H0N09 

302 

H0N09 

303 

HON09 

30A 

H0NO9 

305 

H0N09 

306 

H0N09 

307 

HQN09 

308 

H0N09 

309 

H0N09 

310 

H0N09 

311 

H0N09 

312 

H0N09 

313 

H0NO9 

31A 

H0N09 

315 

H0N09 

316 


Appendix A 


60 

C 


H0N09 

317 



PSIYIG-PSIYG 

H0N09 

318 



PSIYG-PSIGSC-PSIl 

H0N09 

319 



IPIPSIYG.GT.PP) PSIYG-PSIYG-PP2 

H0N09 

320 



IF (PSIYG*LT."PP)PSIYG»PSIYG+PP2 

H0N09 

321 

6t) 


PS1YC«PSIYC+1,4»CT»(PSIY1G^PSIYG) 

H0N09 

322 



DCC(<.)«15**(RGSC~Rl)+17.^PSIY6+PSIYC 

H0N09 

323 



DPL(l)-DCNm 

H0N09 

32A 



DPL(2)«0CN(2) 

H0N09 

325 



DPL(3)»DCNO) 

H0N09 

326 

70 


OPL (^)«DCN(<r) 

H0N09 

327 


C 


H0N09 

328 


c****** HONEYWELL ADAPTIVE CONTROL LAWS 

H0N09 

329 


C 

PITCH AXIS SENSOR INPUTS 

H0N09 

330 


c 


H0N09 

331 

75 


UXI16)«SQRT(XDaTE*XOOTE+YOOTE*YOOTE) 

H0N09 

332 



UX<ll)«ZDOTE 

H0N09 

333 



UX(12)«QI 

H0N09 

33^ 



UX( 13)>THEI 

H0N09 

335 



UX<14)«0C( u 

H0N09 

336 

dO 


UX<15)*0C(2) 

H0N09 

337 



UXU9)»CT 

H0N09 

330 



VBAR*UX<I6)/270* 

H0N09 

339 



UX<10)«VBAR 

H0NO9 

3^0 


c 


H0N09 

3A1 

85 

LATERAL SENSOR INPUTS 

H0N09 

3A2 


C 


H0N09 

3A3 



UX( 37)«PI 

H0N09 

3<»A 



UX( 38 )»R1 

H0N09 

3A5 



0X<39)-PHII 

H0N09 

3A6 

90 


UX(^0)-PSII 

H0N09 

347 


C 


H0N09 

346 


COMMAND INPUTS 

H0N09 

349 


C 


H0N09 

350 



CALL T$IG(C1#C2#UTS) 

H0N09 

351 

95 


VZC-0.4UX(6)*UTS 

H0N09 

352 



THC-THEN4«UX(7)*UTS 

H0N09 

353 



PHlC-0. 

H0N09 

354 



PSIC-0. 

H0N09 

355 


c 


H0N09 

356 

100 

COMMAND DOUBLET 

H0N09 

357 


C 


H0N09‘ 

358 



SDB-0. i 

H0N09 

359 



IFIFSSUOM S0B«UX(17) 

H0N09 

360 



IFIFSSdl) ) S08«-UX(17) 

H0N09 

361 

105 

C 


H0N09 

362 


LOGIC TO APPLY STEP COHMANOS AS PERTURBATIONS 

H0N09 

363 


C 

(CONVERTkD TO RAD* FOR ATTITUDE COMMANDS) 

HO NO 9 

364 


C 


H0N09 

365 



IF(LX(11)> VZC-SOB 

H0N09 

3 66 

110 


1F(LX<12>) THC*TNEN4>SDB«0*017A5 

H0N09 

367 



IF(LX(I3)) PHlCwS0B«O«0l7A5 

H0N09 

368 



IF(LXaA)) P$IC«S0BP0«017A5 

H0N09 

369 



UX(1)»VZC 

H0N09 

370 



UX(2)«THC 

HQN09 

371 

115 


UX (3)»PHIC 

H0N09 

372 



UX(A)-PSIC 

H0N09 

373 


C 


H0N09 

374 


C*A.**A LOGIC TO APPLY GUIDANCE 

H0N09 

375 


C 


H0N09 

376 

120 


1F(,NDT.LX(6)) GO TO "»00 

H0N09 

377 


C 


H0N09 

378 


COMMANDS VZ» THETA#PH1# PSI GO TO UX 1#2#3#<» 

H0N09 

379 


C 


H0N09 

380 



UX(1)«ZD0TGSC 

H0N09 

381 

125 


UX(2)«THE6SC 

H0N09 

382 



UM3)*PH1G$C 

H0N09 

383 



UX (^t)«PSlGSC 

H0N09 

384 



CCNTINU6 

H0N09 

385 


C 


H0N09 

366 

130 


1F(*NQT.FS$(9) ) GO TO 500 

H0N09 

387 



LX (9)-FSS( 12) 

H0N09 

388 



CALL PCMLE 

H0N09 

389 


C 


H0N09 

390 



IF(LX(3)) VBAR-UX(26) 

H0N09 

391 

135 


CALL GAIN2(VBAR) 

H0N09 

392 



IF<LX(A)) CALL GAIn(VBAR) 

H0N09 

393 


C 


H0N09 

394 


PITCH CONTROLLck 

H0N09 

395 


C 


H0N09 

396 


Figure 6. -Continued. 
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Appendix A 


lAO CALL HCON1(DCCU)>OCC(2)) 

occu)-occm«(ux(B)^i* > 

DCC(2)«0CC(2)^(UXf9)>l* > 

5>00 CCNTiNUt 
C^*tt^** LATERAL******* 

1F(FSS(13M CALL HCJ^2(0CC(3)»0CCU)I 
C 

RETURN 

END 


H0N09 

397 

H0N09 

39B 

H0N09 

399 

H0N09 

^00 

NQN09 

<*0l 

HQN09 

A 02 

H0N09 

403 

RTVALT 

4036 

KTVALT 

4037 


SYMBOLIC REFERENCE MAP (R«l) 

ENTRY POINTS 
1 C0NTRL3 


VARIABLES SN TYPE RELOCATION 


13 

AAT 

REAL 


FREQ 

3 

ADT 

REAL 


FREQ 

14 

ALPHAD 

REAL 


VEL 

20 

ALPHAH 

REAL 


HELCOP 

1 

ALT 

REAL 


FREQ 

26 

ANAT 

REAL 


FREQ 

30 

AICT 

REAL 


FREQ 

15 

BAD 

REAL 


VEL 

21 

BETAH 

REAL 


HELCOP 

51 

BLAI 

REAL 

ARRAY 

INERT 

43 

BRI 

PEAL 

ARRAY 

INERT 

35 

BRII 

REAL 

ARRAY 

INERT 

74 

COM 

REAL 

ARRAY 

AR 

124 

C0N2 

REAL 

ARRAY 

AR 

10 

CT 

REAL 


FREQ 

243 

Cl 

REAL 



244 

C2 

REAL 



0 

OAH 

REAL 

ARRAY 

AR 

10 

DAIH 

REAL 

ARRAY 

AR 

20 

DA2H 

REAL 

ARRAY 

AR 

206 

UC 

REAL 

ARRAY 

AR 

0 

OCC 

REAL 

ARRAY 

CONTROL 

44 

DCMAX 

REAL 

ARRAY 

AR 

40 

DCMIN 

REAL 

ARRAY 

AR 

212 

OCN 

REAL 

ARRAY 

AR 

4 

0C1C2 

REAL 


CONTROL 

156 

OHH 

PEAL 

ARRAY 

AR 

162 

OHIH 

REAL 

ARRAY 

AR 

166 

0H2H 

REAL 

ARRAY 

AR 

4 

OLH 

REAL 

ARRAY 

AR 

14 

DLIH 

REAL 

ARRAY 

AR 

24 

DL2H 

REAL 

ARRAY 

AR 

6 

OPL 

REAL 

ARRAY 

CONTROL 

30 

DRLH 

REAL 

ARRAY 

AR 

50 

DRLMAX 

REAL 

ARRAY 

AR 

34 

ORLIH 

REAL 

ARRAY 

AR 

34 

DT 

REAL 


FREQ 

154 

DTRL 

REAL 


AR 

3 

OUDOTH 

REAL 


HELCOP 

0 

OUNOUT 

REAL 

ARRAY 

SPRNT 

5 

OWOQTH 

REAL 


HELCOP 

14 

ENVl 

REAL 


FREQ 

15 

ENV2 

REAL 


FREQ 

16 

ENV3 

REAL 


FREQ 

17 

ENV4 

REAL 


FREQ 

20 

ENV5 

REAL 


FREQ 

21 

ENV6 

REAL 


FREQ 

22 

ENV7 

REAL 


FREQ 

23 

ENV8 

REAL 


FREQ 

24 

ENV9 

REAL 


FREQ 

6 

ERATE 

REAL 


CONST 

155 

FAR 

REAL 


AR 

62 

FI6SC 

REAL 


TRAJECT 

27 

FKT 

REAL 


FREQ 

10 

FSS 

LOGICAL 

ARRAY 

REALTIM 

22 

FXH 

REAL 


HELCOP 

23 

FYH 

REAL 


HELCOP 

24 

FZH 

REAL 


HELCOP 

1 

6 

REAL 


CONST 

63 

GAINGSC 

REAL 


TRAJECT 

216 

GAINH 

REAL 

ARRAY 

AR 

5 

GCT 

REAL 


FREQ 

257 

GN 

REAL 

*UNOEF 


32 

GNOhT 

REAL 


FREQ 

6 

GPT 

REAL 


FREQ 

7 

GRT 

REAL 


FREQ 

33 

GSLT 

REAL 


FREQ 

4 

GST 

REAL 


FREQ 

12 

HRT 

REAL 


FREQ 

11 

HTER 

REAL 


VEL 

1 

Hrf 

REAL 

ARRAY 

MIND 

17 

IDUM 

INTEGER 


SPRNT 

1 

INTA0LS 

INTEGER 


REALTIM 

7 

ITYPE 

INTEGER 


REALTIM 

2 

IVARBUF 

INTEGER 

ARRAY 

REALTIM 

0 

LATE 

INTEGER 


Esr 

41 

LAT6SC 

INTEGER 


TRAJECT 

14 

LATGSN 

INTEGER 


TRAJECT 

1 

LONE 

INTEGER 


fcST 

42 

LONGSC 

INTEGER 


TRAJECT 

15 

lONGSN 

INTEGER 


TRAJECT 

67 

LX 

LOGICAL 

ARRAY 

VAKDAT 

0 

NC 

INTEGER 


CONST 

0 

NW 

INTEGER 


NINO 

253 

OUT 

REAL 

*UN0£F 


11 

PDI 

REAL 


INERT 

11 

POOTH 

PEAL 


HELCOP 

45 

PGSC 

REAL 


TRAJECT 

20 

PGSN 

REAL 


TRAJECT 

14 

PH 

REAL 


HELCOP 

250 

PHIC 

REAL 



51 

PHIGSC 

REAL 


TRAJECT 

25 

PHIGSN 

REAL 


TRAJECT 

1 

PHIH 

REAL 


HELCOP 

1 

PHI! 

REAL 


INERT 

65 

PHI MAX 

REAL 


TRAJECT 

I 

PHIN 

REAL 


NOM 

6 

PI 

REAL 


INERT 

64 

PHAX 

REAL 


TRAJECT 

4 

PP 

REAL 


CONST 

25 

PPH 

REAL 


HELCOP 

7 

PP180 

REAL 


CONST 

10 

PP1801 

REAL 


CONST 

11 

PP2 

REAL 


CONST 

251 

PSIC 

REAL 



52 

PSIGSC 

REAL 


TRAJECT 

26 

PSIGSN 

REAL 


TRAJECT 

2 

PSIH 

REAL 


HELCOP 

7 

PSIHR 

REAL 


VEL 

53 

PSIHRC 

REAL 


TRAJECT 

12 

PSIHRG 

PEAL 


GUIDE 

27 

PSIhRN 

REAL 


TRAJECT 

2 

PSll 

REAL 


INERT 

17 

PSIIHR 

REAL 


HELCOP 

12 

PSIYC 

REAL 


CONTROL 

13 

PSIV6 

REAL 


CONTROL 

241 

PSIY16 

REAL 



11 

PT 

REAL 


FREQ 

31 

PTCT 

REAL 


FREQ 

12 

QDI 

REAL 


INERT 


Figure 6, -Continued. 
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12 

QUOTH 

REAL 


HELCOP 

46 

QGSC 

REAL 


TRAJECT 

21 

OGSN 

REAL 


TRAJECT 

15 

OH 

REAL 


HELCOP 

7 

QI 

REAL 


INERT 

26 

OQH 

REAL 


HELCOP 

15 

OTHCW 

REAL 


WIND 

16 

OVELCW 

REAL 


WIND 

17 

OZCW 

REAL 


WIND 

13 

RDI 

REAL 


INERT 

13 

ROOTh 

REAL 


HELCOP 

2 

RE 

REAL 


CONST 

3 

RE2 

real 


CONST 

47 

RGSC 

REAL 


TRAJECT 

22 

RGSN 

REAL 


tRAJECT 

0 

RGT 

REAL 


FREQ 

16 

RH 

REAL 


HELCOP 

10 

R1 

REAL 


INERT 

2 

RIGT 

REAL 


FREQ 

32 

RLAI 

REAL 

ARRAY 

INERT 

27 

RRH 

REAL 


HELCOP 

27 

RRI 

REAL 

ARRAY 

INERT 

2A 

RRII 

REAL 

ARRAY 

INERT 

23 

RIGSN 

REAL 


TRAJECT 

5 

SCHULER 

REAL 


CONST 

252 

SDB 

REAL 



A6 

SFRI 

REAL 

ARRAY 

INERT 

40 

SFRII 

REAL 

ARRAY 

INERT 

22 

SIGBLAI 

REAL 


INERT 

17 

SIGBRI 

REAL 


INERT 

14 

SI6BR11 

REAL 


INERT 

16 

SI6FRII 

REAL 


INERT 

23 

SIGRLAI 

REAL 


INERT 

20 

SIGRRI 

REAL 


INERT 

15 

SIGRRII 

REAL 


INERT 

21 

SIGSFRI 

REAL 


INERT 

36 

SIGTHCW 

REAL 


WIND 

37 

SIGVCW 

REAL 


WIND 

23 

SIGVRW 

REAL 


WIND 

41 

SIGW 

REAL 


WIND 

40 

SIGZCW 

REAL 


WIND 

24 

SIGZRW 

REAL 


WIND 

25 

SSTHU 

REAL 


WIND 

7 

SSTHWE 

REAL 


EST 

33 

SSVELW 

REAL 


WIND 

6 

SSVELWE 

REAL 


EST 

54 

TEHGSC 

REAL 


TRAJECT 

61 

TGSC 

REAL 


TRAJECT 

30 

TGSN 

REAL 


TRAJECT 

247 

THC 

REAL 



31 

THCW 

REAL 


WIND 

42 

THClW 

REAL 


WIND 

50 

THtGSC 

REAL 


TRAJECT 

24 

THE6SN 

REAL 


TRAJECT 

0 

THEH 

REAL 


HELCOP 

0 

THEI 

REAL 


INERT 

0 

THEN 

REAL 


Non 

20 

THTW 

REAL 


WIND 

5 

Jtik 

REAL 

ARRAY 

WIND 

56 

TIHAGSC 

REAL 


TRAJECT 

57 

TINBGSC 

REAL 


TRAJECT 

60 

TinCGSC 

REAL 


TRAJECT 

35 

TIMEH 

REAL 


FREQ 

55 

TIHGSC 

REAL 


TRAJECT 

0 

UAO 

REAL 


VEL 

6 

UAG 

REAL 


GUIDE 

5 

UAN 

REAL 


NON 

3 

UOOTI 

real 


INERT 

0 

UG 

REAL 


GUIDE 

6 

UH 

REAL 


HELCOP 

2 

UN 

REAL 


NON 

245 

UTS 

REAL 



54 

UW 

REAL 


WIND 

0 

ux 

REAL 

ARRAY 

VAROAT 

1 

VAC 

REAL 


VEL 

7 

VAG 

REAL 


GUIDE 

6 

VAN 

REAL 


NON 

0 

VARCHNG 

REAL 


REALTIH 

176 

VAR12 

REAL 

ARRAY 

AR 

172 

VAR2 

REAL 

ARRAY 

AR 

202 

VAR22 

REAL 

ARRAY 

AR 

242 

VBAR 

REAL 



4 

VOOTH 

REAL 


HELCOP 

4 

VDOTI 

REAL 


INERT 

13 

VEGSN 

REAL 


TRAJECT 

34 

VELCW 

REAL 


WIND 

43 

VEtClW 

REAL 


WIND 

44 

VELDTW 

REAL 


WIND 

21 

VELTW 

REAL 


WIND 

11 

VELW 

REAL 

ARRAY 

WIND 

1 

V6 

REAL 


GUIDE 

7 

VH 

REAL 


HELCOP 

3 

VN 

REAL 


NON 

11 

VNGSN 

REAL 


TRAJECT 

55 

VW 

REAL 


WIND 

246 

VZC 

REAL 



2 

mad 

REAL 


VEL 

10 

WAG 

REAL 


GUIDE 

54 

WAH 

PEAL 

ARRAY 

AR 

7 

WAN 

real 


NON 

5 

WDOTI 

REAL 


INERT 

2 

WG 

REAL 


GUIDE 

10 

WH 

REAL 


HELCOP 

4 

WN 

REAL 


NOH 

64 

WRH 

REAL 

ARRAY 

AR 

45 

WTHCW 

REAL 


WIND 

46 

WTHClW 

REAL 


WIND 

47 

WVELCW 

REAL 


WIND 

50 

WVELCIU 

REAL 


WIND 

26 

WVELRW 

REAL 


WIND 

56 

WW 

REAL 


WIND 

51 

WZCW 

REAL 


WIND 

52 

WZCIW 

REAL 


WIND 

27 

WZRW 

REAL 


WIND 

5 

WIG 

REAL 


CONTROL 

32 

XOOAGSC 

REAL 


TRAJECT 

1 

XODAGSN 

REAL 


TRAJECT 

31 

XD06SC 

REAL 


TRAJECT 

0 

XODGSN 

REAL 


TRAJECT 

67 

XDOHAX 

REAL 


TRAJECT 

3 

XOOTE 

REAL 


EST 

3 

XD0T6 

REAL 


GUIDE 

4 

XOPTGS 

REAL 


VEL 

33 

XDOTGSC 

REAL 


TRAJECT 

2 

XDOTGSN 

REAL 


TRAJECT 

10 

XDTAOE 

REAL 


EST 

34 

XOTAGSC 

REAL 


TRAJECT 

3 

XDTAGSN 

REAL 


TRAJECT 

30 

XOTAVH 

REAL 


HELCOP 

10 

XOTNGSN 

REAL 


TRAJECT 

13 

XE 

REAL 


EST 

43 

XGSC 

REAL 


TRAJECT 

16 

XGSN 

REAL 


TRAJECT 

25 

XLNT 

REAL 


FREQ 

30 

XTHN 

REAL 


WIND 

32 

XVELW 

REAL 


WIND 

4 

YOOTE 

REAL 


EST 

4 

YOOTG 

REAL 


GUIDE 

5 

YDOTGS 

REAL 


VEL 

70 

YDOTGSC 

REAL 


TRAJECT 

71 

YD0T6SN 

REAL 


TRAJECT 

11 

YDTAOE 

REAL 


EST 

12 

YOTEGSN 

REAL 


TRAJECT 

14 

YE 

REAL ' 


EST 

44 

YGSC 

REAL 


TRAJECT 

17 

YGSN 

REAL 


TRAJECT 

10 

ZAA 

REAL 


VEL 

12 

ZAAC 

REAL 


VEL 

13 

ZAAG 

REAL 


VEl 

60 

ZAH 

REAL 

ARRAY 

AR 

3 

ZBAU 

REAL 


VEL 

11 

Z6G 

REAL 


GUIDE 

35 

ZCW 

REAL 


WIND 

53 

ZCIW 

REAL 


WIND 

35 

ZDOGSC 

REAL 


TRAJECT 

4 

ZOOGSN 

REAL 


TRAJECT 

5 

ZDOTE 

REAL 


EST 

5 

ZDOTG 

REAL 


GUIDE 

6 

Z00T6S 

REAL 


VEL 

36 

ZDOTGSC 

REAL 


TRAJECT 


Figure 6. -Continued. 
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5 ZOQTGSN 

REAL 


<rO 2DTLGSC 

PEAL 


6 ZOTIGSN 

REAL 


37 ZGSC 

REAL 


70 ZRM 

REAL 

ARRAY 

EXTERNALS 

TYPE 

ARCS 

GAIN 


1 

HCONl 


2 

PCMLE 


0 

TSIG 


3 

STATEMENT LABELS 



155 AOO 



COMMON BLOCKS 

LENGTH 


EST 

13 


AR 

1A6 


CONST 

10 


CONTROL 

12 


FREO 

30 


GUIDE 

11 


HELCOP 

25 


INkPT 



NOM 

8 


SPRNT 

16 


TRAJECT 

58 


VAROAT 

70 


VEL 

1<» 


RIND 

A7 


REALTIM 

2A 



TRAJECT 

12 

ZDTAOE 

REAL 

TRAJECT 

66 

ZOTMAX 

REAL 

TRAJECT 

2 

2E 

REAL 

TRAJECT 

7 

ZGSN 

REAL 

AR 

22 

ZTH 

REAL 


6AIN2 

HC0N2 

SORT REAL 


200 500 


EST 

TRAJECT 

EST 

TRAJECT 

WIND 


1 

2 

1 LIBRARY 


STATISTICS 

PR06RAH LENGTH BOBS X95 

CM LABELED COMMON LENGTH 1020S 52B 

730006 CM USED 


Figure 6. -Concluded 


Appendix A 


SUfekOUTlNF ‘lC2M 73/7<i QPT-1 


FTN 4,64-<ib2 76/0^/06. att«36«<»7 


1 


5 


10 


lb 


20 


SUBkDUTINc HCUNUUCl#UC2) 

C 

CUhhQN/CI/ XCl»XC2 
COMHQN/VAkDAT/ UX(55 tX(l5) 

LOGICAL LX 
V2-UX( U ) 

Q-UXU2) 

TH»UX(n) 

VZC-UXU ) 

THC-UX(2 > 

DT-0XU9) 

C PITCH CONTROL 

- UCl-XCl ♦ 0X(5U^VZC + UX{52)*0,8*THC 
UC2-XC2 ♦ UX(53)^VZC ♦ UX ( 5<t )*0. 8*THC 
UCl-JCl ♦ UXC^D^VZ ♦ UX(<»2)*Q ♦ UX(^3)*TH 
UC2-UC2 ♦ UX(46)*VZ ♦ UX(<r7)*Q ♦ UX(^8)^TH 
DXC1-UX(^^)*( VZC-VZ) ♦ UX(<»5)*0*8*(T4C-TH) 
0X:2*LX(^9)*<VZC-VZ) ♦ UX(50 )pO*8^(THC-TH) 
C 

XCl* XCl ♦ OXCl+DT 
XC2« XC2 ♦ 0XC2+0T 
RETURN 
END 


H0N09 

7^0 

H0N09 

7<»1 

HQN09 

7A2 

H0N09 

7<»3 

H0N09 

74^ 

H0N09 

7A5 

H0N09 

7A6 

H0N09 

7A7 

H0N09 

7A8 

H0N09 

7A9 

HO NO 9 

750 

H0N09 

751 

H0N09 

752 

H0N09 

753 

H0N09 

7 5A 

H0N09 

755 

H0N09 

756 

H0N09 

757 

H0N09 

758 

rlON09 

759 

H0N09 

760 

H0N09 

761 

H0N09 

762 


SYNaOLIC ktHERENCt NAP (R-U 


ENTRY POINTS 
3 HCONl 


I 


variables sk type 


62 

OT 

REAL 

6A 

DXC2 

REAL 

56 

0 

REAL 

61 

THC, 

REAL 

0 

0C2 

REAL 

55 

VZ 

REAL 

0 

XCl 

REAL 


RELOCATION 


F.P. 

Cl 


SUeRULTINE HCONl 73/7<» OPT-1 


63 

OXCl 

REAL 



67 

LX 

LOGICAL 

ARRAY 

VAROAT 

57 

Tri 

REAL 



0 

UCl 

REAL 


F*P. 

0 

UX 

REAL 

ARRAY 

VAROAT 

60 

VZC 

REAL 



I 

XC2 

REAL 


Cl 


FTN ^m6*MbZ 78/0-^/06. 08.36*<»7 


COMhON BLOCKS LENGTH 
Cl 2 

VAROAT 70 


/ 


statistics 

PPOGRAN LENGTH 
Ch LABELED C3HH0N LENGTH 
730008 CH USiC 


65S 53 

HOB 72 


Figure 7, -Program listing of subroutine HCONl. 
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SueRUlTlN? HC0N2 73/7^ QPT-i FTN 78/0<i/06* 08.36*<»7 


1 


SUBROUTINE HC UN2 < UC3> UC4 ) 

H0N09 

764 


C 


H0N09 

765 



COMMON /C0NST/NC#6#RE#Rt 2# PP#$CHULER#ER ATE# PP IS 0#PP 1801# PP2 

H0N09 

766 


C 


H0N09 

767 

5 


C0MM0N/CI/XC1#XC2#XC3#XC4 

^ON09 

768 



C0NM0N/VARDAT/UX(55)#LX(I5) 

HQN09 

769 



LOGICAL LX 

H0N09 

770 


C 


HON09 

771 



PHIC«UX(3) 

H0N09 

772 

10 


PSIC-UXC4) 

HONOV 

773 



P-UX(37) 

HQN09 

774 



R-UX(38) 

H0N09 

775 



PHI-UX( 391 

HON09 

776 



PSI»UX(40) 

HON09 

777 

lb 


0T-UX(19) 

H0N09 

778 


C 


H0N09 

779 



DUM-PSIC-PSI 

^0N09 

780 



IF (ABS <OUM)«LT.PP) GO TO 100 

H0N09 

781 



1F(DUM.GT.PP.AN0#PSI.LT.0. ) XC3-XC 3- 7* 4PP2 

HOMO 9 

762 

20 


IF<DUM.GT.PP.AN0#PS1.LT*0*I XC4» XC4-50.*PP2 

rl0N09 

783 



IF(DUM#GT«PP«AN0«P$IC#GT»0*)XC3«XC 3-4» 540*84PP2 

N0N09 

784 



IF(0UM.GT.PP*AN0«PSIC*GT*0«) XC4-XC4-34* 3d«0#d «PP2 

HQN09 

785 



1F(OUM#LT.-PP#ANO#PS1*GT#0#) XC3«XC347*4PP2 

H0N09 

786 



IF(DUH.LT#-PP«AN0#PSI*GT*O«) XC4«XC4450«4PP2 

' H0N09 

787 

25 


IF ( OUM.LT*-PP. AND* PS IC •LT.O , } XC 3»XC 3^4«540# 8*PP2 

H3N09 

780 



IF (DUM.LT#-PP.ANO. PS IC*LT#0. ) XC4«XC4+34*38^0* 84PP2 

H0N09 

789 



IF(OUM.GT.PP) DUM-0UM-PP2 

H0N09 

790 



IF(OUM.LT.-PPI 0Urt»0UM^PP2 

H0N09 

791 



100 CONTINUE 

H0N09 

792 

30 

c 


HQN09 

793 



UC3-XC3+15.94*0*8*PHIC+4.5*0,8*PSIC 

H0N09 

794 



UC4«XC4^3*3140*84PHIC*34.38*0*8*PSIC 

H0N09 

795 



UC3»UC3-9*«P-3,*R-23.5*PHI-7**PSI 

H0N09 

796 



UC4»UC442.3754P-24.4R44*94PHI-50*#PSI 

H0N09 

797 

35 

c 


H0N09 

790 



Si*CPHIC-PHI)^.8 

H0N09 

799 



S2*DUM««8 

H0N09 

800 



XC3«KC3+(15#94#S1^4*5*S2)40T 

H0N09 

301 



XC4«XC44(3.31^S1434# 39*S2>*0T 

H0N09 

802 

40 

c 


H0N09 

803 



RETURN 

H0N09 

004 



END 

H0N09 

805 


SYheOlIC R&FERtNCt MAP (R*U 

tHTPY POINTS 
3 HC0N2 


VARIABLES SN TYPE RELOCATION 


203 

DT 

real 



204 

DUN 

REAL 



6 

ERATE 

REAL 


CONST 

1 

G 

REAL 


CONST 

67 

LX 

LOGICAL 

ARRAY 

VAROAT 

0 

NC 

INTEGER 


CONST 

177 

P 

REAL 



201 

PHI 

REAL 



175 

PHIC 

REAL 



4 

PP 

REAL 


CONST 

7 

PPIBO 

REAL 


CONST 

10 

PP160I 

REAL 


CONST 

11 

PP2 

REAL 


CONST 

202 

PSI 

REAL 



176 

PSIC 

REAL 



200 

R 

REAL 



2 

RE 

REAL 


CONST 

3 

RE2 

REAL 


CONST 

5 

SCHULER 

REAL 


CONST 

205 

SI 

REAL 



206 

S2 

REAL 



0 

UC3 

REAL 


F#P# 

0 

UC4 

REAL 


F*P. 

0 

UX 

REAL 

ARRAY 

VAROAT 

0 

XCl 

REAL 


Cl 

1 

XC2 

REAL 


Cl 

2 

XC3 

PEAL 


Cl 

3 

XC4 

REAL 


Cl 

:nline 

FUNCTIONS 

TYPE 

ARGS 








ABS 

REAL 

1 INTRIN 







STATEMENT LABELS 
111 100 

COMMON BLOCKS LENGTH 
CONST 10 

Cl <» 

^ VARDAT 70 

STATISTICS 

PROGRAM LENGTH 2078 13& 

CM LABELED COMMON LENGTH 12AB 64 

" 73000B CM USED 


Figure 8* -Program listing of subroutine HCON2, 
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PROGRAM LISTING OF EXPLICIT-GAIN COMPUTATION 

This appendix consists of the program listings of subroutine GAIN (Fig. 9) 
and GAIN2 (Fig. 10). 


SoPROUTiNt GAIN 

73/7A OPT-1 

FTN A.6^A52 

7B/0<i/06* 

0 d. 36 . 4 7 

1 ' 

SUBROUTINE GAINU) 


H0NO9 

666 


COMMJN/GTAB/GNUA, 13) 


H0N09 

669 


COMMON/ VAR0AT/UX( bb)#LX( 15) 


H0N09 

670 


LOGICAL LX 


H0N09 , 

671 

C 



H0N09 

672 


' TABLE LOOK-UP FOR PITCH AXIS GAIN 


H0NO9 

673 

C 

UX(27)«ACCURACY 


H0N09 

67^ 

c 



H0NO9 

67b 


IS»3 


HQN09 

676 

10 

IF(UX<27).LT.0.6) IS»1 


H0N09 

677 


IF(UX(27).LT.O«J) IS-2 


H0N09 

678 


IF(.N0T.LX<3)) IS-2 


H0N09 

679 


V8«A,*2+2,b 


H0N09 

680 


IVB-VB 


H0N09 

6B1 

lb 

ifuvb.lt*!) IVB«1 


H0N09 

682 


IF(XVB.GT*6) IVB-b 


H0NQ9 

6B3 


I»IVB*6*a$-l)l-l 


H0N09 

68<» 


1FUS.E0.3) I«1 


H0N09 

68 5 


IFUX(IO)) I-l 


H0N09 

686 

20 

UX<A1)«GNU#1) 


H0N09 

687 


UX(A2)«6N<2U) 


H0N09 

688 


UX(A3)«GN(3#I) 


HQN09 

689 


UX(^A)«6N<A»I) 


H0N09 

690 


UX(Ab)-GN(5»l) 


H0N09 

691 

2 b 

UX(bl)-GNUU) 


H0N09 

692 


UX(b2)*GN<7U) 


H0N09 

693 


UX (A6)>GN(6>I) 


H0N09 

69<» 


UX ( A7)-Gi>M9>I } 


HJN09 

69 b 


UX< A6)"GN( 10# I ) 


H0N09 

696 

30 

UX(A9)»GN(I1#I) 


H0N09 

697 


UX( bO)»GN( i2U ) 


H0N09 

698 


UX (b3)-GNU3#I ) 


H0N09 

699 


UX(b^)«GN(lA#I ) 


H0N0<^ 

700 


RtTURN 


H0N09 

701 

3b 

END 


H0N09 

702 


StHBOLlC RfcFERENCE HAP (R-1) 


ENTRY POINTS 
3 GAIN 


VARIABLES 

SN TYPE 

ReiGCATION 






0 

ON 

REAL 

ARRAY GTAB 

7A 

I 

INTEGER 



71 

IS 

INTEGER 


73 

IVS 

INTEGER 



67 

LX 

LOGICAL 

ARRAY VAKOAT 

0 

ux 

REAL 

ARRAY 

VAROAT 

72 

VB 

REAL 


0 

2 

REAL 


F. P. 


CdNMON BLOCKS LENGTH 
GTAb 182 

VaROAT 70 


statistics 

PROGRAM LENGTH 75B 61 

CH LA0FLEO COMMON LENGTH 37<»B 2&2 

730000 CM bStC 


Figure 9. -Program listing of subroutine GAINI. 
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SJB«0LTINE GAIN2 


73/7<t OPT-i 


PTN 78/0<i/06. 08*36. <i7 


1 


5 


10 


15 


20 


SUBROUTINE GA1N2(Z1) 

C0HM0N7VARDAT/ UX ( 55 ) # LX ( 15) 

LOGICAL LX 

UXC^U*-0.25AZ1 

IF(Zl.GE.0.5) UX(<»1)— 0.125 

UX(^2)«- 7.6^1. 5«Zl 

UX(<*3)«-27.+6.*Zl 

UX(<»^)«0.1<»»Z1 

IF(Z1.GE.0.5) UX(4<»)»0.07 

UX(<i5)«l7.5-5.*Zl 

UX(46)«0.3-0.12«Z1 

UX(47) — l.*Zl 

IF<Z1.G6.0.5)UX<<»7) — .5 

UX(<»8)«-18.*Z1 

UX(<»9)— 0.23>0.08AZ1 

UX(50)»2.4>Z1 

UX(51)«UX(A4) 

UX(52)-UX«^5) 

UX(53)«UXf^9) 

UX(5^)«UX(50) 

RETURN 

END 


H0N09 70<» 
H0N09 705 
N0N09 706 
H0N09 707 
H0N09 70d 
H0N09 709 
N0N09 710 
H0N09 711 
rt0N09 712 
H0N09 713 
H0N09 71^ 
H0N09 715 
H0N09 7X6 
H0N09 717 
H0N09 718 
H0N09 719 
H0N09 720 
H0N09 721 
H0N09 722 
H0N09 723 
H0N09 724 
H0N09 725 


SYNBOLIC REFERENCE HAP (R«l) 


ENTRY POINTS 
3 GAIN2 


VARIABLES 
67 LX 
0 Z1 

SN TYPE 
LOGICAL 
REAL 

RELOCATION 
ARRAY VARDAT 

F.P. 

0 UX 

REAL 

ARRAY 

VAROAT 

COHHON BLOCKS 
VAROAT 

LENGTH 

70 






STATISTICS 

PR06RAN LENGTH 
CH LABELED COHHON LENGTH 

738 59 

1068 70 






Figure 10. -Program listing of subroutine GAIN2. 
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APPENDIX C 

FLOW CHART FOR REAL-TIME PORTION OF 
PCMLE PROGRAM 

This appendix consists of the flowchart for subroutine PCMLE (Fig. 11). 



Figure 11. - Flowchart for subroutine PCMLE, 
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Figure 11. -Continued, 















Appendix C 


Subcycle 2; Significance tests and change logic 



Figure 11. -Continued 
























APPENDIX D 

LISTING OF PCMLE REAL-TIME SOFTWARE 


This appendix is comprised of the program listings of PCMLE, FILT, SENS, 
ACUM, FH, and TSIG, as presented in Figures 12 through 17, respectively. 


SuBROUTiNt PCMLE 7i/7<» OPT*i >TN 

7U/04/06. 

o 

CD 

• 

O' 

• 

i 


SUBROUTINE PCMLE 

HQN09 

1213 



COMHON/VAROAT/UX(55)#LX(15) 

H0N09 

1214 



COMMON/OAT/X(5#19)# F < 5# 22) » D< 5# 22 )#E IP# E 2Pf S IGSQ 

H0NO9 

1215 



C0HH0N/SEN$P/XS(4#19)#0S( 5# 4# 22 )# GE (4 )# GL ( 4 )# 6S0£ ( 10) # GSQL (10) 

H0NO9 

1216 



1 #GK(5»4#12) 

H0N09 

1217 



COMhON/lPIC/ NC#NP#aS#CHNGCH#CHSIG#AN$ 

H0N09 

1210 



1 #TPRINT#M0D£#S1GUT0#ZP1#ZP2#ZP3#ZP4#T1ME#RTJS 

H0N09 

1219 



2 #RTJG#THRTJ6#RTJC#THRTJ€#RTJZ#THRTJZ»DTPR1N 

HQN09 

1220 



3 #ZP 1HAX#ZP2MAX#ZP3MAX#ZP4MAX 

HON09 

1221 

10 


4 #ZP1HIN#ZP2MXN>ZP3NIH>ZP4MIN 

H0N09 

1222 



5 #ZP( 5#4)#GSQL0(4)#Y(3)#U(2) 

H0N09 

1223 



6 #XVZ(2)#XQ(2)#XTH(2)#XDa(2)#XDC(2)#PRTIHE 

H0N09 

1224 



C0MM0N/MEAS/YP(3)»UP(2) 

HON09 

1225 



CQHM0N/SUBRH/DTHP#0HP#WHPfClHP#C2HP \ 

HQN09 

1226 

15 


C0HMQN/UTEST/DTU#S£E0#WUT#0UT#SIGUT#UT10#UT20#UTMAX#W2UT#TZWUT# 

H0NO9 

1227 



1 UT1#UT2#GAMUT0 

H0N09 

1220 



C0NM0N/D£V1C£/NR£A0# MFUE 

H0N09 

1229 



C0MH0N/REALT1N/VARCHNG#ZNTABLS#IVARBUF(5)#1TYPE#F$SU6) 

HON09 

1230 



DIMENSION TL(5)#TJ(5)#DZP(4)*SS(5) 

H0N09 

1231 

20 


LOGICAL FSS 

H0N09 

1232 



LOGICAL LX# N0CHC#N0EST#PRINT#MLE 

H0N09 

1233 



INTEGER CHNGCH»CHSIG 

H0N09 

1234 


C 


H0N09 

1235 



EQUIVALENCE (TL( 1)#X(1# 15) )# (TJ( 1)#X(1#14) ) 

H0N09 

1236 



EQUIVALENCE ( R 2# GSQL ( 2 > ) # ( R 3# GSQL ( 3 ) ) # (R4#G$QL(4) > 

H0N09 

1237 



1 #(R6#GS0L(6) )#(R7#GS0L(7))#(R9#GS0L(9)) 

H0N09 

1238 


C 


HQN09 

1239 



TMINCH«2» 

H0N09 

1240 


c 


HQN09 

1241 

30 

c 


H0N09 

1242 


c 

REAL TIME MOOS 

H0N09 

1243 


c 


H0N09 

1244 



N0CHC«LX(1) 

H0N09 

1245 



N0EST«LX(2) 

H0N09 

1246 

35 


PRINT-LX(5) 

H0N09 

1247 



MLE»LX(9) 

H0NO9 

1240 


G 


H0N09 

1249 



SIGUT-SIGUTO 

H0N09 

1250 


c 


H0N09 

1251 

40 


Y(1)»UX(11) 

H0M09 

1252 



Y(2)»UX(X2) 

HOMO 9 

1253 



Y(3)>UX(13) 

H0N09 

1254 



U(1)-UX(14) 

H0N09 

1255 



U(2)»UX(15) 

HON09 

1256 

45 

c 


H0N09 

1257 

c 


H0N09 

1250 



UX( 20)«-*5^.25tJi 

H0N09 

1259 



UX(21)-X(J$#2)/20# 

H0N09 

1260 



UX(22)«X< JS#3)45.73 

HQN09 

1261 

50 


UX(23)«TJ( D/400* 

HQN09 

1262 



UX(24)«TJ(2)/400* 

H0N09 

1263 



UX(25)«T4(3)/400« 

H0N09 

1264 



UX(26)«ZP1 

H0N09 

1265 



UX(2B)*ZP2 

HON09 

1266 

55 


UX(29)«ZP3 

H0N09 

1267 



UX(30)«ZP4 

H0N09 

1260 


Figure 12. -Program listing of subroutine PCMLE. 
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UXl31)-XU$fl) 

H0N09 

1269 


UX(32)>X(JS#<r> 

H0N09 

1270 


C 

H0N09 

1271 

60 

C 

H0N09 

1272 


C OUTPUT VARIABLES IN RUNNING NOOE (6ATCH ONLY) 

H0N09 

1273 


2F(«NnT*PRlNT) GO TO 1100 

H0N09 

127A 


PRTIME-PRTIHEfOTHP 

H0N09 

1275 


IF (PRUNE. LT.TPRINT) GO TO 1100 

HQN09 

1276 

65 

WRITE(HF1LE^630) PRT INE #4$ #S IGSQ 

H0N09 

1277 


taRlTE(NFlLE>631) (TL(I)U-1#NC) 

H0N09 

1278 


WRITE(HFILE#632) ( TJ ( I ) . !■ 1. NC ) 

H0N09 

1279 


WRITE(HF1LE^642) ( SS ( I ) * I- 1# NC ) 

HQN09 

1280 


hRIT£(NFlLE#6^3) ZININ 

NQN09 

1281 

70 

WRITE (NFILE#636) YP( 1 )« YP( 2 ) # YP( 3)#UP( 1 ) «UP ( 2 ) 

H0N09 

1282 


WRITE(NFILE/633) 

H0NO9 

1263 


WRITE(NFILE#634) (X(JSU)> I"1U7) 

HQN09 

1284 


WRITE(NFILE^637) 

HQN09 

1285 


WRITE (NFILE#6 3<>) ( XS ( I )U«1# 11 ) 

HQN09 

1286 

75 

WRITE (NPILE #63^) ( XS ( 2U )# I- 1# U ) 

H0N09 

1287 


s WRITF(NFILE#634) ( XS ( 3/ I ) » I« 1. 11 ) 

H0N09 

1288 


WRITE (NFILE#63^> (XS (^# I ) U«1* 11 ) 

H0N09 

1289 


WRITE(NFILE#636) 

H0N09 

1290 


WR1TE(NF1LE«63A) ( GL (I ) # I* l^NP) 

H0N09 

1291 

ao 

WRITE (NFILE^6 34) (GSQL ( I ) lUO ) 

H0N09 

1292 


IF(NOEST) GO TO 1090 

H0N09 

1293 


IF(T1NE«LT.TNINCH) GO TO 1090 

H0N09 

1294 


Cl • SQRT(C1/0ET) 

H0N09 

1295 


C5 • SQRT(C5/DET) 

H0N09 

1296 

65 

Ce - $0RT(C6/0ET) 

H0N09 

1297 


CIO- SQRTICIO/OET) 

H0N09 

1298 


C2 ■ C2 /(C1*C5 ♦OET) 

H0N09 

1299 


C3 * C3 /(C1*C8 ♦OET) 

HON09 

1300 


• C^ /(Cl^ClO^OET) 

H0N09 

1301 

90 

C6 • C6 /(C5*C8 ♦OtT) 

H0N09 

1302 


C7 • C7 /(C5*C10»0ET) 

H0H09 

1303 


C9 - C9 /(C8«C10*D£T) 

H0N09 

1304 


1092 WRITE(MFILE#6^0) 

HQN09 

1305 


WRITE (NFILE#634) Cl# C2#C3# C^#C5# C6,C7#C8# C9#C10 

H0N09 

1306 

95 

1090 CONTINUE 

HQN09 

1307 


WRlTt(NFILE#639) 2 PI# ZP2# ZP3# ZP<» 

H0N09 

1308 


C 

H0N09 

1309 


606 F0RHAT(6E20.b) 

H0N09 

1310 


630 F0RNAT(////1X#11HTIH£ » #E12.5« 6H JS - #12# 

H0NO9 

1311 

100 

1 <»X#6HSI6SO - #E12.5^2X#2H«4>, /) 

'IQNOG 

1312 


631 F0KNAT(1X#11HTL - #9£12o5) 

H0N09 

1313 


632 FORNAT( IX# IIHTJ • »9E12o5) 

HQN09 

1314 


633 F0RHAT(1HQ#22H$TATE ESTIMATES X(JS> > 

H0N09 

1315 


63<» F0RNAT(10E12.5> 

H0N09 

1316 

105 

636 FORNATC 1H0#16HY1#Y2# Y3#U1#U2 «#6E12.5> 

H0N09 

1317 


637 F0RHATUX#13HSENSiriVlTIES) 

HON09 

1318 


638 F0RNAT(1X#5HGRA0S ) 

HQN09 

1319 


639 FORMAT(1HO#11HZP • #9E12.5) 

H0N09 

1320 


6A0 FORNATC IX# 8HACCURACY ) 

HQN09 

1321 

no 

642 FORMAT ( IX# IIHPR • #9F12.9» 

H0N09 

1322 


643 FQRNAT(1X#11HZ1NIN • #El2.5) 

H0N09 

1323 


1099 CONTINUE 

HQN09 

1324 


TPRINT-TPRINTf DTPRIN 

H0N09 

1325 


1100 CONTINUE 

H0N09 

1326 

U5 

C 

H0N09 

1327 


C TEST SIGNAL 

H0NO9 

1328 


CALL TSIG(ETA1#ETA2#UT) 

H0N09 

1329 


C 

H0N09 

1330 


C 

H0N09 

1331 

120 

C HIGH PASS INPUTS 

HO NO 9 

1332 


CALL FH(Y(l)#YP(l)#XVZ) 

H0N09 

1333 


CALL FH(Y(2)#YP(2)#XQ) 

H0N09 

133A 


CALL FH(Y(3)#YP(3)#XTH) 

H0N09 

1335 


CALL FH(U( 1)#UP< 1)#X0B ) 

H0N09 

1336 

125 

CALL FH(U(2)#UP(2)#XDC ) 

H0N09 

1337 


C 

H0N09 

1338 


C PARALLEL CHANNELS 

H0N09 

1339 


IF(.NOT.HLE) RETURN 

H0N09 

13^0 


TINE-TIME + DTHP 

HQN09 

13A1 

130 

ANS-E1P«ANS * 1. 

H0N09 

1342 


60 TO (1210#1209#1208#1207#1206) # NC 

H0N09 

13A3 


1206 CALL FILT(5) 

H0N09 

1344 


1207 CALL F1LT(4) 

HON09 

1345 


1206 CALL F1LT( 3) 

H0N09 

1346 

135 

1209 CALL FILT(2) 

H0N09 

1347 


Figure 12. -Continued. 
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1210 CALL flLT(l) 

H0N09 

1348 


C 

H0N09 

1349 


C SENSITIVITIES 

H0N09 

1350 


GO TO (1304#1303#1302»1301)» NP 

H0N09 

1351 

1^0 

1301 CALL SENS(«IS»4) 

H0N09 

1352 


1302 CALL SENS<JS#3) 

HON09 

1353 


1303 CALL SENSUS#2) 

H0N09 

1354 


1304 CALL SENS(JS#U 

MONO 9 

1355 


C 

H0N09 

1356 

145 

C LIKELIHOOD ACCUNOLATION 

H0N09 

1357 


CALL ACUN(JS) 

H0N09 

1358 


C 

H0N09 

1359 


900 CONTINUE 

H0N09 

1360 


C 

H0N09 

1361 

150 

C 

H0NO9 

1362 


. C BRANCH TO LOW RATE OPERATIONS 

1363 


C 

H0N09 

1364 


GO TO (910#920#930#940#9501# NODE 

H0N09 

1365 


C 

H0N09 

1366 

155 

C CYCLE 1. HIN-L CHANNEL SELECTION 

H0N09 

1367 


C 

H0N09 

1368 


910 CONTINUE 

H0N09 

1369 


TLMAX • -I.EIO 

H0N09 

1370 


TLHIN • I.EIO 

H0N09 

1371 

160 

TJhlN-l.ElO 

H0N09 

1372 


00 911 I«1»NC 

HON09 

1373 


Sl-TJU) ♦ SIGSa*ANS^F(I,13) 

H0N09 

1374 


TL(n- SI 

H0N09 

1375 


IF(S1.GT*TLMAX) TLNAX-Sl 

H0N09 

1376 

165 

IF(Sl.LT.TLMIN) TLNIN-Sl 

H0N09 

1377 


IF(TJ<I)«GT«TJN1N1 60 TO 911 

HON09 

1376 


TJMIN • TJ(1) 

H0N09 

1379 


JSTENP • 1 

H0N09 

1380 


911 CONTINUE 

H0N09 

1361 

170 

NODE • 2 

H0N09 

1362 


RETURN 

H0N09 

1383 


C 

H0N09 

1384 


C CYCLE 2* SIGNIFICANCE TESTS AND CHANGE LOGIC 

H0N09 

1385 


C 

H0N09 

1386 

175 

920 CONTINUE 

H0N09 

1387 


CHN6CH«0 

H0N09 

1388 


TJE«0«01<»6SgL(l) 

H0N09 

1389 


C 

H0N09 

1390 


IF(TJE*GT*RTJS4<TJ(JS)) GO TO 925 

H0N09 

1391 

180 

SIGSQ* TJ(jS)/(ANS ♦ ANS) 

H0N09 

1392 


IF(SIGSQ«6T«2.) SIGSQ-2* 

H0NO9 

1393 


IF(SIGS0.LTtlt0E-04) S IGSQ-l* OE-04 

H0N09 

1394 


C 

H0NO9 

1395 


925 IF(NOCHC) 60 TO 926 

H0N09 

1396 

185 

IF<TL(JS) - TL(JSTEHP),6T,(THRTJC>RTJC^TJE)) CHNGCH-l 

H0N09 

1397 


C 

H0N09 

1396 


926 continue 

H0N09 

1399 


DO 927 I*1#NC 

H0N09 

1400 


IF(TLU) - TL( JS).LT«(THRTJZ ♦ RTJZ«TJEM ZlNlN«ZPa#l) 

HON09 

1401 

190 

ssm- <TL(j$i- TLan^o*5 

H0N09 

1402 


1FCSS(I)«GT.0») SS(l)«i0« 

H0N09 

1403 


IF(SS(I).LT.-234) SS(I)*-23# 

H0N09 

1404 


1094 SS ( I)-EXP(SS(I) ) 

H0N09 

1405 


927 CONTINUE 

H0N09 

1406 

195 

UX(33)*$I6SQ 

H0N09 

1407 


MODE- 3 

H0N09 

1406 


RETURN 

H0N09 

1409 


C 

H0N09 

1410 


C CHANNEL TRANSFERS 

H0N09 

1411 

200 

C 

H0N09 

1412 


930 CONTINUE 

H0N09 

1413 


IF(CHNGCH.EQ.O) GO TO 939 

H0N09 

1414 


IF(TIHe.GT.TllINCH) GO 10 9301 

H0N09 

1415 


00 9303 I*l#10 

H0N09 

1416 

205 

GSQL<I)«0* 

H0N09 

14X7 


9303 G$QE(I)«0* 

H0N09 

1418 


DO 9302 1«1^4 

H0N09 

1419 


CZ^(U«0* 

H0N09 

1420 


GL(1)«0* 

H0N09 

142* 

210 

9302 GE(I)«0« 

H0N09 

1422 


60 TO 934 

H0NO9 

1423 


9301 CONTINUE 

HQN09 

1424 


C 

H0N09 

1425 


C CHANNEL TRANSFER SEQUENCE INRI 

H0N09 

1426 


Figure 12. -Continued. 
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Zl'y 

C 


SCALE GRAO SO L 

H0N09 

1427 




Sl«l* ^ ZPUS#U/ZP( JSTEHP#U 

H0N09 

1420 




GSQUU- 6SQL(1)«SI«S1 

H0N09 

1429 




GSQU2I- GSQU2)«S1 

H0N09 

1430 




6SQL(3)« GSQL(3)4S1 

H0N09 

1431 

220 



GSQU4)« GSQL(4)«S1 

H0N09 

1432 


C 



H0N09 

1433 




Sl-GSQLU) ♦ GSQLO(I) 

HON09 

1434 




S2«GSQL(5> G$QL0(2) 

HON09 

1435 




S3-GSQLC6) ♦ 6SQLOO) 

H0N09 

1436 

22*) 



S4«GS0L(10I^ GSQL0C4) 

H0N09 

1437 


C 


NEW GRAO L 

H0N09 

1438 




02PU)»-ZP1 + ZP(J$TEMP#l) 

ri0N09 

1439 




0ZP<2)— ZP2 ♦ ZP(JSTENP#2) 

ri0N09 

1440 




0ZPC3) — ZP3 ♦ ZP(J$T£HP»3) 

H0N09 

1441 

230 



DZP<4) — ZP4 ♦ ZP(JSTEHP#4) 

H0N09 

1442 




GLCl)«Si ♦0ZPa)4>6SQLC2>P0ZP(2)^GSQL(3>«0ZP(3)46S0Lf 4)90ZP(4) H0N09 

1443 




GL(2)«G$0L(2)P0ZP<1I^S2 P0ZP(2)^GSQL(6)P0ZP<3)4GSQLC7)«0ZP(4) HQN09 

1444 




GL(3)«6SQL(3)P0ZPa)4^G$QU6)«0ZP(2)«$3 POZPC 3)«GSQL<9)PDZP(4) H0N09 

1445 




GU4)-GSQL(4)P0ZP(1)^GSQL<7I«0ZPI2KGSQL(9)P0ZP(3)«$4 

P0ZP(4) H0N09 

1446 

235 


934 

JS»JSTENP 

H0N09 

1447 


C 



H0N09 

1448 



939 

CONTINUE 

H0N09 

1449 




NODE- 4 

H0N09 

1450 




RETURN 

H0N09 

1451 

2^0 

c 



H0N09 

1452 


c 


CYCLE 4« parameter INCREMENTS 

H0N09 

1453 


c 



H0N09 

1454 



940 

CONTINUE 

H0N09 

1455 




IF(NOEST) GO TO 943 

H0N09 

1456 

2^5 

c 



H0N09 

1457 


c 


NR INCREMENTS 

M0N09 

1450 




R1«GSQL(1) 4- GSQLOU) 

H0N09 

1459 




R5«G$QL(5) ^ G$QL0(2) 

H0N09 

1460 




R8»GSOL(8) ♦ GSQL0(3) 

H0N09 

1461 

250 



RlO»GS0La0)4'G$QLO(4) 

H0N09 

1462 


c 


COFACTORS FOR EXPLICIT INVERSE 

H0N09 

1463 




R61099-R89R10-R9PR9 

H0N09 

1464 




R61097-R6PR10-R9PR7 

H0N09 

1465 




R6978* R6<»R9-R7*R8 

H0N09 

1466 

255 



R3764-R39R7-R6PR4 

M0N09 

1467 




R2754«R2PR7-R5«R4 

H0N09 

1468 




R2653«R2*R6-R5*R3 

H0N09 

1469 




R2X047«R2^RI0 -R4PR7 

H0N09 

1470 




R2937»R2*R9 -R3*R7 

H0N09 

1471 

260 



R2836*R2*R8-R3«R6 

H0N09 

1472 


c 


INVERSE 

H0N09 

1473 




Cl- R5*R01O99 - R6«R61097 ♦ R70R6978 

H0N09 

1474 




C2 — (R2PR81099 - R3PR61097 * R4PR6978) 

H0N09 

1475 




C3-R7PR3764 - R9*R2754 ♦ R10*R2653 

H0N09 

1476 

265 



C4—(R6PR3764 - R8*R2754 ♦ R9PR2653I 

H0N09 

1477 




C5-R1PR81099 - R3P(R3»R10 - R4<«R9) ♦ R4PCR3»R9 - R4PR8) 

HUN09 

1478 




C6«-(R14R61097 - R3PR21047 ♦ R4P(R2«R9 - R4*R6)) 

H0N09 

1479 




C7-R1PR6976 - R34R2937 ♦ R4*R2836 

H0N09 

1480 




C8«RX'MR5*R10 - R7PR7I -R2<‘R21047 ♦ R4«R2754 

H0N09 

1481 

270 



C9 — (P1P(R5PR9 - R6PR7) - R2«R2937 ♦ R49R2653) 

H0N09 

1462 




CIO»RX*IR5*R0 - R6PR6I - R2PR2836 ♦ R3*R2653 

HON09 

1483 




OET- RX^CX ♦ R2*C2 R3*C3 ♦ R4PC4 

H0N09 

1464 




IF(0ET«EQ.0«0) STOP 2X 

H0N09 

1485 


c 


INCREMENTS 

H0N09 

I486 

Zlb 



DZPdl- -(C1P6L(1)K2P6L(2IK39GU3)«C49GLC4M/0ET 

H0N09 

1487 




0ZP(2)- -(C2«6L(X)K5PGL(2)^C6«GL(3UC7PGL(4n/0ET 

H0N09 

1488 




DZP(3)« '*(C39GLa)K6PGL(2KC8«GU3)K9«GL(4H/0ET 

H0N09 

1489 




DZP<4)- -IC4«GL(X>K7«GL(2HC9«GL(3)^Cl0PGL(4n/DET 

HO NO 9 

1490 




UX427)»SQRT(CX/OET) 

H0N09 

1491 

260 


940 

MODE • 5 

H0N09 

1492 




RETURN 

H0N09 

1493 


c 



H0NO9 

1494 


c 


CYCLE 5« parameter UPDATES 

H0N09 

1495 


c 



HQN09 

1496 

265 


950 

CONTINUE 

H0N09 

1497 


c 


UPDATE ZP 

H0N09 

1498 




ZPX- ZP(J$#1) > DZPUI 

H0N09 

1499 




ZP2- ZPU$^2) ^ DZP(2) 

H0N09 

1500 




ZP3- ZP(JS#3I ♦ DZPC3I 

H0N09 

1501 

290 



ZP4- ZP(JS#4) ^ OZP44) 

H0N09 

1502 




IF(ZPX*GT«ZPIMAX) ZPX-ZPINAX 

H0N09 

1503 




IF(ZP2«GT«ZP2MAX) ZP2-ZP2MAX 

H0N09 

1504 




IF(ZP3*GT«ZP3MAX) ZP3-ZP3MAX 

H0N09 

1505 




1F(ZP4*GT.ZP4MAX) ZP4-ZP4MAX 

H0N09 

1506 


Figure 12. -Continued. 
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295 


300 


lF(ZPl«LT«ZPiniN) ZPl^ZPiniN 
IF( ZP2»lT»ZP2NI<4l ZP2-ZP2NIM 
IF (ZP3«LT*ZP3HXN) ZP3-ZP3MIN 
1F(ZP4*LT»ZPAHXH) ZP<i«ZP<»MIM 
C OTHER l/ARXABLES 
MODE* 1 
RETURN 
C 

END 


H0N09 

H0N09 

HON09 

H0N09 

H0NO9 

H0N09 

H0N09 

H0N09 

H0N09 


SYMBOLIC REFEREHCE MAP (R*l) 

ENTRY POINTS 
X PCMLE 


VARIABLES SN TYPE RELOCATION 


5 

ANS 

REAL 


IPIC 

3 

CHN6CH 

INTEGER 


IPIC 

<» 

CHSIG 

INTEGER 


IPIC 

1362 

Cl 

REAL 



3 

CIHP 

REAL 


SU6RH 

1366 

CIO 

REAL 



1367 

C2 

REAL 



4 

C2HP 

REAL 


SUBRH 

1370 

C3 

REAL 



1371 

C4 

REAL 



1364 

C5 

REAL 



1372 

C6 

REAL 



1373 

C7 

REAL 



1365 

C8 

REAL 



1374 

C9 

REAL 



315 

D 

REAL 

ARRAY 

OAT 

1363 

DET 

REAL 



1 

OHP 

REAL 

, 

SUBRH 

114 

OS 

REAL 

ARRAY 

SENSP 

0 

DTHP 

REAL 


SUBRH 

25 

DTPRIN 

REAL 


IPIC 

0 

DTU 

REAL 


UTEST 

3 

OUT 

REAL 


UTEST 

1426 

DZP 

REAL 

ARRAY 


1375 

ETAl 

REAL 



1376 

ETA2 

REAL 



473 

EIP 

PEAL 


OAT 

474 

E2P 

REAL 


OAT 

137 

F 

REAL 

ARRAY 

OAT 

10 

FSS 

LOGICAL 

ARRAY 

REALTIM 

14 

6AMUT0 

REAL 


UTEST 

1004 

G£ 

REAL 

ARRAY 

SENSP 

1040 

GK 

REAL 

ARRAY 

SENSP 

1010 

GL 

REAL 

ARRAY 

SENSP 

1014 

GSOE 

REAL 

ARRAY 

SENSP 

1026 

GSQL 

REAL 

ARRAY 

SENSP 

62 

GSQLO 

REAL 

ARRAY 

IPIC 

1360 

I 

INTEGER 



1 

1NTA8LS 

INTEGER 


REALTIM 

7 

ITYPE 

INTEGER 


REALTIM 

2 

IVARBUF 

INTEGER 

ARRAY 

REALTIM 

2 

JS 

INTEGER 


IPIC 

1404 

JSTEMP 

INTEGER 



67 

LX 

LOGICAL 

ARRAY 

VAROAT 

1 

MFUE 

INTEGER 


DEVICE 

1356 

HLE 

LOGICAL 



7 

MOCE 

INTEGER 


IPIC 

0 

NC 

INTEGER 


IPIC 

1353 

NOCHC 

LOGICAL 



1354 

NOEST 

LOGICAL 



1 

NP 

INTEGER 


IPIC 

0 

NREAO 

INTEGER 


DEVICE 

1355 

PRINT 

LOGICAL 



105 

PRTIME 

REAL 


IPIC 

21 

RTJC 

REAL 


IPIC 

17 

RTJG 

REAL 


IPIC 

16 

RTJS 

REAL 


IPIC 

23 

rtJz 

REAL 


IPIC 

1411 

R1 

REAL 



1414 

RIO 

REAL 



1027 

R2 

REAL 


SENSP 

1423 

R21047 

REAL 



1422 

R2653 

REAL 



1421 

R2754 

REAL 



1425 

R2836 

REAL 



1424 

R2937 

REAL 



1030 

R3 

REAL 


SENSP 

1420 

R3764 

REAL 



1031 

R4 

REAL 


SENSP 

1412 

R5 

REAL 



1033 

R6 

REAL 


SENSP 

1416 

R61097 

REAL 



1417 

R6978 

REAL 



1034 

R7 

REAL 


SENSP 

1413 

R8 

REAL 



1415 

R61099 

REAL 



1036 

R9 

REAL 


SENSP 

1 

SEED 

REAL 


UTEST 

475 

SIGSO 

REAL 


OAT 

4 

SIGUT 

REAL 


UTEST 

10 

SXGUTO 

REAL 


IPIC 

1432 

ss 

REAL 

ARRAY 


1403 

SI 

REAL 



1406 

S2 

REAL 



1407 

S3 

REAL 



1410 

S4 

REAL 



22 

THRTJC 

REAL 


IPIC 

20 

THRTJ6 

REAL 


IPIC 

24 

THRTJZ 

REAL 


IPIC 

15 

TIME 

REAL 


IPIC 

101 

T4 

REAL 

ARRAY 

DAT 

1405 

TJE 

REAL 



1402 

T4MIN 

REAL 



106 

TL 

REAL 

ARRAY 

DAT 

1400 

7\.MAX 

REAL 



1401 

TLMIN 

REAL 



1357 

TMINCH 

REAL 



6 

TPRINT 

REAL 


IPIC 

11 

TZWUT 

REAL 


UTEST 

71 

U 

REAL 

ARRAY 

IPIC 

3 

UP 

REAL 

ARRAY 

HEAS 

1377 

UT 

REAL 



7 

UTMAX 

PEAL 


UTEST 

12 

UTl 

REAL 


UTEST 

5 

UTIO 

REAL 


UTEST 

13 

UT2 

REAL 


UTEST 

6 

UT20 

REAL 


UTEST 

0 

UX 

REAL 

ARRAY 

VaROAT 

0 

varchkg 

REAL 


REALTIM 

2 

WHP 

REAL 


SUBRH 

2 

WUT 

REAL 


UTEST 

10 

W2UT 

REAL 


UTEST 

0 

X 

REAL 

ARRAY 

OAT 

101 

XD8 

REAL 

ARRAY 

IPIC 

103 

XOC 

REAL 

ARRAY 

IPIC 

75 

XQ 

REAL 

ARRAY 

IPIC 

0 

XS 

REAL 

ARRAY 

SENSP 

77 

XTH 

REAL 

ARRAY 

IPIC 


Figure 12/-Continued. 
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1508 

1509 

1510 

1511 

1512 

1513 
I5l<» 
1515 


50 


Appendix D 


73 

xvz 

REAL 

ARRAY 

IPIC 

0 

YP 

REAL 

ARRAY 

MEAS 

11 

ZPl 

REAL 


IPIC 

32 

2P1MIN 

REAL 


IPIC 

27 

ZP2MAX 

REAL 


IPIC 

13 

ZP3 

REAL 


IPIC 

34 

ZP3M1N 

REAL 


IPIC 

31 

ZP4MAX 

REAL 


IPIC 

1361 

ZIMIN 

REAL 




66 

Y 

REAL 

ARRAY 

IPIC 

36 

ZP 

REAL 

ARRAY 

IPIC 

26 

ZPIMAX 

REAL 


IPIC 

12 

ZP2 

REAL 


IPIC 

33 

ZP2MIN 

REAL 


IPIC 

30 

ZP3MAX 

REAL 


IPIC 

14 

2P4 

REAL 


IPIC 

35 

ZP4MIN 

REAL 


IPIC 


EXTERNALS 

TYPE 

ARCS 

ACUM 


1 

FH 


3 

SENS 


2 

TSI6 


3 


EXP REAL 1 LIBRARY 

FILT 1 

SORT REAL X LIBRARY 


STATEMENT LABELS 


1246 

606 

FMT 

NO 

REFS 


12 50 

630 

FMT 

1260 

631 

FMT 

1264 

632 

FMT 




1270 

633 

FMT 

1274 

634 

FMT 

1276 

636 

FMT 




1302 

637 

FMT 

1305 

638 

FMT 

1310 

639 

FMT 




1314 

640 

FMT 

1317 

642 

FMT 

1323 

64 3 

FMT 




0 

900 

INACTIVE 

352 

910 


377 

911 





403 

920 


417 

925 


427 

92 6 





0 

927 


462 

930 


562 

934 





564 

939 


566 

940 


763 

94B 





765 

950 


244 

1090 


0 

1092 


INACTIVE 


0 

1094 

INACTIVE 

0 

1099 


250 

1100 





303 

1206 


305 

1207 


307 

120B 





311 

1209 


313 

1210 


326 

1301 





330 

1302 


332 

1303 


334 

1304 





503 

9301 


0 

9302 


0 

9303 











.OOPS 

LABEL 

INDEX 

FROM-TO 

LENGTH 

PROPERTIES 




105 


P I 


72 

72 

IIB 


EXT REFS 




124 


♦ I 


74 

74 

7B 


EXT REFS 




137 


♦ I 


75 

75 

7B 


EXT REFS 




152 


♦ I 


76 

76 

7B 


EXT REFS 




165 


♦ I 


77 

77 

7B 


EXT REFS 




363 

911 

I 


161 

169 

15B 

OPT 





430 

927 

* I 


188 

194 

27B 


EXT REFS 




470 

9303 

I 


204 

206 

3B 

INSTACK 





476 

9302 

I 


207 

210 

3B 

INSTACK 






COMMON BLOCKS LENGTH 
VARDAT 70 

DAT 316 


COMMON BLOCKS LENGTH 


SENSP 784 
IPIC 70 
ME AS 3 
SUBRH & 

UTEST 13 
DEVICE 2 
REALTIM 24 

STATISTICS 

PROGRAM LENGTH 14413 801 


CM LABELED COMMON LENGTH 24UB 1291 
73000B CM USED 


Figure 12, -Concluded, 
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SUBR'JtTINc FUT 73//<» f;PT»l 


FTN 4».6^439 78/02/1^. 


1 


1C 


lb 


2G 


2£ 


30 


3i 


SU8K0UTINC FUT(K) H0N09 

C0MH0N/0AT/X(5# 19) > F ( 22 ) * 0( 2 2 ) > 6 IP > €2P# SIGSO H0N09 

COMMON/MEAS/rP(3)#UP(2) H0M09 

C K-CHANNH tNCcX H0N09 

C SAVt X H0N09 

X(K»i)«X(K#5 ) H0N09 

X(K,2 l-XlK^b) H0N09 

XC<,3)-X(K»7) H0M09 

X(K»<»)»X(K»« ) H0N09 

C residuals H0N09 

X(R>9)- rP(i) - X(K>2) rt0N09 

X (K» XOl -YP(2) - X(K,3) H0N09 

X(K,11)«YP(3) - X(K,A) H0N09 

C UPOATc X M0M09 

X(K#5)«0(K, 1)*X(K,5) ♦ D(K# 2)tX(K#6) ♦ 0(K# 3)«X(K#7) ♦ H0N09 

1 0(K» A)4X(K>6) + 0(K>17)«UP<1) 4 0<K#X8)4UP(2) ♦ H0NO9 

2 F(K,i)+X(K#9) 4 F(K,5)4X(K,10) 4 F (K# 9)*X(K# 11 ) H0N09 

X(K,6)-D(K^ 5)4X(K#5) 4 D(K/ 6)4X(K#6) 4 0(K# 7)*X(K#7) 4 H0N09 

1 0(K> 9)4X(K>8) 4 0(K>19)4UP(1) 4 0 (K>20)4UP ( 2) 4 H0N09 

2 F<K#2 )*X(K#9) ♦ F(K>6)4X(K#10) 4 F ( K# lO.) 4X ( K» IX ) 40N09 

<<K>7)»0(K, 9)4X(K^5) 4 DIK, 10)4X(K#6) 4 0(K# XI )4X( K> 7) 4 H0N09 

1 0(K#12)4X(K»8) 4 0(K,21)*UP(1) 4 0 ( K#22) 4UP ( 2 ) 4 M0N09 

2 F(K> 3)4X IK#9) 4 F(K>7)4X(K#10) 4 F (K#X1 I 4X ( Krll ) M0N09 

X(K»9)«0(K#13)4X(K#5) 4 0( K# l<i) *X ( K# 6 ) 4 D< K# X9 )4Xi K#7 ) 4 HQN09 

1 0(X»16)4X<K>8) ♦ H0N09 

2 F<K» A»4X(Kf9) ♦ F(K>9)4X(K#10) 4 F (K# 12 )4X (K » 1 1 ) M0NO9 

C (NUKRIHNU) H0N09 

X (Kf 12)«F(K#lA)4X(Kf 9)4X(K»9) 4 F(K» X8 ) ♦X(K#10)*X<K# 10 ) 4 H0N09 

1 F(K,22)4X(K#ll)4X<K,ll) 4 ( F CK# 15 ) 4F ( K# 17)) *X ( K# 9 ) 4X < K, 10 ) 4 H0N09 

2 (F(K#l6)4F(K#20))4X(K,9)4X(Kai) 4 ( F ( K# 19) 4F (K# 21) )*X ( K# 10) * H0N09 

3 XCK^ll) H0N09 

C SUH H0M09 

X(K,l^)«tlP4X(K#.lA) 4 X(K#12) H0M09 

C (NUMRIHnU) 4 LNO£TR H0N09 

X(k#l3)-X(K> 12) 4 3IGSQ4F<K,13) rt0N09 

RETURN H0N09 

END 


SYMBJLIC REFERENCL "AP (R-1) 


ENTRY POINTS 
3 FUT 


VARIA 

"LES 

SN TYPE 

RELOCATION 

315 

0 

REAL 

ARRAY 

DAT 

474 

l2P 

REAL 


FAT 

C 

K 

INTEGER 


F*P. 

3 

UP 

REAL 

array 

PEAS 

0 

YP 

REAL 

ARRAY 

MLAS 


473 

EIP 

REAL 


OAT 

137 

F 

REAL 

ARRAY 

DAT 

475 

SIGSQ 

REAL 


OAT 

0 

X 

REAL 

ARRAY 

DAT 


common blocks length 
JAT 31B 

MFAS 5 


statistics 

PRCGRAM length I72fl 122 

CM LmTELcD common length 503B 323 


Figure 13, -Program listing of subroutine FILT, 


08«28r50 


589 

590 

591 

592 

593 
59A 

595 

596 

597 

598 

599 

600 
601 
602 
603 
60A 

605 

606 

607 

608 

609 

610 
611 
612 

613 

614 

615 

616 

617 

618 

619 

620 
621 
622 

623 

624 

625 
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SUBHUfTINfc SENS 73/7^ 0PT«1 FTN 79/02/1^* 08«26«90 


1 


SUBROUTINE SENS(K#J» 

H0N09 

1679 



C0MM0N/0AT/X(5,19)# F( 5# 22)# D( 5#22 )#€1P#E2P#SI6SQ 

H0N09 

1680 



COMMGN/McAS/YP(3)#UPC2) 

HON09 

1681 



C0MM0N/SEN$P/XS(<»#19)#DS( 5#<i#22)#GE(^)#GU^)»GSQ£U0)#GSQL(10) 

H0N09 

1682 

5 


1 #6K(5#^#12) 

HON09 

1683 


C 

J-PARANETcR INDtX 

H0N09 

168A 


C 

K-CHANNEL INDEX 

H0N09 

1685 


C 

SAVE GRAD X 

H0N09 

1686 



XS(J#1)« XSU#5) 

H0N09 

1687 

10 


XSU#2)« XSU»t) 

H0N09 

1688 



XS(J#3)- XS(J#7) 

H0N09 

1689 



XS(J#A)« XS(J»8) 

H0N09 

1690 


C 

GRAD NU* -H*(GX) 

H0N09 

1691 



XS(J, 9)- -X$(J,2) 

H0N09 

1692 

15 


XSIJ# iOl- -XSIJ^S) 

H0N09 

1693 



XS(J»1U« -XSU#<») 

H0N09 

169A 


C 

GRAO X update 

H0N09 

1695 



XS(J#5)«0(K# 1)AXS(J#5) ♦ D(K# 2)*XS(J#6) ♦ 0<K# 3)PXS(J#7) ♦ 

N0N09 

1696 



1 0(K, 4)*XS<J#8) ♦ OS(K#J# 1)4>XIK#5) ♦ OSCK#J# 2)PX(K#6) ♦ 

H0N09 

1697 

20 


2 OS(K#J» 3)*X(K#7) ♦ OS(K#J# 4)PX(K#6) ♦ OS ( K# 0* 17) «UP ( 1 ) 

N0N09 

1698 



3 4 OS(K#J#18)4>UP(2) ♦ P(K# 1M>X$I4#9) t F(K# 5)*XSU#10) ♦ 

H0N09 

1699 



A F(K» 9)PXS(J*il) 4 GK(K#J# 1 )^XCk#9) 4 GK|K#J# 9)*XIK#10) 

H0N09 

1700 



5 ♦ GK(K#J# 9)AX(K#ll) 

H0N09 

1701 



XSU#6)"0(K# 5)PX$(J»5) 4 D(K# 6)PX$U#6) 4 OIK# 7)«XS(J#7) * 

H0N09 

1702 

25 


1 OIK* 8)*»X$U*8) 4 DSIK#J* 5)«XIK*5) ♦ DSIK#J# 6)PXIK*6) ♦ 

N0N09 

1703 



2 0SIK#J# 7)^XIK#7) ♦ DSIK#J» B)*XIK#8) ♦ OSI K# J#19)*UP ID 

H0N09 

170A 



3 * 0SIK*J#?0)«UPI2) 4 FIK# 2)PXS(J#9) 4 FIK# 6)PX$IJ#10) ♦ 

H0N09 

1705 



A FIK#10)AX$I J*ll) ♦ GKIK*J* 2)AXIK*9) ♦ 6KIK#J# 6)«XIK#10) 

H0N09 

1706 



5 4 GKIK#0#10)tXIK#ll) 

H0N09 

1707 

30 


XSU*7)»0(K* 9)AXSU#5) 4 0IK#10)AXSU#6) 4 0 IK# 11 ) «XS I J# 7 ) 4 

H0N09 

1708 



1 0IK#12)«XSIJ#8) 4 0SIK#J# 9)AX(K#5) 4 OS (K# J# 10) PX (K# 6) 4 

H0N09 

1709 



2 DSIK# J#11)4XIK*7) 4 0SIK#J*12)AXIK*8) ♦ DS IK# J # 21 ) *UPI 1 ) 

H0N09 

1710 



3 4 OSIK# J#22)*UPI2) 4 FIK# 3)4XSU#9) 4 FIK# 7)PXSIJ#10) 4 

H0N09 

1711 



A F(K*11)4XSU*11) 4 GKIK*J* 3)4XIK#9) 4 GKIK#J# 7)«XIK*10) 

H0H09 

1712 

35 


5 4 GKIK# J# 11)4XIK#11) 

H0N09 

1713 



XStJ#6)-0IK# 13)4X$IJ#5) 4 DIK#1A)«XS(J#6) 4 OlK# 19 )«XS 14# 7) 4 

H0N09 

171A 



1 0tK#16)4XSI4#8) 4 0$IK#4#13)4XIK#9) 4 OSt K# J# lA ) 4X IK#6 ) 4 

H0N09 

1715 



2 0SIK#4#15)4XIK#7) 4 D$ I K# 4# 16) 4X I K# 8) 

H0N09 

1716 



3 4 FIK# A)4XS<4#9) 4 FIK# 8)4XSI4#10) 4 

N0N09 

1717 

MO 


A FIK#12)4XSI4#11) 4 GKIK#4# A)4XIK#9) 4 6KIK#4# 6)4XIK#I0) 

H0NO9 

1718 



5 4 6KIK#4#12)4XIK#11) 

H0N09 

1719 



RETURN 

H0N09 

1720 



END 

H0N09 

1721 



SYMBOLIC 

REFERENCE 

MAP <R«1) 


FNTRV 

POINTS 




3 

SENS 




VARIABLES SN 

TYPE 

RELOCATION 

315 

0 

REAL 

ARRAY 

OAT 

A73 

clP 

REAL 


OAT 

137 

F 

REAL 

ARRAY 

DAT 

lOAO 

GK 

REAL 

ARRAY 

SENSP 

101 A 

GSQE 

REAL 

ARRAY 

SENSP 

0 

J 

INTEGER 


F.P, 

A75 

SIGSQ 

REAL 


DAT 

0 

X 

REAL 

ARRAY 

OAT 

0 

YP 

REAL 

ARRAY 

MFAS 

COMMON 

BLOCKS 

LENGTH 




OAT 

318 




MEAS 

5 




SENSP 

78A 




114 

OS 

REAL 

ARRAY 

SENSP 

474 

E2P 

REAL 


OAT 

1004 

GE 

REAL 

ARRAY 

SENSP 

1010 

GL 

REAL 

ARRAY 

SENSP 

1026 

GSOL 

REAL 

ARRAY 

SENSP 

0 

K 

INTEGER 


F.P. 

3 

UP 

REAL 

ARRAY 

MEAS 

0 

X$ 

REAL 

ARRAY 

SENSP 


STATISTICS 

PROGRAM length 
CM LA3ELE0 COHMTN LENGTH 


^026 258 

2129B 1107 


Figure 14. -Program listing of subroutine SENS. 


53 



Appendix D 


SUPRJOTINf 4C'Jh 0PT"l FTN 


1 SUBROUTINE ACUh<K) 

CGHHaN/0AT/X(5>19), F(5>22)f 0(5# 22)#EIP#E2P«SI6$Q 

CCir-K0N/SeNSP/XS(4,19)#DS( 5#4#22)*GE(4)#6L(4)#GSOEaO)#G$QLaO> 
1 #GK(5#A#12) 

5 C K-CHANNEt INDICATOR 

Tll»XSa#9)^F(K#l4) ♦ XS(1#10)«P(K#17) ♦ X$(1#UI«P(K#20) 
Ti2«XS(l#9)*F(K#15) ♦ XS(1# 10 )«F(K# 18 ) ♦ X$(l# intf(K#2U 
m«XS(l#9)«F(K#16) 4 XS(1#10)«F(K#19I * XSa#Xll«F fK#22 ) . 
T21«XS(2#9)4F(K#UI ♦ XS(2#10)«F(K#17) ♦ X$(2»U)*F(K»20) 

1C T22«XS(2#9)*F(K#15) ♦ XS(2# 10 )*F(K# 18) ♦ XS (2 # 11 ) (K# 21 ) 

r23-XS(2#9)AF(K#16) ♦ XS(2# 10 )^F(K# 19) 4 XS (2# 11 )« F (K# 22 ) 
T3l-XS( 3#9)4F(K,14) ♦ XS(3# 10)«F (K# 17) ♦ XS ( 3#1 1 ) «F (K#20 ) 
T32«XS(3»9)4F(K#15) ♦ XS (3 # 10 )«F (K# 18 ) ♦ X$ (3 # 11 )<»F (K# 21 ) 
r33«XS(3#9)4F(K#16) 4 XS(3#10)«F(K#19) 4 XS(3#11)«F(K#22) 

It* T4l«XS(4#9)4F(K#l<») 4 XS(4#10)«F(K# 17) 4 XS (4#11) 4F (K#20) 

T*»2«XS(‘i#9)4F(K#l5) 4 XS (^ # 10 )4F ( K# 18 ) 4 XS (^*11 ) «F (K# 21 ) 
TA3«XS(<»>9)4F(K#16) 4 XS(4#10)4F(K#19) 4 XS ( ^# 11 )♦ F (K^ 22 ) 

C 

C $«(GRAO NU)(RI)(NU) 

20 C 

S1«T114X(K#9) 4 T124X(K#10) 4 T134X(K#11) 

S2*r?l^X(K#9) 4 T224X(K#10) 4 T234X<K>H) 

$3«T314X(K,9) 4 T32*X(K#10) 4 T334X(K#11) 

54»TA14X(K#9) 4 T424X(K#10) 4 T434X(K#11) 

25 GE(1)«E2P*(GE(1) - SI) 4 $1 

GE(2)«E2P^(Gfc(2) - S2 ) 4 S2 
GE (3)"£2P4(GE(3) • S3) 4 S3 
6E(A)»E2P4(GE(4) - S^) 4 S4 
C 

30 GL(i)» E1P4GL(1I 4 GEU) 

GL<2)- E1P4GL(2) 4 GE(2) 

GL(3)- E1P4GL(3) 4 GE(3) 

GL(A)« E1P4GL(^) 4 GE(^) 

C 

35 C S - (GRAD NU)(RI)(GRAD NU) 

r 

Sll-T114XS( lf9) 4 T12^XS(1#10) 4 T134XS(1#11) 

S12»T114XS(2#9) 4 T12*XS(2#10) 4 T134XS(2#11) 

S13«T11*XS(3#9) 4 T124XS(3#10) 4 T13«X$(3#11) 

40 S14«T114XS(4#9) 4 T12*XS(4#10) 4 T134XS(4#11) 

S22*T214XS(2#9) 4 T224XS(2#lO) 4 T234XS(2#11) 

S23-T214XS(3»9) 4 T22<>XS (3 # 10 ) 4 T23«XS(3#11) 

S24«T214XS(4#9) ♦ T224X$(4#10) 4 T230XS ( 4# 11) 

S33-T31»XSn#9) 4 T 324XS(3#10) 4 T334XS(3#11) 

45 S34.T314XS(4#9) ♦ T324X$(<i# 10) ♦ T334XS(4#ll) 

T<il4XS(4, 9) 4 T424XS(4#10) 4 T434XS(4#11) 

GS0E(l)«t2P*(&S06(l) - Sll) 4 Sll 
&SQF(2)-F2P*(GSQ6C2) - $12) 4 $12 
GSQt (3)«E2?4(GS0e(3) • Sl3) 4 Sl3 
5C GSCE (4)*c2P4(6S0n(4) - Sl4) 4 $14 

GSQt(5)*E2P4(GSOE(5) -- S2?) 4 $22 
GS0E{6)»E2P4(6SQ£(6) - S23) 4, $23 
GS?e (7)*E2P4(6SCE(7) - S24) 4 S2^j 
GS3E(8)-t2P*(GS06(8) - S33) 4 $33 
55 G$06 (9)«E2P* (GS0F(9) - $34) 4 $34 

GSOE (10)*E2P4(&SQe( 10) • S44) 4 $<»<» 

C 

GSOL(l) • flP4GS0l(l) 4 6$0E(1) 

GSQL<2) - eiP4GS0L(2) 4 GSQE(2) 

6C GS0L(3) • 61P4GSQL(3) 4 GSQfcO) 

GS0LC4) • eiP»6$Qt(4) 4 GSQE(4) 

GS9L(5) - E1P4GS0L(5) 4 GSQE(5) 

GS0L(6) • elP4GS0L(6) 4 G$0E(6) 

GS0L(7) • EIP4GSOU7) 4 GS0E(7) 

45 GSOL(B) • E1P4G$QL(8) 4 GS0E(8) 

G$CL(9) • tlP*6SQL(9) 4 G$0E(9) 

GSOLUO)* ElP4GSQt(lO)4 GSOEdO) 

RETURN 

END 


78/02/14* 08*28*50 


H0N09 

81 

H0N09 

82 

H0N09 

63 

H0N09 

84 

H0N09 

65 

H0H09 

86 

H0N09 

07 

H0N09 

68 

H0N09 

69 

H0N09 

90 

H0N09 

91 

H0N09 

92 

H0N09 

93 

H0N09 

94 

H0N09 

95 

HON 09 

96 

HON 09 

97 

H0N09 

96 

H0H09 

99 

H0N09 

100 

H0N09 

101 

H0N09 

102 

H0N09 

103 

H0N09 

104 

H0N09 

105 

H0N09 

106 

H0N09 

107 

H0N09 

108 

H0N09 

109 

H0N09 

1X0 

H0N09 

111 

H0N09 

112 

H0N09 

113 

H0N09 

114 

H0N09 

115 

H0N09 

116 

H0N09 

117 

H0N09 

116 

H0H09 

119 

H0N09 

120 

H0N09 

121 

H0N09 

122 

H0N09 

123 

H0N09 

124 

H0NO9 

125 

H0N09 

126 

HON 09 

127 

H0N09 

128 

H0N09 

129 

H0N09 

130 

H0N09 

131 

H0N09 

132 

H0N09 

133 

H0N09 

134 

HON09 

135 

H0N09 

136 

H0N09 

137 

H0N09 

136 

H0N09 

139 

HON09 

140 

H0N09 

141 

H0N09 

142 

HON09 

143 

H0N09 

144 

H0N09 

145 

H0N09 

146 

H0N09 

147 

H0N09 

146 

H0N09 

149 


SYMBCLIC REFERtNCe MAP (R«l) 

fNTRY POINTS 
?. A CUM 

Figure 15. -Program listing of subroutine ACUM, 
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VAMABlHi SN TYPf «E10CATI0N 


316 

3 

REAL 

ARRAY 

DAT 

114 

OS 

REAL 

ARRAY 

SENSP 

473 

EIP 

REAL 


CAT 

474 

E2P 

REAL 


OAT 

137 

F 

REAL 

ARRAY 

DAT 

1004 

GE 

REAL 

ARRAY 

SENSP 

1C4(. 

OK 

REAL 

ARRAY 

SENSP 

1010 

GL 

REAL 

ARRAY 

SENSP 

1C14 

.SjE 

REAL 

ARRAY 

SENSP 

1C26 

GSQL 

REAL 

ARRAY 

SENSP 

0 

K 

INTEGER 


F.P* 

475 

SIGSC 

REAL 


DAT 

344 

Si 

REAL 



350 

Sll 

REAL 



2'jX 

S12 

REAL 



352 

$13 

REAL 



363 

514 

REAL 



345 

S2 

REAL 



3?4 

S22 

REAL 



355 

$23 

REAL 




S24 

REAL 



346 

S3 

REAL 



^>7 

S33 

REAL 



360 

S34 

REAL 



347 

04 

PEAL 



361 

$44 

REAL 



33C 

ni 

REAL 



331 

T12 

REAL 



332 

T13 

REAL 



333 

T2l 

REAL 



334 

T22 

REAL 



335 

T23 

REAL 



336 

T31 

REAL 



337 

T32 

REAL 



340 

T33 

real 



341 

T41 

REAL 



342 

T42 

REAL 



343 

T43 

REAL 



0 

X 

REAL 

ARRAY 

OAT 

0 

XS 

REAL 

ARRAY 

SENSP 


COMMON BLOCKS LENGTH 
OAT 31P 

JENSP 


STATISTICS 

FPCGfrAM LENGTH 362B 2^2 

CM labeled CUHMJN LENGTH 21169 1102 


Figure 15. -Concluded. 


SUBROiJTINC FH 

73/74 0PT*1 

FTN 4.64439 

78/02/X4. 

Oe.28.50 

1 

SUBROUTINE Fh<U>Y,X) 


H0N09 

567 


COhMON/SUBRH/ DT »D#W#C1#C2 


H0N09 

568 


DIMENSION X<2) 


HON09 

569 

C 

FILTER 


H0N09 

970 

5 C 

Y/U- (S*SI/(S*S ♦ 2*0*W4S ♦ «♦«) 


H0NQ9 

571 


2 Y«U - X(2»*C2 - XCUACl 


H0N09 

572 


X(l)»xm ♦ 0T*XC2) 


H0N09 

573 


Xf2)«Xf2) ♦ OT*Y 


H0N09 

574 


RETURN 


H0N09 

575 

10 

END 


H0N09 

576 


SYMBOLIC REFERENCt MAP (R«ll 

ENTRY POINTS 
3 FH 

VARIABLES SN TYPE RELOCATION 


3 

Cl 

REAL 


SUBRri 

4 

C2 

REAL 

SUBRH 

1 

0 

REAL 


SUBRH 

0 

DT 

REAL 

SUBRH 

0 

u 

REAL 


F.P* 

2 

W 

REAL 

SUBRH 

0 

X 

REAL 

ARRAY 

F.P. 

0 

Y 

REAL 

F.P. 


STATEMENT LABELS 

C 2 INACTIVE 

COMMON BLOCKS LENGTH 
SU8RH t 

STATISTICS 

PROGRAM LENGTH 160 14 

CM LABELED COMMON lENGTH 6B 5 


Figure 16. -Program listing of subroutine FH. 
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SObhr.jTiNt rSIG 7<»/7<i QPT»1 


FTN 4.6+439 77/12/14. 18.57.02 


SUHEOUTINE TSIG<E7A1.LTA2#UT) 

H0N09 

1720 

UT/TnPUT • S^SIGUT/<S*S ♦ 2*DtT*W0T*S ♦ WUT*WUT) 

H0N09 

1721 

INPUT • RANDOM NO >»1TH C MEAN UMFOPM -M.? TO ).5 

H0NC9 

1722 

UT1.UT2 STATES IN SECOND ORDER EUTEft 

H0N09 

1723 

CrMMaN/UTeST/DT#SIE0>WUT,DUT#SI6UT»lJTlM»UT20#UTMAX,W2UT>TZWUT> 

H0N09 

1724 

1 LT1.UT2#6AMUT0 

H0N09 

1725 

2 SE£C-AMDD(312t.4SlE0# 34359738368. ) 

H0N09 

1726 

S«U.291C38304567g-104SE60 - C.5 

H0N09 

1727 

S-W2UT4UT1 ♦ TZWUT+UT2 ♦ GAMUTO*S IGUT+S 

HON09 

1728 

UTl-UTl ♦ 0T+UI2 

H0N09 

1729 

UT2-UT2 ♦ S 

H0N09 

1730 

UT-UT2 

H0NO9 

1731 

IH ABS(UT).66.UTMAX) UT«S IGN(UTMAX.UT) 

H0N09 

1732 

St ED-AMOIX 3125. ♦SLEO> 34359738368.) 

H0N09 

1733 

S»O.291t-333C4567E-10*SbeD - 0.5 

H0N09 

1734 

ETAl-S 

HQN09 

1735 

S ttn-AMuD( 3125. *S£ED# 34359738368.) 

H0N09 , 

1736 

S*o.29103d3O4567=-104SEE0 - C.5 

H0N09 

1737 

ETA2-S ' 

H0N09 

1738 

RETURN • 

H0N09 

1739 

END 

H0N09 

1740 


SYMBOLIC REFERqNCt MAP (R«l) 

ENTRY POINTS 
3 TSIG 


VARIABLES 

SN TYPc 

RELOCATION 

0 

DT 

PEAL 

UTEST 

0 

ETAi 

real 

F.P. 

14 

6AMUT0 

REAL 

UTcST 

1 

SEED 

REAL 

UTcST 

a 

TZWUT 

REAL 

UTEST 

7 

utmax 

real 

UTEST 

5 

UTlw 

RcAL 

UTEST 

t 

UT20 

REAL 

UTEST 

IG 

I.2UT 

REAL 

UTcST 


INLIM fUNCTICNS 

1 YPL 

ARGS 


ABS 

R:tAL 

i 

INTRIN 

SIGN 

REAL 

2 

INTRIN 


statement labels 

0 2 inactive 


3 

OUT 

REAL 

UTEST 

0 

ETA? 

R^AL 

P*P* 

46 

S 

PEAL 


4 

STGUT 

REAL 

UTEST 

i) 

UT 

REAL 

F.P. 

12 

UT1 

REAL 

UTEST 

13 

UT? 

REAL 

UTEST 

2 

WUT 

REAL 

UTEST 


AMHO PEAL 2 INTRIN 


COMMON BLOCKS LE>IGTH 
U1 EST U 


statistics 

peCGPA« It NOTH 
CM LABELED COMMON LENGTH 


47B 39 

15B 13 


Figure 17, -Program listing of subroutine TSIG. 
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APPENDIX E 


FLOW CHART FOR INITIALIZATION PORTION 
OF PCMLE SOFTWARE 


This appendix consists of the flowcharts for subroutine NRTIC (Fig. 18) 
and MODEL (Fig. 19). 



Figure 18. -Flowchart for subroutine NRTIC. 
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Appendix E 



Figure 19. -Flowchart for subroutine MODEL. 
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APPENDIX F 

LISTING OF PCMLE INITIALIZATION SOFTWARE 


This appendix contains the program listings of subroutines NRTIC, 
DISC, DIAK, CAL, MP, FHIC, and INPT, as presented in Figures 
through 27, respectively. 


SUBROUTINE NRTIC 73/7<i ORT«l 


FTM R.6*<>52 78/0'./06. 


1 


b 


10 


15 


20 


25 


30 


35 


^0 


<»5 


C 


C 


SUBkOUTlNE NRTIC H0N09 

COHMON/PRSWTCH/Ii>SCH H0N09 

C0HH0N/CI/XCl#XC2#XC3#XC^ ri0N09 

COHMON/OFILE/HF rtON09 

C0HH0N/VARDAT/UX<55)#LXC15) H0N09 

COMMUN/S£NSP/XS( A#19)#DS< 5# 22 ) > GE ( A GL ( A )# GSOE (10) # 6SQL <10 ) H0N09 

1 #6K(5#A#12) H0NO9 

C'JMH0N/0AT/X(5#19># F ( 5# 2 2 ) f D (5 #22 )» fc 1P> c 2P* S I GSQ H0NO9 

COHhOH/lPIC/ NC#NP# J$#CHNGCH#CHS1G#ANS HONOO 

1 #TPRINT#H00E#SI6UT0#2P1#ZP2»ZP3#ZPA,TIME#RT JS H0N09 

2 #RTJG#THRTJG#RTJC#THRTJC#RTJZ#THRTJZ#OTPRIN H0N09 

3 #ZPIHAX#ZP2NAX#ZP3MAX#ZPAMAX H0N09 

A #ZP1HIN#ZP2HIN#ZP3M1N#ZPAH1N HQN09 

5 #ZP( 5»A)#GSQL0(A)#y(3)#U(2) HON09 

6 #XyZ(2)#XQ(2)#XTH(2)#XD0(2)#XDC(2)#PPTIHE H0N09 

C0MH0N/MEAS/yP(3)#UP(2) H0N09 

C0HHUN/SU8R/ NX #NR#NN# NRH#E E # 1TER#0T» OETR H0NO9 

1 #DUHH(292) H0N09 

CUHMON/SU8RH/OTHP#ONP#WHP#CIHP#C2HP riON09 

CCMMUN/UTEST/DTU#S6EO#WUT#OUT#SIGUT#UT10#UT20#UTHAX#W2UT#TZWUT# H0N09 
1 UT1#UT2#GAHUT0 H0N09 

COHHON/OE VICE/NREAO# NFILE N0N09 

COHHON/NOH/THEN# PHIN#UN# VN# WN# UAN# VAN# WAN H0N09 

DIMENSION 6N(2#10) H0N09 

INTEGER CHN6CH#CriSI6 H0N09 

LOGICAL LX H0N09 

N0N09 

CALL RTHQLOl H0N09 

REWIND 10 H0N09 

READ (10) XS#0$»Gc#GL#6SQE#GSQL#GK H0N09 

REAO(IO) X#F#0#E1P#£2P#SIGSQ H0N09 

RE AO (10) NC#NP#JS#CHNGCH#CHSI6#ANS#PRTIM£ HQN09 

1 #TPRINT#Ha0£#SIGUT0#ZPl#ZP2#ZP3#ZPA#TINE#RTJS H0N09 

2 #RTJG#TNRTJG#RTJC#THRTJC#RTJZ»THRTJZ#DTPRIN H0N09 

3 ,ZP1MAX#ZP2HAX#2P3HAX#ZPAMAX H0N09 

A #ZP1MIN#ZP2HIN»ZP3HIN#ZPAMIN N0N09 

5 #ZP>GS0L0#y#U#XVZ#XQ#XTH#X0B#XDC H0N09 

PEAD(IO) 0THP#0HP#W*(P#C1HP#C2HP#0TU#S£6D#WUT#0UT#SISUT HON09 

1 >UT10#UT20#UTMAX#W2UT»TZWUT#UT1#UT2#GAMUT0 ri0N09 

REAO(IO) NREAD#NFILE H0N09 

READ(IO) UX#LX H0NO9 

HON09 

XC1»0. H0N09 

XC2-0. H0N09 

XC3-0* H0N09 

XCA«0* H0N09 

CALL RTSRT H0N09 

RETURN H0N09 

FNO HQN09 


SyMBOLlC REFERENCE MAP (R«l) 


r 

ENTRy 

POINTS 










1 

NRTIC 










VARIABLES 

SN TyP£ 

RELOCATION 








5 

ANS 

RcAL 

IPIC 


3 

CHNGCH 

INTEGER 


IPIC 


A 

CNSIG 

INTEGER 

IPIC 


3 

CIHP 

REAL 


SUBRH 



C2hP 

REAL 

SUBRN 


315 

D 

REAL 

ARRAy 

OAT 


Figure 20, -Program listing of subroutine NRTIC. 


MODEL, 

20 


0b*36»A/ 

1103 

116A 

U65 

1166 

1167 

1168 

1169 

1170 

1171 

1172 

1173 
117A 

1175 

1176 

1177 

1178 

1179 

1180 
1181 
1182 
1183 
118A 
1185 
1136 

1187 

1188 

1189 

1190 

1191 

1192 

1193 
119A 

1195 

1196 

1197 

1198 

1199 

1200 
1201 
1202 
1203 
120A 

1205 

1206 

1207 

1208 

1209 

1210 
1211 
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Appendix F 




7 

DcTR 

PEAL 


ll<» 

OS 

REAL 

ARRAY 

0 

DThP 

RE U 


0 

ort 

REAL 


3 

OUT 

PEU 


A73 

EIP 

PEAL 


137 

F 

PEAL 

ARRAY 

lOOA 

GE 

REAL 

ARRAY 

1010 

6L 

REAL 

ARRAY 

lOlA 

GSOE 

REAL 

ARRAY 

tz 

GSOLO 

PEAL 

ARRAY 

b 

ITkR 

INTtGEP 


67 

LX 

LOGICAL 

ARRAY 

1 

MFILE 

INTEGER 


0 

NC 

INTtGEP 


1 

NP 

INTEGER 


0 

NRL-AO 

INTEGER 


0 

NX 

INTEGER 


105 

PKTIME 

PEAL 


17 

PTJG 

PEAL 


23 

RTJ2 

PEAL 


<.75 

SIGSQ 

PEAL 


10 

SIGUTO 

PcAL 


22 

THFTJC 

PdAL 


2<» 

TrtPTJZ 

PEAL 


6 

TPPINT 

PEAL 


71 

U 

PEAL 

• ARRAY 

2 

UN 

REAL 


7 

UTMAX 

REAL 


5 

UTIO 

REAL 


6 

UT20 

REAL 


6 

VAN 

PEAL 


7 

NAN 

PEAL 


<» 

nN 

PEAL 


10 

W2tT 

REAL 


0 

XCl 

PEAL 


2 

XC3 

REAL 


101 

XOB 

REAL 

ARRAY 

7 5 

XJ 

PEAL 

ARRAY 

77 

XTH 

REAL 

ARRAY 

66 

Y 

REAL 

ARRAY 

36 

2P 

REAL 

ARRAY 

26 

ZPIMAX 

HEAL 


12 

ZP2 

RcAL 


33 

ZP2MIN 

PEAL 


30 

ZP3MAX 

peal 


16 

2P6 

REAL 


35 

ZP6MIN 

REAL 


FILE NAMES 

MODE 



TAFEIO 

UNFMT 


EXTERNALS 

TYPE 

ARGS 


RThOLDl 


0 

COMMON 

BLOCKS 

length 



PRSWTCH 

1 



Cl 

6 



OFILE 

1 



VAKDAT 

70 



SENSP 

786 



OAT 

310 



IPIC 

70 



meas 

5 



SU6R 

300 



SUBPH 

5 



UTEST 

13 



DlVICE 

2 



NOM 

e 


STATISTICS 




PUJOKAi LfcNGTH 

LABkLEO COf^MON LENGTH 
730008 CM USED 


SUBR 

1 

OHP 

SENSP 

6 

DT 

SU8RH 

25 

OTPRIN 

UTEST 

10 

OUHh 

UTcST 

6 

EE 

OAT 

676 

E2P 

OAT 

16 

GAMUTO 

SENSP 

1060 

GK 

SENSP 

167 

GN 

SENSP 

1026 

GSQL 

IPiC 

0 

IPSCH 

SUBK 

2 

4S 

VAROAT 

0 

MF 

DEVICE 

7 

MODE 

IPIC 

2 

NN 

IPIC 

1 

NR 

DEVlCt 

3 

NRM 

SUBR 

1 

PHIN 

IPIC 

21 

RTJC 

IPIC 

16 

RTJS 

IPIC 

1 

SEED 

OAT 

6 

SISUT 

IPIC 

0 

THEN 

IPIC 

20 

THRTJG 

iPlC 

15 

TIME 

IPIC 

11 

TZWUT 

IPIC 

5 

UAN 

NOM 

3 

UP 

UTEST 

12 

UTl 

UTEST 

13 

UT2 

UTtST 

0 

UX 

NOM 

3 

VN 

NOM 

2 

WHP 

NJM 

2 

WUT 

UTEST 

0 

X 

Cl 

1 

XC2 

Cl 

3 

XC6 

IPIC 

103 

xoc 

IPIC 

0 

xs 

IPIC 

73 

xvz 

IPIC 

0 

YP 

IPIC 

11 

ZPl 

IPIC 

32 

ZPIMIN 

IPIC 

27 

ZP2MAX 

IPIC 

13 

ZP3 

IPIC 

36 

ZP3N1N 

IPIC 

IPiC 

31 

ZP6HAX 


RTSRT 


213B 139 

3055B 1581 


REAL 


SUBRH 

REAL 


SUBR 

REAL 


IPIC 

REAL 

ARRAY 

SUBR 

REAL 


SUBR 

REAL 


OAT 

REAL 


UTEST 

REAL 

ARRAY 

SENSP 

REAL 

^UNOEF 


REAL 

ARRAY 

SENSP 

INTEGER 


PRSwTCH 

INTEGER 


IPIC 

INTEGER 


OFILE 

INTEGER 


IPIC 

INTEGER 


SUBR 

INTEGER 


SUBR 

INTEGER 


SUBR 

REAL 


NOM 

REAL 


IPIC 

REAL 


IPIC 

REAL 


UTtST 

REAL 


UTEST 

REAL 


NOrt 

REAL 


IPIC 

real 


IPIC 

REAL 


UTtST 

REAL 


NOM 

REAL 

ARRAY 

McAS 

REAL 


UTEST 

REAL 


UTEST 

REAL 

ARRAY 

VAROAT 

REAL 


NOM 

REAL 


SUBRH 

REAL 


UTEST 

REAL 

ARRAY 

OAT 

REAL 


Cl 

REAL 


Cl 

REAL 

ARRAY 

IPIC 

REAL 

ARRAY 

SENSP 

REAL 

ARRAY 

IPIC 

REAL 

ARRAY 

NEAS 

REAL 


IPIC 

REAL 


IPIC 

REAL 


IPIC 

REAL 


IPIC 

REAL 


IPIC 

REAL 


IPIC 


\ 


0 


Figure 20 


Concluded, 



Appendix P 


S'JPP-IUTIN 

L 

MODEL 7Z/r^ OPT-1 FTN A,64<»39 

7«/02/lA. 

09«26< 

1 


SUBROUTINE npoEL( ZP i / ZP2 » Z Pd » > S 1 GVZ # S 1 GO # S I GTH , 

H0NO9 

835 



1 SICXG#SIGZG#ICH,NTER) 

H0N09 

836 


c 

ICH IS CHANNEL INDICATOR bScO FOR DEFINING STORED DATA 

HON09 

837 



CaMM0N/VAP0AT/bX(50)#LX(50) 

H0N09 

838 



C0MMnN/DAT/X(5^19)» F ( 5#22 P( 5# 22 )# El P# E 2P,SIGS0 

H0N09 

839 



C0MM0N/SENSP/XS(^,19)#0S( 5> » 22 ) » GE ( <» ) > GL ( 4 ) , GSOE <10 ) » 6S0L ( 10) 

^0N09 

8A0 



1 »GK(5#<»fl2) 

H0N09 

8A1 



C0NHaN/SUBRH/0THP#0riP#WHPfClHP#C2HP 

H0N09 

8A2 



COMMON /SUBR/ 

H0N09 

8A3 

IC 


1 NX»NR#NN#NRH#Eb#ITER#OT»0£TR 

H0N09 

8AA 



2 >A(^><i)>GH(A>6)>H(A>^)>».<A>6)»B<^*2) 

H0N09 

8A5 



3 #AKH(^»<»)#Gf1KW<A^6)#KF(4#4)#P(^#A)#Pl( A#^)»R<<i#4) 

H0N09 

8A6 



A •RKAt A)#HP<A>A)>0<A# A)» WW(A#A)>GMW(A>A)/SS< <»#A)#SC<A) 

H0N09 

8A7 



CQMMON/OLVICt/NREAOfHFILE 

H0NO9 

6A8 

15 


DIMENSION 00<22)#F0(22)#6IC< A#A) 

H0N09 

8A9 


c 


H0N09 

850 



LOGICAL LX 

H0N09 

851 



REAL KF 

H0N09 

852 



DATA L2/6HA /#L1/6HA - KH/ 

H0N09 

853 

2C 


DATA NPD/22/#NPF/22/ 

H0N09 

85A 



DATA 0VBAR/*05/»DC2/.01/»0C3/«01/#0C4/*01/ 

H0NO9 

855 


c 


M0N09 

856 



00 60 I«1#NX 

H0N09 

657 



DO 60 J-1»NR 

H0N09 

858 

2 5 


oO GIC(I»J)«0. 

H0N09 

859 



CALL INPT(3IC»A, A#NREAD) 

H0N09 

860 



CALL HPINXf NR>GIC>NX#AHKFIC»HFILE) 

H0N09 

661 


c 

LOOP COMPUTES NOMINAL FIRST»THEN SENSITIVITIES 

H0N09 

662 



00 999 LQ0P*1#5 

H0N09 

863 

30 

c 


H0N09 

86A 


c 

DEFINE PARAMETERS 

ri0N09 

865 



V8AR«ZP1 

H0N09 

866 



C2-ZP2 

H0N09 

867 



C3-ZP3 

H0N09 

868 

35 


CAAZPA 

H0N09 

669 



SO TO (100#200# 300^ A00#50U)f LOOP 

H0N09 

670 



500 OELT-OCA 

H0N09 

871 



IP«A 

HON09 

8 72 



CA«CA ♦ OCA 

H0N09 

873 

<»0 


GO TO 100 

H0N09 

67A 



AOO 0ELT-DC3 

H0N09 

875 



IP*3 

H0N09 

876 



C3«C3 ♦ 0C3 

H0N09 

877 



GO TJ 100 

H0N09 

878 

^5 


300 DELT-0C2 

H0N09 

879 



IP-2 

H0N09 

880 



C2-C2 ♦ 0C2 

H0N09 

681 



GO TO 100 

rtON09 

862 



200 oelt-dvbar 

H0NC9 

683 

50 


TP-1 

H0N09 

88A 



^8AR*V8AR + OVBAR 

H0N09 

885 



100 CONTINUE 

H0NO9 

886 


c 


H0N09 

887 



VBS-VBAR*V8AR 

H0N09 

688 

55 


ACl>U--.018-*03A*VB$ 

H0N09 

889 



A<1#2)-,0A8*VBS 

iiON09 

890 



A(1,3)-2.A 

H0N09 

691 



A<1#A)— '32t2 t lA. + VBS 

H0N09 

692 



A (2f 1 )-• 13* (\/BAR-#75) /0«75 

H0N09 

893 

60 


IF(V8AR*LT.0.0) A ( 2 » 1 ) -0 •! 6* < *75 AVB AR ) /O* 75 

H0N09 

89A 



A(2#2)— .5 

HdN09 

895 



MZp 3)-0, 

H0N09 

696 



A(2f A)— 180*»VBAR 

H0N09 

897 



A(3#U«0. 

H0N09 

696 

65 


A<3» 2)-,02*VBAR 

H0N09 

899 



A<3,3)— 1*5 + CA 

H0N09 

900 



A(3#A)-A*AV8AR 

H0N09 

901 



IF< VBAR,LT,0*0) A ( 3 , A ) — A* ♦ VB AR 

H0N09 

902 



A(A#l)-0« 

H0N09 

903 

7C 


A(A>2)»0* 

H0NO9 

90A 



A(A,3)-1. 

H0N09 

905 



A<A# A)-0* 

H0N09 

906 



3(1«1)-0«12 

HQN09 

907 



B(l»2)-0. 

MON09 

908 

75 


3(2#l )-*7*VBAP 

H0N09 

909 



3<2#2)— 7.8 - 3.*VBAR 

H0N09 

910 



8<3» I)-. 35 ♦ .12*V8AR ♦ C3 

HON09 

911 


Figure 21 . -Program listing of subroutine MODEL. 
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B(3»2 |«,2^*V3AR ♦ C2 

H0N09 

912 


y (A,U«0# 

H0N09 

913 

yc 

•3(<»»2)f0. 

H0NO9 

919 


r 

H0N09 

915 


00 550 

H0N09 

916 


00 550 J*1>NX 

H0N09 

917 


550 

N0N09 

9ia 

y s 

y(i#2 >fi. 

H0N09 

919 


H ( 2 # 3 ) • U 

H0N09 

920 


^(3, <.)•!• 

H0N09 

921 


C 

H0N09 

922 


OJ 50 I«l»NX 

H0N09 

923 

9 0 

00 5C 

H0N09 

929 



H0N09 

925 


5C GM(r,J)«0. 

H0N09 

926 


. c 

H0N09 

927 


SMC1,U«SURT<2.MSIGXG 

H0N09 

92$ 

9 5 

3H<2,2)*S0RT(2*)»SIGZG 

H0N09 

929 


6M(3t6>-S0RT<2« 02*S I GXG+* 6l*S IGZG ) / A ( 3# 3 ) 

H0N09 

930 


C 

H0N09 

931 


CALL OISC(NTER) 

H0N09 

932 


C 

H0N09 

933 

IOC 

W< Ir 3)«SI0VZ 

H0N09 

939 


/<(2#<»)-SIG0 

H0N09 

935 


W(3,5)*SIGTH 

H0N09 

936 


c 

H0N09 

937 


C DEFINE. WORKING DATA ARRAYS 

H0N09 

936 

10 5 

ODU) -Ad^l ) 

H0N09 

939 


0D(2)«A(1#2) 

H0N09 

990 


DDm«A(l,3) 

H0N09 

991 


00(^l«Aa»9) 

H0N09 

992 


00(5»-A(2# 11 

H0N09 

993 

lie 

DD(o|«A(2»2) 

H0N09 

999 


00(7)«A(2# 3) 

H0N09 

995 



H0N09 

996 


D0(9l-A(3>n 

H0N09 

997 


OD(iO)«A<3»2) 

H0N09 

996 

115 

0D( IU»A( 3# 3) 

H0N09 

999 


00( 12)-A(3, <t) 

H0N09 

950 


00(13)-A(^>1) 

H0N09 

951 


00<14)«A(4,2) 

H0N09 

952 


00(15)«A(<»> 3) 

H0N09 

953 

12C 

00(16)«A(<»»<») 

H0N09 

959 


DD(17)>B(1> U 

H0N09 

955 


0D<ld)«B(l*2) 

H0N09 

956 


00( 19)*B( 2» 1 ) 

HQN09 

957 


00(20)>e(2>2) 

H0N09 

958 

\>b 

00(21 )*B( 3f 1 ) 

rt0N09 

959 


l)P<22)-B(3r2) 

H0N09 

960 


C 

H0N09 

961 


C •>»3VIDF OlAK STARTING GAINS 

H0N09 

962 


DU 150 I«UNX 

H0N09 

963 

130 

JU 150 J»1>NR 

H0N09 

969 


150 KF ( I, J ) .f,IC ( 1 » J J 

H0N09 

965 


r 

H0N09 

966 


CALL OIAK 

H0N09 

967 


C 

H0N09 

966 

135 

FD<1)-kF(UX) 

H0N09 

969 


F0{2)-KF(2»l) 

H0N09 

970 


FD(3)"KF(3, 1) 

HQN09 

971 


FO(^) -KF(4» 1 ) 

H0N09 

972 


F0(5)«KF(1#2) 

H0N09 

973 

1<»C 

FD(b)«KF(2»2) 

H0N09 

979 


F0(7)«KF(3>2I 

H0N09 

975 


FD(B)«<F(4,2 » 

H0N09 

976 


FC(9)»KF( 1, 3) 

H0N09 

977 


F0(10)»KF(2#3) 

H0N09 

978 

1A5 

FOdl )«KF (3» 3) 

^0N09 

979 


F0d2 »*KF(<», 3) 

H0N09 

980 


IF(OETR.GT.O.O) F 0 ( 1 3 ) - ALOG ( Of TR ) 

H0N09 

981 


FDd<t»«RId#U 

HON09 

962 


FDd5)«RId»2) 

H0N09 

963 

15C 

FD( 16)-PId/3) 

H0N09 

989 


FD(l7l»RI(2d) 

H0N09 

965 


FDda )»RI(2f 2) 

H0N09 

986 


F0( 19J«RI (2»3) 

H0N09 

987 


F0(20)-Pl ( 3#n 

H0N09 

966 

155 

F0(21)»RI(3»2I 

H0H09 

969 


*=D(22)«R1 ( 3» 3) 

H0N09 

990 


r 

H0N09 

991 


Figure 21. -Continued. 
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IHLOOP-Nc-l) G3 TO 900 

H0N09 

992 


C 


HCN09 

993 

16C 


30 115 I-1#NP0 

N0N09 

999 


115 

oncH-n»oo(u 

H0N09 

995 



00 120 I-l-NPF 

H0N09 

996 


120 

MICH#n«FO(I) 

H0N09 

997 


C 


H0N09 

998 

165 


IF(LXCS)) 30 TO 999 

H0N09 

999 



WRIT6(HFILf:,n99) ICH 

•40N09 

1000 



»»PITE<MFUc#l2Cn SIGV7fSlG0#SIGTH^SIGXGf SI6Z& 

H0N09 

1001 



DO 1209 I-WNX 

HQN09 

1002 


1209 

o(I#U • SOPKPII,!)) 

H0N09 

1003 

1 /o 


(HFUE» 1202) (P< I* n> I«1*NX) 

H0N09 

1009 



WRITE<HMLE#1205) 

H0N09 

1009 



DO 1208 I-UNX 

H0N09 

1006 


1208 

WRITF(HFILc>1203) < A ( I# J )/ J-1#NX) 

H0N09 

1007 


C 

CALL PQLES<^>A>RrSS#h^DT#t2) 

H0N09 

1000 

175 


WRITE(NFILc#1205) 

H0N09 

1009 



DO 1207 I-l»NX 

H0N09 

1010 


1207 

WRITE(HFILE-120A) C AKH( 1 - J W J •! * NX ) 

N0N09 

1011 


C 

CALL POLfcS(A#AKH#A#SS#H#DT#Ll) 

H0N09 

1012 


1199 

FORMAT! IHl# 19HFUTFR DATA ICH • #I2> 

H0N09 

1013 

130 

1201 

fORMATdX dlHSTATS IN -9E12-5) 

H0N09 

1019 


1202 

FORHATdX dlHSTATS OUT •9E12.5) 

N0N09 

1015 


1203 

FORMATdXdlHA #<i612-5) 

H0H09 

1016 


I20<t 

FORMATdX, llrtAHH #<iE12.5) 

H0N09 

1017 


1205 

FORMAT! IHO) 

H0N09 

1018 

185 

C 


H0N09 

1019 



30 TO 999 

HQN09 

1020 


C 


H0N09 

1021 


900 

CONTINUE 

H0N09 

1022 


C 

SENSITIVITIES 

H0N09 

1023 

19C 


00 910 I«1#NP0 

H0N09 

1029 


910 

OSdCH,IP#I)» !00!n - D<ICH»I))/0£LT 

H0N09 

1025 


C 

GAIN SENSITIVITIES 

H0N09 

1026 



JO 920 I»l#12 

H0N09 

1027 


920 

GKdCH-Iodl-IFOm - FdCH,in/OELT 

NQN09 

1026 

195 

999 

CONTINUE 

H0N09 

1029 


C 


H0N09 

1030 



RETURN 

N0N09 

1031 



END 

H0N09 

1032 


SYMBOLIC R«:F6R£NCt HAP (R-i) 

ENIRY POINTS 
3 MCOFL 


rARlA:^! 

. :S 

SN TYPE 

R 

ELKATION 






It 

A 

REAL 

ARRAY 

SUdR 

190 

AKH 

REAL 

ARRAY 

SUBR 

i3C 

3 

REAL 

ARRA7 

SUflP 

3 

CIHP 

REAL 


SU8RH 

t>77 

C2 

REAL 



9 

C2HP 

REAL 


SUBRH 

700 

C3 

REAL 



701 

C9 

REAL 



315 

0 

REAL 

ARRAY 

DAT 

590 

0C2 

REAL 



591 

DC 3 

REAl 



592 

0C9 

REAL 



/05 

DO 

REAL 

ARRAY 


702 

OCLT 

REAL 



7 

DETF 

REAL 


SU0R 

1 

OHP 

REAL 


SUBRH 

119 

OS 

REAL 

ARRAY 

SENSP 

6 

DT 

REAL 


SUBR 

0 

OTHP 

REAL 


SU6RH 

537 

DVBAR 

REAL 



9 

EE 

REAL 


SUBR 

973 

EIP 

REAL 


OAT 

979 

c2P 

REAL 


DAT 

137 

F 

REAL 

ARRAY 

DAT 

733 

FO 

REAL 

ARRAY 


1009 

6E 

REAL 

ARRAY 

SENSP 

761 

GIC 

REAL 

ARRAY 


1090 

GK 

REAL 

ARRAY 

SENSP 

loic 

GL 

REAL 

ARRAY 

SENSP 

30 

6H 

REAL 

ARRAY 

SUBR 

160 

GMKW 

REAL 

ARRAY 

SUBR 

910 

6HW 

REAL 

ARRAY 

SUBR 

1019 

GSOE 

Rt AL 

ARRAY 

SENSP 

1026 

GSQL 

REAL 

ARRAY 

SENSP 

60 

H 

»E AL 

ARRAY 

SUBR 

330 

HP 

REAL 

ARRAY 

SUBR 

673 

I 

INTEGER 



0 

ICH 

INTEGER 


F.P* 

703 

IP 

INTEGER 



5 

ITER 

INTEGER 


SUBR 

679 

J 

INTEGER 



210 

KF 

REAL 

ARRAY 

SUBR 

67 5 

LOOP 

INTEGER 



62 

LX 

LOGICAL 

ARRAY 

VARDAT 

5 39 

LI 

INTEGER 



533 

L2 

INTEGER 



1 

MFUE 

INTEGER 


DEVICE 

2 

NN 

INTEGER 


SUBR 

535 

NPO 

INTEGER 



536 

NPF 

INTEGER 



1 

NR 

INTEGER 


SUBR 

0 

NREAD 

INTEGER 


DEVICE 

3 

NRK 

INTEGER 


SUBR 

0 

NTER 

INTEGER 


F.P# 
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0 


INTtGM 


iUBR 

230 

P 

REAL 

ARRAY 

SUBR 



REA'l 

ARRAY 

SU9R 

J&O 

0 

REAL 

ARRAY 

SU8R 

^ TO 



ARRAY 

SUBR 

310 

RI 

REAL 

ARRAY 

SU8R 


SC 

HfAL 

ARRAY 

SUBR 

0 

SIGO 

REAL 


F.P. 

<.7*. 

SIGSO 

REAL 


DAT 

0 

SIGTM 

REAL 


F.P. 

0 

SIGVZ 

REAL 


F.P. 

0 

SIGX6 

REAL 


F.P. 

C 

SIGZG 

REAL 


F.P. 

A30 

S$ 

real 

ARRAY 

SUBR 

0 

ux 

REAL 

ARRAY 

VAROAT 

676 

VBAR 

REAL 



70^. 

VBS 

REAL 



100 

W 

REAL 

ARRAY 

SUBR 

2 

(^hp 

REAL 


SUBRH 

370 

WW 

REAL 

ARRAY 

SUBR 

0 

X 

REAL 

ARRAY 

TAT 

0 

xs 

REAL 

ARRAY 

SENSP 

0 

ZPl 

REAL 


F*P. 

0 

ZP2 

REAL 


F.P. 

0 

ZP3 

REAL 


F*P. 

0 

ZPA 

REAL 


F.P. 


RX T cKNALS 


TYPE 

AR3S 












A log 


PEAL 

1 

LI3RARY 



OIAK 




0 




DISC 



1 




INPT 



4 




iP 



6 




SORT 

REAL 


1 

LIBRARY 

STATtMeMT LABhLS 












C 

50 





0 

60 






77 

100 

0 

115 





0 

120 






0 

150 

73 

2 CO 





66 

30C 






61 

400 


500 





0 

550 






457 

900 

C 

^10 





0 

920 






512 

999 

60 7 

1 199 


fmt 



613 

1201 

FHl 





6X7 

1202 

623 

1203 


F»1T 



627 

1204 

FMT 





633 

1205 

0 

1207 





0 

1206 






0 

1209 

LOOPS 

LABEL 


Index 

FROM 

-TO 

IFNGTH 

PROPERTIES 







16 

oO 


I 

23 

2y 

13B 


NOT 

INNER 






23 

60 


j 

24 

25 

2B 

INSTACK 








3? 

999 

* 

LJap 

29 

195 

4608 


EXT 

REFS 

NOT 

INNER 




160 

550 


I 

82 

64 

138 


NOT 

INNER 






15l> 

550 


j 

63 

34 

28 

INSTACK 








166 

5C 

* 

I 

89 

92 

13B 


NOT 

INNER 






172 

50 


J 

90 

92 

3B 

INSTACK 








26& 

150 

* 

1 

l?9 

131 

13B 


NOT 

INNER 






271 

150 


J 

130 

131 

30 

INSTACK 








350 

115 


I 

160 

161 

33 

INSTACK 








360 

120 


I 

162 

163 

38 

INSTACK 








371 

1209 

* 

I 

168 

169 

78 


EXT 

REFS 






AD3 


* 

1 

170 

173 

108 


EXT 

REFS 






417 

1208 

* 

I 

172 

173 

178 


EXT 

REFS 

NOT 

INNER 




422 


* 

J 

173 

173 

LOB 


6XT 

REFS 






440 

120 7 

* 

I 

176 

177 

178 


EXT 

REFS 

NOT 

INNER 




443 


* 

J 

177 

177 

108 


EXT 

REFS 






467 

910 


I 

190 

191 

5B 

INSTACK 








504 

920 


i 

193 

194 

58 

INSTACK 








ruhncN 

BLOCKS 


length 













VARDAT 


ICO 













OAT 


316 













SENS? 


784 













SUBPH 


5 













SUBR 


300 













DEVICE 


? 












static 

ncs 














PROG 

PAH length 


10J6R 

542 









LABELED COHMCN LE^s»H 


274S8 

1509 










Figure 21, -Concluded. 
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iU«RQOri^*lE JISC 

73/7^ DPT-1 FTN 

76/02/14. 

OB .28 

1 


SURPTUTIMfi i)ISC(NTFR) 

H0N09 

503 



CCMMJN /SbdR/ 

H0N09 

504 



1 M(»NJ<#NN>NRM#£e# ITtR»OT#OFTR 

HON09 

505 




H0N09 

506 

1 


3 #AKH<4,4),GMKW(4f6) »KF<4#4)#P(4f 4>,P1(4,4»,R(4#4) 

H0NO9 

507 



4 >Rn4»4)«HP(4#4)#Q<4#4)»WM(4*4) *GHW(4*4)>SS(4»4)»SC(4) 

H0N09 

508 



R 1- A L K F 1 

H0N09 

509 



DATA NU/2/ 

H0N09 

510 


C 

?# PI# k# kX tSfO AS LURKING ARRAYS 

H0N09 

511 

IC 

C 

A#B AND on ARE replaced WITH DISCRETE VERSIUN 

H0N09 

512 




H0N09 

513 


r 

P1»I ♦ (A*on ♦ ( (A^DT)^*2)/2 ♦ • . 

H0N09 

514 


f 


H0N09 

515 


r 

P-DT ♦ (U*0T)**2)/2 ♦ . • 

H0N09 

516 

15 

c 

*i«P 4B 

H0N09 

517 


c 

:>h»P*G«/S0RT(01 ) 

H0N09 

518 


c 


H0N09 

519 



DU 1 I»1#NX 

H0N09 

520 



DO 1 J«1#HX 

H0N09 

521 

20 


P<I# J)#0* 

H0N09 

522 



R(I#J)*0. 

H0N09 

523 



?1(I» J)«0. 

H0N09 

524 



Rn#n»i. 

H0N09 

525 


X 

Pi(i#n»i. 

H0N09 

526 

25 


on 100 n«i#NTtP 

H0N09 

527 



on 10 T»1#NX 

HQN09 

528 



DO 10 J*1#NX 

H0N09 

529 


10 

RKI# J)«Pn#J»*DT/FCOAT(Il) 

H0N09 

530 



DO 11 I*1#NX 

H0N09 

531 

3C 


00 11 J»1#NX 

H0N09 

532 


11 

p<i# j)«p(i#j) 4 Rin#j) 

HQN09 

533 



00 12 I-1#NX 

HON09 

534 



DO 12 J*l#NX 

H0N09 

535 




H0N09 

536 

35 


i)0 12 K*1,SX 

H0N09 

537 


12 

R(I#J)«R(I*4) ♦ A(I#K)«RI(K#J) 

H0N09 

530 



DJ 13 I«1#MX 

H0NO9 

539 



DU 13 J*1#NX 

H0N09 

540 


13 

PKl# J)-P1(I#JI > R(1#J) 

H0N09 

541 

<»0 

IOC 

continue 

H0N09 

542 


c 


H0N09 

543 



DO 15 I»l#NX 

H0N09 

544 



DO 15 J«1#NX 

H0N09 

545 


15 

Ad# j )«pi(r# j) 

H0N09 

546 

<*5 


31»1./SQRT(0T) 

H0N09 

547 



orj 16 I-1#NX 

H0N09 

548 



DO 16 

H0N09 

549 



GHKW( I# J) «D« 

H0N09 

550 



Oil 16 K«1#NX 

H0N09 

551 

50 

16 

GhKW(I» J)«GMKW(I#J) ♦ P( I#K)AGM<K» J) 

H0N09 

552 



DO 10 I»1#NX 

H0N09 

553 



DU lb J*l»NN 

HQN09 

554 


18 

GH(I, J)«GMKW(I# J)*S1 

H0N09 

555 



DO 17 I»1#NX 

H0N09 

556 

55 


DO 17 J«1,NU 

H0N09 

557 



^1(T# JI-0* 

H0N09 

558 



DC 17 K-1#NX 

H0N09 

559 


17 

Pld# J)«PU I ♦ J) ♦ P(I#KM8(K# J) 

H0NC9 

560 



00 19 I«l#NX 

H0N09 

561 

60 


DU 19 J»l#NiJ 

H0N09 

562 


l<i 

3(I#J)-Pin#J) 

H0N09 

563 



R=TURN 

H0N09 

564 



f NT 

H0N09 

565 



SYMBOLIC 

RbFERENCt 

MAP (R -I ) 






ENTRY 

POINTS 








3 

DISC 








VARIABLES $N 

TYPE 

RELOCATION 






10 

A 

REAL 

ARRAY SUBR 

140 

AKri 

REAL 

ARRAY 

SUBR 

130 

6 

REAL 

ARRAY SUSP 

7 

DETR 

REAL 


SUBR 

6 

OT 

R2AL 

SUBR 

4 

Efc 

REAL 


SUBR 

30 

GP 

REAL 

ARRAY SUBR 

160 

GNKW 

REAL 

ARRAY 

SUBR 

410 

GMW 

REAL 

ARRAY SUBR 

60 

H 

REAL 

ARRAY 

SUBR 


Figure 22 . -Program listing of subroutine DISC. 


65 



Appendix F 


3i0 

HP 

pPAt 

ARRAY 

5UBR 

Ztt 

II 

INTEGER 



265 

J 

integer 



210 

KE 

REAL 

ARRAY 

SOBR 

1 

HR 

INTEGER 


SU6R 

C 

HT'-.R 

INTr r,PP 


P.P 

0 

NX 

iNTEGcR 


SL3R 

2 5C 

PI 

REAL 

ARRAY 

SOBR 

2 7c 

9 

REAL 

ARRAY 

SUBR 

^5C 

SC 

REAL 

ARRAY 

SU«R 

2 70 

SI 

REAL 



370 

W w 

REAL 

ARRAY 

SUBP 

EXTERNALS 

TYPE 

APG$ 



SOPT 

REAL 

1 LIBRARY 

inline 

FONCTIONS 

TYPE 

APGS 



PLOAT 

REAL 

X INTRIK 


26<» 

I 

INTEGER 



5 

ITER 

INTEGER 


SUBR 

267 

K 

INTEGER 



2 

NN 

INTEGER 


5UBR 

3 

NPM 

INTEGER 


SUBR 

262 

NU 

INTEGER 



230 

P 

RIAL 

ARRAY 

SUBR 

350 

0 

REAL 

ARRAY 

SUBR 

310 

R1 

REAL 

ARRAY 

SUBR 

^30 

SS 

REAL 

ARRAY 

SUBR 

100 

w 

REAL 

ARRAY 

SUBR 


STATEMCMT LABELS 


0 

1 

0 

10 

0 

12 

0 

13 

C 

16 

0 

17 

0 

19 

0 

ICO 


0 11 
0 15 
0 16 


LOOPS 

LABEL 


INDEX 

fron-Tq 

It 

1 

* 

I 

IB 

29 

23 

1 


J 

19 

29 

35 

100 

♦ 

II 

25 

AO 

3t 

10 


I 

26 

2H 

^5 

10 


J 

27 

2B 

55 

11 

A 

I 

29 

31 

62 

11 


J 

30 

31 

72 

1? 

A 

I 

32 

36 

73 

12 

A 

J 

33 

36 

105 

12 


K 

35 

36 

12C 

13 

A 

I 

37 

39 

125 

13 


J 

3P 

39 

137 

15 

A 

I 


A<» 

1^3 

15 


J 

A3 

AA 

156 

16 

A 

I 

96 

5-J 

157 

16 

A 

J 

a7 

50 

171 

16 


K 

99 

50 

20A 

18 

A 

I 

51 

53 

210 

le 


J 

52 

53 

22C 

17 

A 

I 

5A 

5B 

221 

IT 

A 

J 

55 

56 

233 

17 


K 

57 

56 

2^6 

19 

A 

I 

59 

61 

252 

19 


J 

60 

61 

corrcN 

iiCCK^ 


La noth 




SUBP 


300 




Length 

PROPERTIES 


166 

5B 

INSTACK 

NOT 

INNER 

lOlB 


NOT 

INNER 

16B 

9B 

INSTACK 

NOT 

INNER 

lAB 

38 

INSTACK 

NOT 

INNER 

253 


NOT 

INNER 

210 

AB 

INST3CK 

NOT 

INNER 

lAB 

38 

INSTACK 

NOT 

INNER 

138 

33 

instack 

NOT 

INNER 

258 


NOT 

INNER 

218 

AB 

INSTACK 

NOT 

INNER 

138 

A8 

INSTACK 

NOT 

INNER 

25B 


NOT 

INNER 

219 

AB 

INSTACK 

NOT 

INNER 

138 

33 

INSTACK 

NOT 

INNER 


'*>TATI JTICi 

PRr;GkAf< Length 2 73B 1«7 

Ci lABEL'^O CLf^PTN ILNGTM <»5<i3 300 


Figure 22 , -Concluded. 
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SURPIHT INc 

OlAK 73/7<* ( PT*i FTN 

78/02/l<»* 

08*26 

1 

Sl3R0uriNt OUK 

H0N09 

363 


CuMf*GN / 5U3R / 

H0N09 

384 


1 NX,NR>NNfNRM#Ec9 ITERaOT#OETR 

H0N09 

385 



H0N09 

306 

5 


H0N09 

387 


^ >Ri ( A,4)#SS<<»#4)#SC(4) 

H0N09 

388 


REAL KP 

H0N09 

389 

C 


H0N09 

390 


DO iOl I«l#NX 

H0N09 

391 

1C 

00 101 J«1»NX 

H0N09 

392 


101 Ra>J)»0. 

H0N09 

393 

c 

6MW^> WW^ 

HON09 

394 


JO 1 I*1>NR 

H0N09 

395 


00 2 J*1,NX 

H0M09 

396 

15 

S«0, 

H0N09 

397 


OJ 3 K-l#MN 

H0N09 

398 


3 S ■ S ♦ GM< J,K) Awd *K) 

H0N09 

399 


? GMW(J»T) * S 

H0NO9 

400 


00 1 J*I»NR 

H0N09 

401 

2C 

S-0* 

H0S09 

402 


00 K-l#NN 

H0N09 

403 


^ S • S ♦ W(I »K)*W( J# K1 

H0N09 

404 


WW(I»J) • S 

H0N09 

405 


1 rfW(J>U « s 

H0N09 

406 

2 5 

IT « 0 

H0N09 

407 



H0H09 

406 


800 CONTINUE 

H0N09 

409 


tl • 0* 

H0N09 

410 


IT • IT ♦ 1 

M0N09 

411 

30 C 


H0N09 

412 

C 

AKH • A - Krt# GNKk • G^1 - KU# PI • GMKW GHKW* 

H0N09 

413 


00 10 I»1#NX 

H0N09 

414 


00 12 J«l#NX 

H0N09 

415 


S • A(I»J) 

H0N09 

416 

35 

00 l<t Kml,HR 

H0N09 

417 


1^ S • S - KFn>K)*H(K» J» 

H0N09 

410 


12 AKH(I»J) • S 

H0N09 

419 


03 16 J-1#MN 

H0N09 

420 


S • GHU,J) 

H0N09 

421 


DC 13 K«1»NR 

HON09 

422 


18 S • S - KF( I#K)AW(K# J) 

H0N09 

423 


16 jMKWdfJ) • S 

HQN09 

424 


on 10 J-Id 

H0N09 

425 


j « 0 • 

H0NO9 

426 


00 20 K«1,NN 

H0N09 

427 


20 S • ^ 4- GMK'<^ d» K ) AGHKmI J>K ) 

H0N09 

428 


R1 ( d J) • S 

fl0N09 

429 


10 RKJd) • S 

H0N09 

430 

C 


H0N09 

431 

5C C 

P • AKH P AKH» ♦ GMKW GHKW* 

H0N09 

432 


Call cal 

H0N09 

433 

c 


H0NU9 

434 

c 

KF • (APrt^ ♦ GNW^MHPH^ 4 WW*J-1 

HQN09 

435 


00 36 I-1>NR 

H0N09 

436 

5 5 

DO 32 J-UNX 

HCN09 

437 


S»C, 

H0N09 

438 


DO 3^ K*1#NX 

H0N09 

439 


3^ S • S 4 Hd ,K )4Pl (K# J) 

H0N09 

440 


32 HPd,J) • S 

H0N09 

441 

60 

00 3b J«ld 

H0N09 

442 


S • WW(UJ) 

H0N09 

443 


00 38 K*1.NX 

H0N09 

444 


36 $ • S 4 HP( 1 #K) ) 

H0N09 

445 


Rld»j) • s 

HON09 

446 

65 

36 «I { j, n • s 

HON09 

447 

C 


H0N09 

440 


IFINR .GT. 1) GO TO 333 

H0N09 

449 


OCTK • RU1#1) 

HON09 

450 


«Id#l) • IWRIdfl) 

H0N09 

451 

7C 

GO TO 33A 

H0NO9 

452 


333 CONTINUE 

HON09 

453 


Of; 39 I«l#4P 

H0M09 

454 


DP 39 J-WNR 

H0N09 

455 


39 R(dJ)«Kl(dJ) 

H0N09 

456 

75 C 

CALL TOINVR (If J»NR#NP#RI#NPM, SC>D6TP) 

H0N09 

457 

C 

EXPLICIT 3X3 INVERSE 

H0N09 

458 


OfcTP«R< l» U*(P( 2»2MR(3»3)-ft( 3,2)*R<2»3) ) 

H0N09 

459 


Figure 23. -Program listing of subroutine DIAK. 
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I 






H. 




IOC 


lot; 


lie 


lit; 


1 - »(2»l)*(Rn»?M«(3#3l-Rn#21*R(l,3H 

2 ♦ i>(3,l)*(Rn#2)*R(2#3)-R(2#2)*R(l#3)) 
IF(D€TP«tQ*0*0) STOP 31 

R1 (1/ l)-(R(2f 2) *H<3#3)-R(3f 2»*ft(2f 3M/D£TR 
RI(l#2)-(R(3»2l*P(l>3)-R<l»2MR(3#3))/DeTR 
RI(1» 3)«(R<1,2)*R(2»3|-R(?f 2MR(l» J)1/0ETR 
RI (2fl)-(R(3> l)^R(2#3)-R(2»l)»R(3f 3))/0ETR 
RM2#2»»(P( l#l)*Rn#3)-R(l#3)^R(3#lU/OETP 
Hr(2»3)-(R(2»lMR(l>3KR<X#i)*R<2#3U/0ETR 
l»«(R(2»n*R<3#2HRn#l)*R(2>2n/0eTR 
•U (3f 2J»(P(3»U*P(l»2)-k(l» l)*R(3»2))/OETR 
Rl(3» 3)«(R(1» l)*R(2»2)-R(2»l)*Ra#2n/0ETR 
C 

33A CONTINUE 
C 

30 ^6 J-ltNX 
00 30 J»I»NX 
S - P(1»J) 

iF(s .eo, 01 or to 3 i 

S • ABSUPK WJ)-S)/S) 

IMFl ♦LT. S) f-1 - S 
31 P(I>JI • PUltJ) 

30 ^U»I) • 

30 <»2 J-1>NR 
S • GMwn,J) 

DO <»<» K«i,NX 

<»<i S • S. ♦ An#K)*MP(J»K> 

<r? SC<J) • S 
00 A6 J-1»NR 

s«o, 

Oo ^8 K«l,NR 

<»B S - S ♦ SC(KMPI(K^ J) 

^6 KF(lfJ) • S 
C 

IF(IT .Lc# 1) CO TO 800 
IF(£l .LT, EE) GO TO 220 
IF( IT .GT. IT:-P) go to 220 
GO TO »00 
220 CONTINUE 
C 

RETURN 

END 


H0N09 

HON09 

H0N09 

H0N09 

H0N09 

H0N09 

HON09 

H0N09 

H0N09 

H0N09 

H0N09 

H0N09 

H0N09 

H0N09 

H0NO9 

H0N09 

H0N09 

H0N09 

HON09 

H0N09 

HON09 

H0N09 

(40N09 

H0N09 

H0N09 

H0N09 

H0N09 

H0N09 

H0H09 

H0N09 

H0N09 

H0N09 

H0N09 

rt0N09 

H0N09 

H0N09 

H0N09 

H0N09 

H0N09 

H0N09 

H0N09 

H0N09 


SYHeOUC RfcFERFNCE NAP (W-1) 


ENTRY POINTS 
1 niAK 


VARIA9LFS 

SN TYPt 

RELOCATION 

10 

A 

real 

ARRAY 

SUBR 

130 

3 

REAL 

ARRAY 

SU8R 

6 

OT 

PEAL 


SUBR 

453 

bl 

REAL 



160 

GNKW 

REAL 

ARRAY 

St'BR 

6C 

H 

REAL 

ARRAY 

SUBR 

44t 

I 

INTEGER 



c 

ITER 

INTEGER 


SUBR 

451 

K 

INTEGER 



? 

NN 

INTEGER 


5U3R 

3 

NRH 

INTEGER 


SU3R 

230 

P 

REAL 

ARRAY 

SUBR 

350 

•3 

RFAL 

ARRAY 

SUBR 

310 

RI 

REAL 

ARRAY 

SU3R 

45C 

SC 

REAL 

ARRAY 

SUBR 

IOC 

w 

real 

ARRAY 

SUBR 

EXTERNALS 

TYPE 

ARCS 



CAL 


0 


INI INE 

FUNCTIONS TYPE 

ARCS 



ARS 

REAL 

1 INTPIN 


140 

AKH 

REAL 

ARRAY 

SUBR 

7 

OtTR 

REAL 


SUBR 

4 

EE 

REAL 


SUBR 

30 

GH 

REAL 

ARRAY 

SUBR 

410 

GKW 

REAL 

ARRAY 

SUBR 

330 

HP 

REAL 

ARRAY 

SUBR 

452 

IT 

INTEGER 



44 7 

J 

INTEGER 



210 

KF 

REAL 

ARRAY 

SUBR 

1 

NR 

INTEGER 


SUBR 

0 

NX 

INTEGER 


SUBR 

250 

PI 

RFAL 

ARRAY 

SUBR 

270 

R 

REAL 

ARRAY 

SUBR 

450 

S 

REAL 



430 

$S 

REAL 

ARRAY 

SUBR 

370 

ww 

REAL 

ARRAY 

SUBR 


Figure 23, -Continued, 


460 

461 

462 

463 

464 

465 

466 

467 

468 

469 

470 

471 

472 

473 

474 

475 

476 

477 
476 

479 

480 
401 
462 
483 
404 

485 

486 

487 
486 

489 

490 

491 

492 

493 

494 

495 

496 

497 

498 

499 

500 

501 
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STATineNT LABfLS 
0 1 

0 

C Lh 
0 20 
0 32 

0 3P 
0 

0 lol 
3<»3 33<f 


LOOPS 
3 
iC 
17 
2C 
30 3 

I 

5>3 

10 

7t 12 

no i*i 

123 16 

135 le 
150 10 

160 20 

200 36 

201 32 

211 34 

224 36 

236 33 

261 39 

265 39 

344 46 

354 30 

366 42 

40C 44 

412 46 

42C 48 


IND6K 

♦ I 
J 

♦ I 

♦ J 

K 

♦ J 

K 

♦ I 

♦ J 
K 

♦ J 

K 

♦ J 

K 

♦ I 

♦ J 

K 

♦ J 

K 

♦ I 
J 

♦ I 

J 

♦ J 

K 

♦ J 

K 


LABEL 

101 

101 

1 


COMMON blocks LtNGTH 
SUBR 300 


C 2 
0 lo 
0 16 
0 30 

0 34 

0 29 

0 46 

442 220 

72 800 


FROM 

-TO 

LENGTH 

PROPERTIES 


9 

11 

13B 


NOT 

INNiR 

10 

11 

26 

INSTACK 



13 

24 

526 


NOT 

INNER 

14 

13 

22 3 


NOT 

INNFR 

IS 

17 

48 

INSTACK 



19 

24 

233 


NOT 

INNER 

21 

22 

49 

ikstack 



32 

4S 

1008 


NOT 

INNER 

33 

37 

243 


NOT 

INNER 

3S 

3b 

4B 

INSTACK 



36 

42 

246 


NOT 

INNER 

40 

41 

40 

INSTACK 



43 

4b 

238 


NOT 

INNER 

45 

46 

4B 

INSTACK 



54 

65 

538 


Ncr 

INNER 

55 

59 

228 


NOT 

INNER 

57 

5d 

40 

INSTACK 



60 

65 

250 


NOT 

INNER 

62 

63 

46 

INSTACK 



72 

74 

138 


NOT 

INNER 

73 

74 

3B 

INSTACK 



93 

110 

70B 


NOT 

INNER 

94 

100 

108 

OPT 



101 

105 

23B 


NOT 

INNER 

103 

104 

48 

INSTACK 



106 

IID 

20 B 


NOT 

INNER 

1U8 

109 

4B 

INSTACK 




0 

0 

0 

361 

0 

0 

0 

260 


STATISTICS 

PROGRAM LENGTH 456B 302 

CM LABELED COMMON LENGTH <»54p 300 


Figure 23, -Concluded. 


3 

12 

18 

31 

36 

42 

48 

333 
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! f’t :i I iNi. 

CAL 

73/7<» CM«1 

FTN A#6fA39 

7S/02/1A. 

06.28< 

1 



subpdjtinc cal 


H0N09 

181 




COHMtlN / SUSP / 


H0N09 

182 




1 NX,NR/KN#NPM,ec^ iter,ot#d;:Tp 


H0N09 

183 




2 >A{^,5),GM{<»»6)>H(A,^|#W(4»6)^B<<»>2) 


HQN09 

18A 




3 #AKH{^f )#GNKW(^>t) 

H0N09 

185 




<» »»I(<»>AWhP<A><»),w(<i#4)#Wk(^#<»UGMW(^#<i)>SS(4a4)»SC(^) 

H0N09 

166 




PEAL KF 


H0N09 

187 


c 




H0N09 

168 


c 

•« 

- AKF 


H0N09 

189 

10 



DO 10 I-1»NX 


H0N09 

190 




00 10 J«1,NX 


H0N09 

191 



10 

• AKHINJl 


40N09 

192 




IT • 0 


H0N09 

193 


c 




H0N09 

19A 

15 


POO 

CONTINUE 


H0N09 

195 




El • C. 


H0N09 

196 




11 • IT +1 


H0N09 

197 


c 




H0N09 

196 


c 

Pi 

• 3 PI 04 + PI 


H0N09 

199 

ZC 



00 20 I«1#NX 


H0N09 

200 




OC 2G J«1>NX 


H0N09 

201 




S-0. 


H0N09 

202 




OQ 2<i K»1»NX 


H0N09 

203 




S - S ♦ Q(I#K)*PUK#J) 


H0H09 

20A 

25 


20 

SS(I,J) « S 


H0N09 

205 




00 26 I -It NX 


H0N09 

206 




OU 26 J-I»NX 


H0N09 

207 




S-0* 


H0N09 

208 




00 23 K-ltNX 


H0N09 

209 

30 


2P 

$ • S ♦ SS( ltK)*Q(J,K) 


H0N09 

210 




S3 - PKitJ) 


H0N09 

211 




IF(53 .60* 0) GO TO 30 


H0H09 

212 




S3 - ABS($/S3) 


H0N09 

213 




IF(tl .LT* S3) 51 • S3 


H0NO9 

21A 

3 5 


30 

Pm#J» • PUItJ) ♦ S 


H0N09 

215 



26 

PKjf n • PKit j) 


HON 09 

216 


* 




HQN09 

217 


c 

Q 

•QQ 


HON09 

218 




S3 • 0« 


H0N09 

219 

<»C 



O.J 6C i*ltNX 


H0N09 

220 




00 AO J-ltNX 


H0N09 

221 




S-0. 


H0N09 

222 




)J A 2 K-l,N< 


H0N09 

223 



^2 

S • S ♦ 0(1 t'^)*C(Kt J) 


H0N09 

22A 



hO 

SS (It J) • i 


H0N09 

225 




DJ A 5 I-1#nX 


H0N09 

226 




00 AA J-ltNX 


H0N09 

227 




J(ItJ) - SS(ltJ) 


H0N09 

228 



^5 

S3-S?4SS(Itl) 


H0N09 

229 

5 0 

c 




H0N09 

230 




IF (IT .Li. 1 ) GC TO 800 


H0N09 

231 




IF(53 .6T. NX) GO TO 230 


H0N09 

232 




IF(tl .LT. EE) GO TC 220 


H0N09 

233 




IF(IT .GT. ITcP) GO TO 220 


H0N09 

23A 

5 5 



3C Tn 800 


H0N09 

235 


c 




H0N09 

236 



220 

CJNl INUF 


H0N09 

237 




RETURN 


H0N09 

238 



230 

CONTINUE 


H0N09 

239 

60 



;top ai 


H0N09 

2A0 




FNO 


H0N09 

2A1 


SYMBOLIC kEFF.PENCi: MAP («•!) 

lntpy plums 

1 CAL 


VARIABLfcS SN lYPE RELOCATION 


10 

A 

REAL 

ARRAY 

SUBR 

lAO 

AKH 

REAL 

ARRAY 

SUBR 


X 30 

fl 

REAL 

ARRAY 

StBR 

7 

DETR 

REAl 


SUBR 

' 7 = 

6 

Ot 

PEAL 


SUBR 

A 

E5 

REAL 


SUBR 

173 

cl 

Rt AL 



30 

GM 

REAL 

ARRAY 

SUBR 


160 

GMKW 

REAL 

ARRAY 

S08R 

AlO 

GNW 

REAL 

ARRAY 

SUBR 


60 

H 

REAL 

ARRAY 

SUBR 

330 

HP 

REAL 

ARRAY 

SUBR 


170 

I 

INTEGER 



172 

IT 

INTEGER 





Figure 24, -Program listing of subroutine CAL. 
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1 

IT'^R 

INTEGER 


5UBR 

1T5 

K 

INTEGER 



2 

NN 

INTcGEP 


SUBR 

3 

NRM 

INTEGER 


SUBR 

230 

P 

REAL 

ARRAY 

SUBR 

3i>0 

0 

' real 

ARRAY 

SUBR 

ilO 

PI 

real 

ARRAY 

SUBR 

‘♦Ijb 

ii- 

real 

ARRAY 

SUBR 

17t 

S3 

real 



3 70 

VW 

REAL 

ARRAY 

SUBR 


in 

J 

INTEGER 



210 

KF 

REAL 

ARRAY 

SUBR 

1 

NP 

INTEGER 


SUBR 

0 

NX 

INTEGER 


SUBR 

250 

PI 

REAL 

ARRAY 

SUBR 

270 

R 

REAL 

ARRAY 

SUBR 

17A 

$ 

REAL 



A30 

SS 

REAL 

ARRAY 

SUBR 

100 

u 

REAL 

ARRAY 

SUBR 


INLINE’ fUNCTnNS TYPt 
ABS BfAL 


ABGS 

1 INTRIN 


STAfE^tNT LABtlS 


C 

10 



0 

20 



0 

24 

c 

26 



0 

28 



76 

30 

0 




0 

42 



D 

44 

0 




166 

220 



167 

230 

1 7 

BOO 









L OOPS 

LABEL 

INDEX 

FR^OM-TO 

length 

PROPERTIES 




3 

10 

♦ I 

10 12 

138 


NOT 

INNER 



7 

10 

J 

11 12 

3B 

IhSTACK 





22 

2C 

♦ 1 

20 25 

25B 


NOT 

INNER 



23 

20 

4 J 

21 25 

22B 


NOT 

INNER 



33 

2A 

K 

23 24 

4B 

INSTACK 





50 

26 

♦ 1 

26 36 

40B 


NOT 

INNER 



52 

26 

♦ J 

27 36 

33B 


NOT 

INNER 



62 

26 

K 

29 30 

4B 

INSTACK 





112 


♦ 1 

40 45 

25B 


NOT 

INNER 



113 

hO 

♦ J 

41 45 

226 


NOT 

INNER 



123 

42 

K 

43 44 

4B 

INSTACK 





IhO 

45 

♦ I 

46 49 

168 


NOT 

INNER 



lAA 

44 

J 

47 48 

3B 

TNSTACK 






COMHGN BLOCKS 
SUBR 


LENGTH 

300 


STATISTICS 

PROGRAM LENGTH 2018 129 

CM LAdcLEO COMMON LENGTH 300 


Figure 24 . -Concluded, 


SuBRlUTI Me 

MP 

73/74 OPT-1 

FTN 4 ,6+ 4 39 

78/02/14. 

0B.2B 

1 


SUBROUTINE MP ( N N2 # A# NM, L A BE L# NO ) 


H0N09 

1034 



DIMFNSICN A(NM,U 


H0N09 

1035 



IF(Nl,EQ.l) GO TO 20 


H0N09 

1036 



IF(N2.E0,U go TO 30 


H0N09 

1037 

5 


WRITE (NO# 600) LABEL 


H0N09 

1038 


600 

F JPHAT ( IHO# 7HMATRIX #A4) 


H0N09 

1039 



00 103 I«1#N1 


H0N09 

1040 


lOU 

4RITE(N0#601) ( A ( 1# J } # J » 1# N2 ) 


H0N09 

1041 


601 

FCRMAT(lX#i0eil*4) 


H0N09 

1042 

10 


RFTUPN 


H0N09 

1043 


cC 

rfPI TF. (N0#610) LABEL 


H0N09 

1044 


610 

FORMATdHO# 12HR0W VECTOR #A4) 


HQN09 

1045 



WPTTt(N0»601) (A( 1# I )#I-l#N2) 


H0N09 

1046 



RETURN 


H0N09 

1047 

15 

30 

WRITE (N0#620» LABEL 


H0N09 

1048 


620 

FORMATUHO# 12HC0L VECTOR #A4) 


H0N09 

1049 



WRITE (NG# 601 ) (A<1#1 )#!•!# Nl) 


H0N09 

1050 



RETURN 


H0N09 

1051 



= ND 


HON09 

1052 


SY>(B11IC R6(=E«tNCE MRP (R*l) 

Figure 25. -Program listing of subroutine MP. 
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fMKT POINTS 
? MP 


VAKIABLtS SN TYPt RELOCATION 


0 

A 

REAL 

ARRAY 

F.P. 


150 I 


INTEGER 



lt»l 

J 

INTEGER 




0 LABEL 

INTEGER 


F#P 

0 

NH 

INTEGER 


F.P# 


0 NO 


INTEGER 


F.P 

c 

N1 

INTEGER 


F*P* 


0 N2 


INTEGER 


F.P 

statmht labels 









Ab 

2C 



6A 

30 




0 

100 

lo2 

600 

PHI 


113 

601 

FMT 



122 

6X0 

lAC 

620 

FMT 









L?OPS 

LABEL 

INDEX 

FRON-TO 

1 ?NGTH 

PROPERTIES 





2 Sf 

100 

♦ I 

7 3 

200 


EXT 

REFS 

NOT INNER 



30 


♦ J 

8 8 

IIB 


tXT 

REFS 




?2 


♦ I 

13 n 

103 


EXT 

REFS 





bTAmilCS 

PRC'iPA^ I -NGTh 1703 120 


Figvire 25. -Concluded. 


StBP'^LT INE 

FHIC 73/7A OMT-1 

FTN A.6+A39 

78/02/lA. 

08*26*50 

1 

SUBRjUTIM FHIC(L>Y,X) 


H0NO9 

578 


CLMMCN/SUBRH/ CT #n>W#Cl^C2 


H0N09 

579 


0 I ME NS ION <(2> 


H0N09 

560 


1 C2-2*»D*»^ 


H0N09 

56X 

5 

Cl*W*W 


H0N09 

58^ 


X(1)«U/C1 


H0NO9 

563 


X<2)»0. 


HON09 

56A 


Y-0* 


H0N09 

585 


RfTURN 


H0N09 

586 

10 

END 


H0N09 

567 


SYHBOIIC REFERhNCe P*P (R»l) 


ENTRY P.JINT5 
3 FHlC 


VARIABLES 

SN TYPE 

kcLOr aTion 





3 Cl 

REAL 


SUBRH 

. A 

C2 

REAL 

SU8RH 

1 D 

REAL 


SUBRH 

0 

or 

REAL 

$UBRH 

0 'J 

REAL 


p.p. 

2 

w 

REAL 

SU8RH 

0 X 

REAL 

ARRAY 

F.P. 

0 

Y 

REAL 

F«P* 

STATEMENT 

LAPELS 







0 1 

inactive 







CCMMCN BLOCKS LENoTH 







SlJMPH 5 







STATISTICS 








PRCSPAM 

1 FNGTH 

15P 

13 





CM labeled cumm-jn length 


s 






Figure 26. -Program listing of subroutine FHIC. 
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SL.:^CQtTINF IdPT 73/7<» C'PT-l 


FTN 78/02/l<»* 08.26.50 


1 

SUPfiOUTlNc INPT(A#II#JJ#N) 

H0N09 

T35 


i)lHcN$iaNA{II*JJ)»I0(5I,J0(5WY0(5J 

HON09 

736 


2 FORMAT! 5(2I2#E12. 5) ) 

H0N09 

737 


1 RtAC(N, 2M IG(L ) » JDUI#YO(L) »L«1#5) 

H0N09 

738 

5 

IF(I0<1))3>10#3 

H0N09 

739 


3 00 5 L«l>!? 

H0N09 

7^0 


iFdou n<»* Lp^ 

H0N09 

7kl 


A I*IOU ) 

40N09 

742 


J-JO(U 

H0N09 

743 

1C 

A(I,J)«YO(L) 

H0N09 

744 


GO TO 1 

H0N09 

745 


10 CONTINUE 

H0N09 

746 


RETURN 

H0NO9 

747 


cNO 

H0N09 

748 


SY^'60LIC REFERENCE HAP (R«l) 

ENTRY POINTS 
3 INPT 


INTEGER 
INTEGER 

INTEGER ARRAY 

INTEGER 

REAL ARRAY 

STATEMENT LABELS 


17 

1 



53 

2 FMT 

C 

4 

INACTIVE 

0 

5 

lonps 

LABEL 

INDEX 

FROM-TO 

LENGTH 

PROPERTIES 

22 


♦ L 

4 4 

120 

EXT REFS 

37 

5 

♦ L 

6 10 

lOB 

OPT EXITS 


VAklA0LtS 

SN TYPE 

RELOCATION 



0 

A 

real 

ARRAY F.P. 

70 

I 

72 

in 

INTEGER 

ARRAY 

J 

II 

71 

J 

INTEGER 


77 

JO 

0 


INTEGER 

F.P* 

67 

L 

0 

N 

INTEGER 

F.P, 

104 

YO 


0 

50 


F#Pt 


3 

10 


INACTIVE 


STATISTICS 

PRCGRAH LENGTH 


1218 01 


Figure 27. -Program listing of subroutine INPT. 
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APPENDIX G 

SAMPLE RUN OF PCMLE INITIALIZATION SOFTWARE 

This appsndix consists of £i sample run of PCMLE initialization software, 
as presented in Figure 28. 


1* T 

6 420 3 

2= F 

4 311 

3= F 

4S F 

bs F 

.0313 

,5000 

6.0000 

l.^JbOO 

4,0000 

5.0000 

.6000 

2,0000 

• 7000 


4,0000 

,0026 

.0175 

l.OOOO 

1,0000 

.0100 

3,2200 

.1000 

3.2200 . 

• 1000 

5.0000 

b.OOOO 

5.0000 

5,0000 


1.0000 

-l.OOOO 

-1.0000 

-l.OOOO 


.0010 

.1000 

.00 30 

.0X00 


l.ooon 

,5000 

0.0000 

0,0000 

0.0000 

-.0700 

,0532 

.0190 

o.oooo 

0.0000 

-.0500 

.0087 

-.0208 

0,0000 

0.0000 

-.0450 

-.2368 

.2707 

0.0000 

0.0000 


SAMPLE TIME= ,3l2<iCE-01 


6* F 7= F 8* F 9a F F 

0,0000 0.0000 10.0000 

13,8000 .2500 

f 


matrix KFIC 

• UOOt-Ol ,4600E^O? -,U70E^00 
.aooot^oi .8500t + 0? -,U70E401 

0, ,A600E+nn 0. 

0, .6800E-01 .7o00E-0a 


filter UaTA 
stats in 
stats out 


A 

a 

A 

A 


AKH 

AKH 

AKH 

AKH 


ICH a 1 
,<^0 00 0t.f>l 

,39599E+01 


,9983eE+n0 
. I8589t-0P 
.573A9E-06 
, 600b7L-0rt 


,99B3Bt.00 

,185B9L-02 

♦57349L-06 

,600b7t-n8 


,P60e0E-o2 

,87400F*00 


•Ib082E-02 

,98448t^00 

,60577F-03 

•95934£-Ub 


-,22l30t-0l 

,939g8L*00 

.5734BE-03 

,16701t-04 


, 17500e-01 
,5281 lt-02 


•64419E-01 
-,8600TE«01 
•95474L+00 
.3052 rt-oi 


-,B932bE40l 
-, 14170E*02 
.17594E+00 
-. 15542E-01 


, 1 OOOOF.*Oi 

.1143BF-02 


-,56B2Bt>00 
-.55854E+01 
•12039E+00 
, 10019E^01 


-,12767E^01 

-,bl650E*01 

,UR67F>00 

,99769E^00 


. 1 OOOOE^Ol 


CHANNEL MODELS 1 
D array 
,99B38E*00 
,A5082E*02 
.64419E-01 
-,56820E+OO 
,10509E-O2 
,98448E^00 
-,86007E-01 
-,55854E^0I 
,57349E-U6 
,60577E-03 
,95474E+00 


GRAD D(VBAR) 

-,2l654E-0? 

.30709L-02 

,)3622L-01 

.BOepSE^OO 

.74332E-02 

-,277b4E-04 

-,B5836E-01 

-,55897t401 

,29bl9E-05 

,60618E^03 

,l85b2E-02 


GRAD D(C2) 

0 , 

O, 

O, 

0 , 

0, 

0 . 

0. 

0, 

0, 

0, 

0 . 


GRAD D(C3) 

0 , 

0. 

O. 

O, 

O. 

0, 

0 . 

0 , 

0 , 

O. 

0 , 


6PA0 0(C4) 

.10623E-09 

•22323E-06 

.I0444E-02 

.44503E-04 

-.52473E-10 

^.l3938E-06 

-,88999E-^03 

-,27807E-04 

.59353E-08 

.94146E-05 

•29821E-01 


Figure 28. -Sample run of PCMLE initialization software. 
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GRAD K<VHW) 

grad K(C2) 

grad K(C3) 

6RA0 K(C4) 

•20919E-OI 

0. 

0. 

.30195E-02 

-t205<f4E-02 

0, 

0. 

,A0ll9E-02 

♦29707E-U4 

0* 

0. 

-,25065E»05 

-.43420E-Ob 

0. 

0. 

.1320TE-05 

.21^7lE+02 

0. 

0, 

-,83095E*0X 

.12«J8E^02 

0. 

0. 

-.139»2E*02 

.25734E-01 

0. 

0. 

,16329E*00 

.96002E-03 

0, 

0. 

.86257E-03 

♦21U40E^01 

0. 

0, 

.2H81E*00 

-.26237E^00 

0* 

0. 

,7*300E-01 

t41O06E-O4 

0. 

0. 

-»1A309E-04 

-.50711E-03 

0. 

0. 

.18I97E-03 


MATRIX KfIC 
-.2000E.01 

".6770E+0? 

^•4000E^OO 

•3500E-01 

•036OE^n? 

-.U90E401 

0« 

•3000E+00 

0* 

0. 

•8000E-0 J 

•roOOE-02 


FILTER 

STATS 

uata 

IN 

ICH = 2 

.^♦OOOOE^Ol 

•26080E-02 

•I7500t-01 

•lOOOOE^Ol •lOOOOE^Ol 

STATS 

OUT 

•^7797t*01 

•89407E400 

,474352-02 

.118392-02 

A 


•99917E+O0 

•36194E-03 

•592?7E-0l 

-.095O2E^OO 

A 


-.1B597E-02 

•98449E^00 

-•430032-01 

-.27907E*0l 

A 


-.2862U-0b 

.30188E-03 

• 940lU*OO 

.60425E-01 

A 


-.299e7t-0fi 

•47714E-05 

.304272-01 

.lOOlOE^Ol 

AKH 


• 99917E. + 00 

•36462E-0I 

•867232401 

.685052^00 

AKH 


-.18597E-02 

•93773E+00 

-•9299^2*01 

-.26069E*0l 

AKH 


-,2862It-06 

•26394E-03 

•21595£*00 

•59747E-01 

AKH 


-•29987E-08 

.78I17E-05 

-.14755E-01 

•99643E^00 


CHANNEL models 2 

D ARRAY 

GRAD D<VBAR) 

GRAD D(C2) 

grad D(C3) 

grad D(C4) 

.99917E+00 

111442-02 

0. 

0, 

-.518702-10 

.38194E-03 

,150272-02 

0, 

0, 

.107232-06 

,592272-01 

,701992-02 

0, 

0, 

•908362-03 

-,095O2E*OO 

,460352*00 

0. 

0. 

.214542-04 

-,10b97E-O2 

,743952-02 

0, 

0, 

.131222-10 

,90449E*OO 

-,185852-04 

0. 

0. 

-.34037E-O7 

-,430032-01 

-,856642-01 

0. 

0, 

-,444562-03 

-,27907E^01 

-,558872*01 

0, 

0, 

-.69694E-05 

-,206212-06 

,686912-06 

0, 

0. 

-.29588E-O0 

,301802-03 

,60398£-03 

0. 

,0. 

,468702-05 

,940UE^OO 

,186592-02 

0. 

0. 

.296232-01 

,60425E-01 

,119952*00 

0. 

0, 

.937902-03 

-.29907E-U0 

,719712-08 

0. 

0. 

-.23259E-10 

,47714E-05 

,954462-05 

0. 

0. 

.493242-07 

.30427E-01 

,196592-04 

0. 

0. 

.471172-03 

.lOOlOE^Ol 

,190002-02 

0, 

0, 

.90687E-O5 

,4l72lE-02 

,163572-03 

0, 

,100682-02 

,450782-05 

,11^17E-03 

,404142-05 

•1OO60E-O2 

0, 

.10566E-O5 

,10660E-01 

,212842-01 

0, 

-.450552-03 

-.146332-05 

-,28844E400 

-,932942-01 

-.450552-03 

0. 

-.603222-06 

,127412-01 

,366662-02 

0. 

.30427E-01 

• J9730E-03 

,52256E-02 

,720162-02 

.304272-01 

\ 

0, 

.811482-04 

F array 
-,36061E-01 
,46759E-01 
,179422-04 
-,30403E-05 






Figure 28. -Continued. 
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-, 8613 lE* 0 l 

, 92562 E *01 

. 73216 E *00 

. 45182 E -01 

-, 16801 E «01 

-, 18383 E *00 

. 67764 E -03 

. 45285 E -02 

-, 15703 E *02 

, 59 S 36 E- 0 l 

-, 59256 E *00 

,n 0 S 8 E -01 

-. 59256 E +00 

, 34135 E *05 

., 71031 E *02 

.U 056 E -01 

-, 71031 E »02 

, 32506 E *04 

,J 2 Q 39 E «00 

. 60057 E -08 

, 9 *i 634 E- 0 ‘> 

,30b27E-Ol 

* 10019 E ^01 

, 42 u 91 E -02 

-• 737 UE -04 

. 21330 E -01 

-, 33 b 03 E +00 

*60«164E-02 


1 1866 E 400 

0 - 

0 * 

afl 750 E -02 

31227 E -07 

0 , 

0 . 

. • 46627 E -10 

95666 E -05 

0 . 

0 . 

• 98962 E -07 

j 9 b 83 t -04 

0 - 

o« 

• 47320 E -03 

X 8864 E -02 

0 . 

0 • 

•X 9730 E -04 

235 B 2 K -03 

0 • 

.X 06 X 2 E -02 

. 47077 E -05 

3 Xbbbt -03 

.X 0612 E -02 

0 . 

•X 88 S 8 E -05 

2 X 228 E- 0 X 

0 . 

*, 90055 E ~03 

-. 29389 E -05 

93399 E- 0 X 

-♦ 90055 E -03 

0 . 

-.U 790 E -05 

36854 E -02 

0 , 

. 30527 E- 0 X 

•X 988 XE -03 

, 7 X 967 t --02 

. 30527 E- 0 X 

0 . 

. 79376 E -04 


F AHHAY 
. 23638 E -01 
. 44 A 96 E-OI 
, 32284 E -04 
-• 71301 E-O 5 
. 89969 E ^01 
« 14064 E ^02 
, 77880 E ^00 
, 46070 E-U 1 
, 70844 E *00 
*-, 42041 E ^00 
. 71363 E -03 
, 42 i/ 8 E -02 
-., 15 S 37 E *02 
, 59601 E-O 1 
-. 90503 E *00 
. 25283 E -01 
-, 90 b 03 E 400 
• 28840 E *05 
-. 69719 E ^02 
. 25283 E -01 
-»697X9E^02 
. 32516 E ^04 


6RAD K(VBH) 
. 10498E>00 
-.IOJ^ 42 E-Ol 
.30290E-04 
-•U735E-04 
, 66867 E ^01 
, 15404 E +02 
.19O18E-01 
,29000E-02 
,30O52E*00 
*«66896E^OO 
.22958E-03 
-,6(»070E-03 


GraO K(C 2 ) 

0 . 

0 . 

0, 

0. 

0 . 

0 . 

0 . 

0, 

0 . 

0 . 

0 . 

0, 


6 RA 0 K(C 3 ) 

0, 

0. 

0. 

0. 

0. 

0* 

0. 

Ot 

0* 

Of 

0. 

0, 


ORAD K(C4) 
-♦17602E-02 
.14262E-02 
-•26322E-05 
•82656E*06 
.77239E401 
-.84387E401 
• 18726 E 400 
•10183E-02 
-.15239E>00 
•37392E-01 
-t64499E-04 
♦73678E-04 


Figure 28- -Continued 
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MATRIX 

kfic 


0* 

0. 

0* 

0* 

0, 

0. 

0* 

0. 

0* 

0* 

.3000E-01 

«lOOOE^OO 


r 


filter Data 

STATS IN 

ICH = 3 

•40000L401 

•26080E-02 

•17500E-01 •lOOOOE^Ol •lOOOOE^Oi 

STATS OUT 

•3336/E401 

•93778E^00 

• 642IU-02 •12007E-02 

A ■ 

•99944L*00 

0 • 

.S0O37E-OI -•I0060Et0I 

A 

-•55797E-02 

t98450E+00 

-.17806E-03 .28149E-02 

A 

u. 

0. 

•96231E^00 0, 

A 

0. 

0* 

•30657E-01 .lOOOOE^Ol 

AKH 

•99944E400 

.54521E-0I 

-•17043E-01 -.51563E^00 

AKH 

-.55797E-02 

♦93290E^00 

♦98787E-03 -.35408E-OI 

AKH 

0, 

•34183E-09 

•13628E^00 -•37691E-03 

AKH 

0, 

-.69I70E-06 

~ • 1 3294E:-0 1 • 99535E ♦ 0 0 


CHANNEL MODELS 3 

D array 

GRAD D(VBAK) 

GRAD 0(C2) 

GRAD D(C3) 

GRAD D(C4) 

.99944E+00 

-,b3327t-04 

0. 

0. 

Oo 

0. 

,94184E«04 

o« 

0. 

0, 

.b8037E-01 

.37844E-03 

o« 

Oo 

o98lS4E«03 

-.10060E401 

.25862E-01 

o« 

0. 

Oo 

-•55797E-02 

,74398E-02 

Oo 

Oo 

Oo 

.98450E400 

-.1I306E-05 

Oo 

0. 

Oo 

-♦17806E-03 

-,86089t-0t 

0, 

0. 

-.19209E-05 

.28149E-02 

",55efa3E40I 

Oo 

0. 

Oo 

0« 

-•16JUE-05 

Oo 

0. 

Oo 

0* 

,60837t^03 

Oo 

Oo 

Oo 

,9623IE400 

M9029E-02 

Oo 

Oo 

.30077E-01 

0. 

,l2255E*0n 

Oo 

Oo 

0. 

0. 

-,16859E*07 

Oo 

Oo 

Oo 

0, 

,959I7E-05 

Oo 

Oo 

Oo 

,306S7E-0l 

,l9952E-04 

Oo 

Oo 

.47600E-03 

.lOOOOE+Ol 

.19275E-02 

Oo 

0. 

0. 

.40764E-02 

.12059E-03 

Oo 

.99460F-03 

o35272E-05 

tl8897E-04 

•22784E-03 

.99460E-03 

0. 

•2O3B0E-O6 

-.,1U29E-04 

.21417E-01 

Oo 

-o19397E^05 

-.50808E-08 

-.24186E^00 

-,93050E-01 

-ol9397E-05 

0. 

-•2932BE-09 

.10092E-01 

.368S9E-Q2 

Oo 

.30657E-01 

•15670E-03 

.58249E-03 

,72821E-02 

o30657E-0l 

Oo 

o90441E-05 


F array 
-•5462lE-0i 
,51597E-01 
-,34l83C-09 
.69170E-06 
.75080E-01 
-•U6S9E-02 
.82604E*00 
• 439S>1E-01 
-..49033E^00 
.38303E-U1 
•37691E-03 

.46487E-02 

-.15204E402 

.59244E-01 

•&2226E-04 

-•22586E-02 

,5222bE-04 

.20820E^0b 

-,57504E*02 

-.22586E-02 

-,57584E402 

.32502E^04 


Figure 28, -Continued. 
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GRAD K(VBR) 

GraU K(C2) 

6HA0 K(C3) 

6RAD K(C4) 


0. 

0. 

-.17092E-0T 


0* 

0. 

.13S19E»08 

.36441E-04 

0. 

0. 

,10749E»09 

-,372J5E-05 

0, 

0. 

.13869E-09 


0. 

0, 

♦10643E-01 

.n^^5E*U^ 

0^ 

0. 

.26081E-03 

,45135t>02 

0. 

0. 

,19843E»00 

.43564E-U2 

0. 

0, 

.35J76E-02 

-.?4b07£40l 

0. 

0. 

-.11378E-03 

-•24I81E^00 

0. 

0, 

.75688E-05 

,67689E-03 

0, 

0. 

-.39295E-04 

,52702E-04 

0. 

0. 

.15764E-05 


Figure 28. -Concluded. 
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APPENDIX H 

PROGRAM LISTING OF MODEL REFERENCE ALGORITHM 

This appendix consists of the program listings of subroutines RFMOD 
(Fig. 29), INTGRAL (Fig. 30), and RMIC (Fig. 31). 


SUBROUTINE RFNOO 73/74 0RT<1 FTN 4.6«439 76/02/1$. 07.44.44 


1 


SUBROUTINE RFNOO 

H0N09 

1564 



CONNON/NROAT/AC 3f 3)»B( 3#2I#DX (39I*X U3I#U(2 ) 

H0N09 

1565 



1# n»T2#T3#r4#T5#T6#T7*T8#T9#T10*Tll#THRi#THR2#THR3 

MONO 9 

1566 



2# GU«612#G21#022 

H0N09 

1567 

5 


CONNON/VAROAT/ UX(50)#LXC50I 

HQ NO 9 

1568 



CONNON/IFIC/OUHI59)#XVZC2l#XQf2)#XTH(2>#XOB(2|>XOC<2) 

H0N09 

1569 



LOGICAL LX 

H0N09 

1570 


C 


H0NO9 

1571 



V2»uxan 

H0N09 

1572 

10 


Q«UXU2) 

H0N09 

1573 



TH«UX(13) 

H0N09 

1574 



uai^uxu) 

H0N09 

1375 



U(2)«UX(2I 

H0N09 

1576 



0T-UXU9) 

H0N09 

1577 

15 


CALL FHCQ#Q#XQ) 

H0N09 

1578 



CALL FH<VZ#VZ#XVZ) 

HON09 

1579 



CALL FH(TH#TH»XTH) 

H0N09 

1580 



CALL FH<UU)#UU)»XDB) 

H0N09 

1581 



CALL FH(U(2)«U(2)>XDO 

H0N09 

1582 

20 

C 


H0N09 

1583 



EVZ«X(1)-VZ 

H0N09 

1584 



EQ»X(2)*Q 

H0N09 

1585 



ETH-XC3I-TH 

H0H09 

1586 



UX(24)-EVZ/10« 

H0N09 

1587 

25 


UXC25)*EQ95«73 

H0N09 

1588 



UXC26)«ETH*5.73 

H0N09 

1589 



IFIEQ*GT*1.) EQ«i* 

H0N09 

1590 



IF(EQ«LT«-‘1«I EQ«<»1* 

H0N09 

1591 



IFfETH.6T«l«) ETH«1. 

HQN09 

1592 

30 


IF(ETH,LT.-l.) ETH— X. 

H0NO9 

1593 



IF<EVZ«6T«10*I EVZ«X0o 

H0N09 

1594 



IF(EVZ*LT«-10») EVZ»*10. 

H0N09 

1595 



1F(ABSCEVZ).LT«THRU EVZ-0* 

H0N09 

1596 



IF(ABSCEQI.LT*THR2I EQ-0* 

H0N09 

1597 

35 


IF(ABS(ETH)«LT*THR3) eth«o« 

H0N09 

1598 



Tl«UX<34)/a« * U(I)9U<X)4U(2M>U(2n 

H0N09 

1599 



T5-UX05) 

H0N09 

1600 



T6-UX06I 

H0N09 

1601 



T7«UX(37> 

H0N09 

1602 

40 


TXX-UXOB) 

H0N09 

1603 



P1X«T190,25 

H0N09 

1604 



P22«T190*X944 

H0N09 

16Q5 



P23»T140*06333 

H0N09 

1606 



P33»Tl*l*4167 

H0N09 

1607 

45 


FAC»l. + PUPEVZ*EVZ4P224EQPE042«PP239£Q*ETH4P33*ETHPETH 

H0N09 

1608 



FAC1«$QRT(FAC) 

H0N09 

1609 



$1«P22PEQPVZ ♦ P23PETHPVZ 

H0N09 

1610 



S2-P229EQ9Q ♦ P2396TH9Q 

H0N09 

1611 



S3-P229EQ4TH ♦ P234ETH4TM 

H0N09 

1612 

50 


$4«P229EQ4UC1> 4 P239ETH9U(ll 

H0N09 

1613 



$5«P22«EQ9Ui2) ♦ P239ETH*4U(2) 

H0N09 

1614 


C 


H0N09 

1615 



DX(4)-(GU4P1I4EVZ«VZ ♦ GX29SU4T2 

H0N09 

1616 



DXC5)«(G119P119£VZ9Q ♦ 6129S2)9T3 

H0N09 

1617 

55 


0X(6I»(GU9PUPEVZPTH ♦ GX2PS3I9T4 

H0N09 

1618 



0XC7)»(GXX9P119EVZ9Um ♦ G124S4)4T5 

H0NO9 

1619 



DX(8l«(Gil*P119EVZ9UC2) ♦ 6129S5I«T6 

H0N09 

1620 



0X(9I«(G2X9P1XPEVZ4VZ ♦ 622«$l)PT7 

H0N09 

1621 



DX(10I«IG219PX19EVZ9Q 4 G229S2)9T8 

H0N09 

1622 


Figure 29. -Program listing of subroutine RFMOD, 
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60 


3 . 

DXai)«(S21PPll6EVZ*TH ♦ S22O$3m0 

H0N09 

1623 




DXa2)«(G216Pll6EVZ6U(l) ♦ G226SA)6T10 

H0N09 

1626 




0X<13)*CG216P116EVZ6UC2> ♦ G22^S5)6T11 

H0N09 

1625 


C 

? 

3 

H0N09 

1626 




DO 100 I«l#3 

H0N09 

1627 

65 



s-o« 

H0N09 

1626 




DO 101 J»l,3 

H0N09 

1629 



101 

S»S ♦ A(1#4)6X(J) 

H0N09 

1630 




DO 102 4«1>2 

H0N09 

1631 



102 

s»s ♦ ea#4)«u(4) 

H0N09 

1632 

70 


100 

ox(n«s 

H0N09 

1633 




CALL INTGRAL(0X#X#0T»TINE#13#1I 

H0N09 

1636 


C 


OUTPUT GAINS 

H0N09 

1635 




00 200 !»1»10 

H0N09 

1636 



200 

UXU04il«X(3^I) 

H0N09 

1637 

75 



ux( 20 )*um/io. 

H0N09 

1638 




UX<21)«U(2)*5.73 

H0N09 

1639 




UX(22)»X(U/10, 

H0N09 

1660 




UX(23)-X(3)*5.73 

H0N09 

1661 




UX(27I«XI13) 

H0N09 

1662 

80 



UXC26)«0X(7) 

H0N09 

1663 




UXC29)-0X(8) 

H0N09 

1666 




UX(30I*DX(12) 

HQN09 

1665 




UX(31}»DXf 131 

H0N09 

1666 




UXI32)*FAC 

H0N09 

1667 

85 



RETURN 

H0N09 

1660 




END 

H0N09 

1669 


SYMBOLIC REFERENCE HAP (R*l) 


ENTRY POINTS 
I RFNOD 


variables 

SN TYPE 

RELOCATION 

0 

A 

REAL 

ARRAY 

NRDAT 

355 

DT 

REAL 



17 

DX 

REAL 

ARRAY 

NROAT 

360 

ETH 

REAL 



365 

FAC 

REAL 



123 

Gll 

REAL 


NRDAT 

125 

G21 

REAL 


HROAT 

376 

I 

INTEGER 



62 

LX 

LOGICAL 

ARRAY 

VAROAT 

362 

P22 

REAL 



366 

P33 

REAL 



375 

S 

REAL 



370 

S2 

REAL 



372 

S6 

REAL 



356 

TH 

REAL 



121 

THR2 

REAL 


HROAT 

377 

TINE 

REAL 



116 

TIO 

REAL 


NROAT 

106 

T2 

REAL 


NRDAT 

110 

T6 

REAL 


NROAT 

112 

T6 

REAL 


NROAT 

116 

T8 

REAL 


NRDAT 

103 

U 

REAL 

ARRAY 

NROAT 

352 

VZ 

REAL 



101 

X08 

REAL 

ARRAY 

IPIC 

75 

XQ 

REAL 

ARRAY 

IPIC 

73 

XVZ 

REAL 

ARRAY 

IPIC 


11 

B 

REAL 

ARRAY 

NROAT 

0 

DUN 

REAL 

ARRAY 

IPIC 

357 

£0 

REAL 



356 

EVZ 

REAL 



366 

FACl 

REAL 



126 

G12 

REAL 


NROAT 

126 

G22 

REAL 


NRDAT 

376 

4 

INTEGER 



361 

Pll 

REAL 



363 

P2 3 

REAL 



353 

Q 

REAL 



367 

SI 

REAL 



371 

S3 

REAL 



373 

$5 

REAL 



120 

THRl 

REAL 


NRDAT 

122 

THR3 

REAL 


NRDAT 

105 

T1 

REAL 


NRDAT 

117 

Til 

REAL 


NRDAT 

107 

T3 

REAL 


NRDAT 

111 

T5 

REAL 


NROAT 

113 

T7 

REAL 


NRDAT 

115 

T9 

REAL 


NRDAT 

0 

UX 

REAL 

ARRAY 

VARDAT 

66 

X 

REAL 

ARRAY 

NROAT 

103 

XOC 

REAL 

ARRAY 

IPIC 

77 

XTH 

REAL 

ARRAY 

IPIC 


externals 

TYPE 

ARGS 

FH 


3 

SORT 

REAL 

1 LIBRARY 

INLINE FUNCTIONS 

TYPE 

ARGS 

ABS 

REAL 

1 INTRIN 


INT6RAL 


STATEHENT LABELS 

0 100 0 101 0 102 

0 200 


Figure 29, -Continued, 
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LOOPS 

LABEL 

INDEX 

FRON-TO 

LENGTH 

PROPERTIES 

216 

100 ♦ 

I 

66 70 

27B 

I 

226 

101 

J 

66 67 

6B 

INSTACK 

235 

102 

J 

68 69 

6B 

INSTACK 

252 

200 

1 

73 76 

38 

INSTACK 

CONNON 

BLOCKS 

LENGTH 





HRDAT 

87 





VAROAT 

100 





IPIC 

69 





STATISTICS 

PROCRAH L6NGTH ^OOB 256 

Cn LABSLEO CONNON LENGTH AOOB 256 


NOT INNER 


Figure 29. -Concluded. 


SUBROUTINE INT6RAL 


73/76 ORT«l 


FTN 6«66639 


7B/02/15. 07*66.66 


10 


15 


20 


25 


30 


35 


60 


65 


SUBROUTINE INTGRAL (OX# X#OT>TINE# NX# HINTG) 

H0N09 

750 

CONNON/NRINT/XO<50)>NIHTGP.NX2#BO#B1#CO#C1#C2#TINEO 

H0N09 

751 


H0N09 

752 

THIS SUBROUTINE INTEGRATES OX TO GET X« 

H0N09 

H0N09 

753 

756 

OXNENSION OXm# XUI 

H0N09 

755 


H0N09 

756 

INITIALIZATION 

H0N09 

757 


HON09 

758 

IFCNINTG .GT« 0) GO TO 10 

H0N09 

759 

NX2 • 26NX 

H0N09 

760 

BO • 3«060T/2«0 

H0N09 

761 

Bl • -OT/2.0 

HQN09 

762 

CO • 23.060T/12.0 

H0N09 

763 

Cl •-16»0*0T/12#0 

H0N09 

766 

C2 * 5.0*DT/12»0 

H0N09 

765 

HINTGP • 1 

H0N09 

766 

TINEO • TIHE 

H0N09 

767 

DO 5 I • 1# NX 

HON09 

768 

xom - xm 

H0N09 

769 

RETURN 

H0N09 

770 


HON09 

771 

STARTUP PROCEDURE 

H0N09 

772 


H0N09 

773 

IFININTGP - NINTG) 20* 60* 90 

H0N09 

776 

GO TO (30* 60)* NINTGP 

H0NO9 

775 

DO 35 I • 1* NX 

H0N09 

776 

TEHP • 0T*DX(I) 

H0N09 

777 

Xni ■ X(I) * TEHP 

HON09 

778 

0X(NX4>I) • DX(I) 

H0N09 

779 

GO TO 50 

H0N09 

780 

DO 65 I - 1* NX 

H0N09 

781 

TEHP • B060XCI) ♦ B160X(HX6l) 

H0N09 

782 

X(I) • X(I) - TEHP 

H0N09 

783 

DX(NX26I) • 0X(NX6I) 

H0N09 

786 

DX(NX61) - Dxm 

H0N09 

785 

HINTGP • HINTGP ♦ 1 

H0N09 

786 

TIME • TIHE - OT 

H0N09 

787 

RETURN 

H0N09 

788 


H0N09 

789 

GO TO (87* 70* 80)* NINTG 

H0N09 

790 

DO 75 I • 1* NX 

H0N09 

791 

X(I) » XO(I) 

H0N09 

792 

TEHP • OX(I) 

H0N09 

793 

DX(I) • OX(NX^I) 

HON09 

796 


Figure 30. -Program listing of subroutine INTGRAL. 
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7$ 0XCNX4II • TENP 
GO TO 07 

60 00 80 I • 1$ NX 

xm • xocn 

00 TENP • DX(I) 

OXU) • 0X(NX2«I) 

80 DXINX2^X» • TENP 
si HIMTGP • HINT6P ♦ 1 
TINE - TIMED 

00 90 GO TO (100# 200# 3001# NINTG 

C 

C EULER INTEGRATION (FIRST ORDER) 

C 

100 00 100 I • 1# NX 

60 TEMP • OT^DX(I) 

150 X(I) - X(I) ♦ TEMP 

Gd TO 1000 
C 

C HARRY BEAN INTEGRATION (SECOND ORDER) 

60 C 

200 DO 250 I • 1# NX 

TEMP - BO»OX(I) ♦ B1*0X(NX^I) 

X(I) • X(I) ♦ TEMP 
250 0X(NX>I) « OX(I) 

70 GO to 1000 

C 

C GUNTER STEIN INTEGRATION (THIRD ORDER) 

C 

300 DO 350 I • 1# NX 

75 TENP - C0*DX(I) ♦ C1*0X(NX+I) ♦ C2*0X(NX2^I) 

X(l) • X(I) ♦ TEMP 
DX(NX2>I) • DX(NX^I) 

350 DX(NX«1) • 0X(I) 

C 

80 C INCREMENT TINE 

C 

1000 TIME • TIME ♦ DT 
RETURN 
END 


H0N09 

795 

H0N09 

796 

HON09 

797 

H0N09 

796 

H0N09 

799 

H0N09 

800 

H0N09 

801 

H0N09 

802 

HQN09 

603 

H0N09 

804 

H0N09 

805 

H0N09 

806 

H0N09 

807 

H0H09 

806 

HQN09 

809 

HON09 

810 

H0N09 

811 

H0N09 

812 

H0N09 

813 

H0N09 

814 

H0N09 

810 

H0N09 

816 

H0N09 

817 

H0N09 

818 

H0N09 

819 

H0N09 

620 

H0N09 

821 

H0N09 

822 

H0N09 

823 

H0N09 

824 

H0N09 

825 

H0N09 

826 

H0N09 

827 

H0N09 

828 

HQN09 

829 

H0N09 

830 

H0H09 

831 

H0N09 

832 

H0N09 

833 


CARD NR, SEVERITY DETAILS DIAGNOSIS OF PROBLEM 

26 I AN IF STATEMENT MAY BE MORE EFFICIENT THAN A 2 OR 3 BRANCH COMPUTED GO TO STATEMENT, 

1 AN IF STATEMENT NAY BE MORE EFFICIENT THAN A 2 OR 3 BRANCH COMPUTED GO TO STATEMENT, 

99 I AN IF STATEMENT MAY BE MORE EFFICIENT THAN A 2 OR 3 BRANCH COMPUTED GO TO STATEMENT. 


SYMBOLIC REFERENCE MAP (R-1) 


ENTRY POINTS 
3 INTGRAL 


variables 

SN TYPE 

RELOCATION 

64 

80 

REAL 

MRINT 

66 

CO 

REAL 

MRINT 

70 

C2 

REAL 

MRINT 

0 

DX 

REAL 

ARRAY F.P, 

0 

NINTG 

INTEGER 

F,P* 

0 

NX 

INTEGER 

F.P# 

246 

TENP 

REAL 


71 

TIMED 

REAL 

MRINT 

0 

XO 

REAL 

ARRAY MRINT 


65 

B1 

REAL 


MRINT 

67 

Cl 

REAL 


MRINT 

0 

DT 

REAL 


F.P. 

240 

I 

INTEGER 



62 

MINT6P 

INTEGER 


MRINT 

63 

NX2 

integer 


MRINT 

0 

TIME 

REAL 


F.P. 

0 

X 

REAL 

ARRAY 

F.P. 


STATEMENT LABELS 
0 0 
40 30 

0 45 

116 70 

0 80 
161 100 
0 200 
232 1000 


34 

10 

0 

20 

INACTIVE 

0 

35 

61 

40 


101 

00 

106 

60 


0 

75 

132 

80 


146 

87 

151 

90 


0 

150 

174 

200 


211 

300 

0 

300 


Figure 

30. -Continued. 
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LOOPS 

LABEL INDEX FROH-TO 

LENGTH 

PROPERTIES 

30 

5 I 

19 20 

2B 

INSTACK 

56 

35 I 

27 30 

6B 

INSTACK 

71 

65 I 

32 36 

7B 

INSTACK 

126 

75 I 

62 66 

68 

INSTACK 

161 

85 I 

66 52 

6B 

INSTACK 

167 

150 I 

59 61 

6B 

INSTACK 

202 

250 I 

66 69 

68 

INSTACK 

221 

350 I 

76 78 

lOB 

OPT 


COHMON BLOCKS LENGTH 
NRINT SB 

STATISTICS 

PROGRAM LENGTH 2A7B 167 

CN LABELED CONNON LENGTH 72B 56 


Figure 30. -Concluded. 


SUBROUTINE RNIC 73/76 OPT«I 


FTN A.66639 78/02/15. 07.66.66 


1 


5 


10 


15 


20 


25 


SUBROUTINE RHIC 
C 

CONNON/VARDAT/ UX C 50 )» LX (50) 

CONNON/NROAT/A( 3 , 3)«B( 3» 2 ) #0X ( 39)# X (13)#U(2 ) 

1# T1*T2#T3#T6#T5#T6*T7#T8#T9#T10#T11#THR1#THR2#THR3 

2# 611#612#G21#622 

LOGICAL LX 
C 
C 

DO 20 I*l#3 
DO 10 4•lp^ 

10 A(I»J)«0. 

DO 15 J»l#2 
15 B(1#J)»0. 

20 CONTINUE 
A(1#D— 2* 

A(2#2)«*6. 

A(2#3)»-7. 

A(3>2)- 1. 

B(l#l)*2. 

B(2#2)«7. 

DO 30 I-l#39 
30 DX(I)«0. 

DT»UX(19) 

CALL 1NTGRAL(0X#X#DT#0*#13#0) 

RETURN 

END 


H0N09 1651 
H0N09 1652 
HONDO 1653 
H0N09 1656 
H0N09 1655 
H0N09 1656 
H0N09 1657 
H0N09 1658 
H0N09 1659 
H0N09 1660 
H0N09 1661 
H0N09 1662 
HON09 1663 
H0N09 1666 
H0N09 1665 
H0N09 1666 
H0N09 1667 
H0N09 1668 
H0N09 1669 
H0N09 1670 
H0N09 1671 
H0N09 1672 
HQN09 1673 
H0N09 1676 
H0N09 1675 
H0N09 1676 
H0N09 1677 


SYMBOLIC REFERENCE HAP (R«l) 

ENTRY POINTS 
1 RHIC 


VARIABLES SN TYPE RELOCATION 


0 

A 

REAL 

ARRAY 

NROAT 

72 

DT 

REAL 



123 

Gll 

REAL 


HRDAT 

125 

G21 

REAL 


HROAT 

70 

I 

INTEGER 



62 

LX 

LOGICAL 

ARRAY 

VAROAT 

121 

THR2 

REAL 


HRDAT 

105 

T1 

REAL 


HROAT 
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